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1 Introduction 



The study of non-relativistic two-body systems consisting of a heavy quark- 
antiquark (QQ) pair has a long tmdition since the discovery of the J/tp and the 
work of Appelquist and Politzcip which had shown that, owing to asymptotic 
freedom, non-relativistic quantum mechanics should apply to a good approxi- 
mation to heavy QQ systems. Thus, heavy quarkonium systems should be pre- 
dominantly Coulomb-like bound states sharing many features of positronium 
and an ideál tool studying the interplay of perturbative and non-perturbative 
aspects of quantum chromodynamics (QCD). 

The realization of this idea has proven a difficult task. Not even to the 
first approximation does the Balmer spectrum of QED give a good quantita- 
tive description of the charmonium or bottomonium spectrum. It was found 
that perturbation theory completely fails to describe the long distance part 
of the heavy quark potential. Subsequently, potential models emergedp which 
incorporated a linear rising potential at long distances in accordance to the 
idea of quark confinementEl and a Coulomb-like potential at short distances in 
accordance to perturbation theory. Potential models were highly successful in 
describing the observed spectra and transition rates, but they were not derivcd 
from first principles in QCD. As such they cannot be used for quantitative tests 
of QCD. 

On the other hand, it was found that non-perturbative effects in heavy 
quarkonia, where the quark mass m is sufficiently large, such that m 3> mv 3> 
mf ^ ^ Aqcd , v being the quark velocity and Aqcd the typical hadronization 
scale, can be described by an operátor product esijaiision in terms of local 
and gauge-invariant quark and gluon condcnsateslj'EHa Only the top-antitop 
quark systém truly belongs into this category. For bottomonia, the ground 
states and, maybe, low radiál excitations might belong to this category too, 
but mesons containing charm do certainly not. However, even for the pertur- 
bative contributions for systems in this category, it was not straightforward 
to systematically determine higher orders corrections, either due to ultraviolet 
(UV) divergences in approaches starting from the non-relativistic approxima- 
tion, or due to the amount oflechnical difficulties in the relativistic approach, 
the Bethe-Salpeter formalismD 

A new conceptual approach using an effective field t 
Caswell and LepageĚI and Bodwin, Braaten and Lepage. 
that, for the description of quarkonia, QCD can be reformulated in terms of an 
effective non-renormalizable Lagrangian built from quark and gluonic field op- 
erators describing fluctuations below m, whereas hard effects associated with 
the scale m are incorporated in the coeíHcients of the operators. The theory 



eory was proposed by 
I These authors found 
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is called "non-relativistic QCD" (NRQCD). At present, NRQCD is one of the 
standard instruments to describe quarkonium production and decay properties 
at coUider experiments. In NRQCD, production and decay of quarkonia is de- 
scribed as a sum of terms each of which is a product of an process-dependent 
hard coefficient and a universal operátor matrix element. The hard coefficient 
can be computed perturbatively, whereas the matrix element has to be de- 
termined from fits to experimental data. Since I will not touch this aspect 
of quarkonium physics in this review I refer to Refs. NRQCD also made 
lattice computations of quarkonium systems possible, since by integrating out 
hard fluctuations the required lattice siae became small enough to be manage- 
able for present day computing power.til 

For the computation of the perturbative contributions of heavy quarko- 
nium systems with m ^ mv 3> mv^ 3> Aqcd, NRQCD provided a prescrip- 
tion to deal with the UV divergences that arise in relativistic corrections to the 
non-relativistic Schródinger equationEJ However, NRQCD is not yet a satis- 
factory theory for systematic perturbative computations of heavy quarkonium 
properties, since no separation of non-relativistic fluctuations is carried out. 
The consequences are that NRQCD becomes inconsistent in dimensional reg- 
ularization, that there is no consistent power counting in v, and that it is not 
possible to sum large QCD logarithms of v. Consistent computations are only 
possible in cutoff schemes, which makes analytic QCD calculations complicated 
and practically impossible at higher orders. 

These deficiencies of NRQCD were addreaaed in two new effective 
theories callad "potential NRQCD" (pNRQCD)^! and "velocity NRQCD" 
(vNRQCD)ly The new effective theories are more complicated than NRQCD 
and based on a complete separation of the modes that fluctuate in the various 
momentům regions that exist for m 3> mv ^ mv'^ 3> Aqcd- Both theories 
are formulated in dimensional regularization. The separation is achieved by 
a strict expansion in energy and momentům components that are small in a 
given region. For the computation of an asymptotic expansion in v of Feynman 
diagrams involving noiLJ:elativistic QQ pairs this method has been called the 
"threshold expansion" Ej 

In this review I give an introduction to the perturbative aspects of heavy 
quarkonium systems with m ^ mv ^ mv^ ^ Aqcd using effective theory 
methods. The progress on the conceptual as well as on the calculational level 
in the last few years has been remarkable. This period was (at least for me) an 
exciting time, since the transition from having methods, which worked well but 
would need to be extended for new applications, to having theories has been 
achieved. Certainly, the development is not completed, but the path to go on is 
by far clearer now than leťs say five years ago. Although the application of the 
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new developments is restricted to top and some aspects of bottom physics, the 
ncw knowledge will certainly help also for a better quantitative understanding 
of Systems where Aqcd is of order mv'^ or larger. 

The review is divided into two parts. Sections cover the conceptual 
developments and Secs. |^-|Tl| the practical applications. Each section has becn 
written in a self-contained way. In Sec. || the perturbative aspects of NRQCD 
are reviewed. In particular, it is shown why NRQCD fails to be consistent when 
dimensional regularization is used. As mentioned before, the phenomenolog- 
ical application of the NRQCD factorization formalism to charmonium and 
bottomonium production and decay is not covered. In Sec. ^the method of re- 
gions and the threshold expansion are discussed. An introduction to pNRQCD 
and vNRQCD is given in Secs. ^ and ^, respectively. Since the concepts used in 
the construction of both theories have some subtle differences and since both 
theories appear not to be equivalent, I have devoted Sec. |^ to a comparison 
of pNRQCD and vNRQCD to visualize the differences. The important role 
of quark mass definitions in perturbative computations of heavy quarkonium 
properties is reviewed in Sec. |^. In Sec. ^ it is shown, how the total heavy 
quark production cross section at next-to-next-to-leading order in fixed order 
and renormalization-group-improvedperturbation theory is computed, and the 
differences between both types of perturbative expansions is discussed. Section 
^ covers the application of these computations to top quark pair production 
close to threshold at the Linear CoUider and their impact on the prospect of 
measurements of the top quark mass and other top quark properties. In this 
section I concentrate on applications of the total cross section and leave many 
other interesting aspects unmentioned. The applications of the total cross sec- 
tion calculations in bottom mass extractions from non-relativistic sum rules 
for T mesons are reviewed in Sec. Sec. ^ discusses applications of pertur- 
bative computations of the quarkonium spectrum. Conclusions and an outlook 
are given in Sec. 



2 Effective Theories I: NRQCD 

The effective theory NRQCD was proposed by Caswell and LepagJ§ and 
Bodwin, Braaten and LepageO to describe non-relativistic QQ dynamics and 
quarkonium production and decay properties. At present it is one of the major 
Instruments for the theoretical analysis of quarkonium production and decay 
data at coUider experiments. 
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2.1 Basic Ideas 

The basic idea behind the construction of NRQCD is that all hard fluctua- 
tions with frequencies of order m are integrated out. The eíTective theory is 
then built from fields for the heavy quark and the hght degrees of freedom, 
which describe non-relativistic fluctuations below the hard scale. The various 
non-relativistic scales mv, mv^ and Aqcd are not discriminated. The theory is 
regulated by a momentům space cutoff A that is of order m. Dimensional regu- 
larizationieads to inconsistencies in NRQCD if the power counting of Ref. |^] is 
employedll3 Therefore, a consistent renormahzation group scahng of NRQCD 
operators is difficult to conatruct. From NRQCD an effective Schrodinger for- 
mahsm can be constructedcl that contains instantaneous QQ potentials. 

2.2 Cutojf Regularization 

The NRQED Lagrangian readsii 

D2 cr B 



-Cnrqcd = -0' < iDo +7; \- Cl— ^ + C2 g- 



2m 



(D-E-E-D) cr.(DxE-ExD) , 

+^^^ 8^^^^ + 8^^ + ^^+(^^^) 



4to 
d2 9' 1 
'3m4 2 
1 

"4 



(^Va2X*) ix' <J2CTÍ-^Dyi;) + h.c. 
C^^G^, , (1) 



where and x are two-component quark and antiquark PauU spinors; íDq — 
ido — gAo and žD = iV + gA are the time and space components of the 
gauge covariant derivative Í3^, = and = ^e^^^G^^ are the chromo- 
electric and chromo-magnetic components of the gluon field strength tensor, 
igG'^'^ = [D^,D'^]. Operators involving hght quark are not displayed, color 
indices are suppressed, and m is the heavy quark pole mass. At Born Icvcl the 
effective Lagrangian is obtained by simply expanding the vertices of fuU QCD 
in the non-relativistic Umit. Bcyond Born level the short-distance coefíicients 
Ci and di encode the effects from momenta of order m, which are integrated out. 
They are determined from matching NRQCD amplitudes to QCD amplitudes. 

From the above Lagrangian one may in principle derive exphcit Feynman 
rules and carry out bound statě computations. Questions concerning the power 
counting are relevant only in so far as one needs to include a sufficient number 
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of operators for a computation at a certain order. In particular, according to 
the power counting rules in Ref. at leading order one has to include the 
quark kinetic term ^ because Dq ~ ^ ~ mv'^, and the correct leading order 
quark propagátor reads 

i 

kn - ^ + i5 

This is the approach of lattice calculations, where the dynamics is computed 
in a non-perturbative manner. In a cutoff scheme a diserimination between 
different non-relativistic fluctuations with soft (~ mv) or ultrasoft (~ mv'^) 
frequencies is not mandatory conceptually, since the cutoff automatically in- 
cludes only the relevant fluctuations. For analytic bound statě calculations, 
however, the present formalism is still too cumbersome to be used in practical 
calculations because the contributions from the Coulomb potential exchange 
need to be summed to all orders, whilc sublcading tcrms can be treated as 
perturbations. The NRQCD Feynman rules do not yet provide a systematic 
distinction of the various contributions. Therefore, to carry out analytic bound 
statě calculations, the NRQCD framework needs to be rewritten in terms of a 
Schrodinger theory which contains potentials and radiation effects and which 
can be used to carry out perturbation theory for bound statesul In fact all ana- 
lytic NRQCD bound statě calculations were carried out in such a Schrodinger 
theory. 

The conceptually most,developed approach to such a Schrodinger theory 
was suggested by LabelleyJ The most efficicnt gauge for this program is the 
Coulomb gauge. In the Coulomb gauge, the longitudinal gluon (the time com- 
ponent of the vector potential) has an energy- independent propagátor, 
This means that the interaction associated with the exchange of a longitudinal 
gluon corresponds to an instantaneous potential. The power counting of di- 
agrams containing instantaneous potentials is particularly simple, because an 
instantaneous propagátor has no particle pole (and no i(5-prescription) and the 
scale of k is just set by the average three-momentum of the quarks ~ mv in the 
bound statě. So one readily obtains the well known result that the Coulomb 
potential is the leading order interaction. On the other hand, the transverse 
gluon (the spatial component of the vector potential) has an energy-dependent 
propagátor of the form 



(kl - k2 + i5) \ k2 



In this case, the propagátor has a particle pole, and fco and k can be of order 
mv (soft) or mw^ (ultrasoft). Also the regime with fco ~ mw^ ^ k ~ mv is 
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described by the propagátor. This has the consequence that NRQCD diagrams 
containing transverse gluons involvc contributions from soft and ultrasoft scales 
and, therefore, do not contribute to a unique order in w. In a cutoff regulariza- 
tion scheme this is a priori not problematic for the consistency of the theory, 
but it makes analytic calculations unnecessarily comphcated, because there is 
no transparent power counting. Therefore one does not gain too much in sepa- 
rating only the hard fluctuations. Labelle suggested to also separate exphcitly 
the soft and ultrasoft fluctuations. A diagram containing a quark-antiquark 
pair and a transverse gluon contains a loop integrál of the generic form 

(2) 

{ko + § ~ ^ +^S][ko - § + ^ -^S][k|-k^ +^S]' 

and p ~ mv stand for some external center-of-mass energies 
and momenta. Carrying out the fco-integration by contours, which picks up 
the gluon and quark poles, one obtains terms of the form 

(3) 



where I have dropped a factor 2in/{E — ^^^^ — h iS) for the quark-antiquark 
propagation. From this, one can see that different contributions arise depend- 
ing on whether the scale of k is set by p ~ mv or by _B ~ p^/m ~ mv'^. The 
effects associated with the latter scale are the retardaidon effects, which are for 
example responsible for the Lamb shift in hydrogentžl The contributions from 
both momentům regions will not^cpntribute to the same order in v, i.e. the 
power counting is broken. LabcUetZl suggested to reformulate NRQCD in such 
a way that the contributions associated with the different scales are coming 
from separate diagrams. This is achieved by Taylor-expanding the NRQCD 
diagrams containing Eq. (|^) around k mv and around k ~ mv'^. The pre- 
scription anticipated one of the crucial ingredients of the threshold expansion 
(Sec. |) and the effective theory vNRQCD (Sec. |). It is a generalization of 
the multipole expansion of vertices in p^/m and by now generally called "mul- 
tipole expansion" in the literatuře. One finds that the lowest order term of the 
expansion around k ~ mv gives a contribution of order 



j3,, J i^^)^ 

|kp (mvY 



d3k_i_^^^^^. (4) 



It originates from the region fco ~ ^ ~ "iw^ in Eq. (g) . Higher order terms in 
the expansion also contain contributions from the region fco ~ k ~ mv. The 
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lowcst order term of the expansion around k ^ mv'^ gives 

/ d-^k 2 7 ^ - mv^ . 5 

i |k| [f - - |k|] (m«2)2 

This contributions comes from the region kg ^ k ^ mv^ . The resuhs from 
Eqs. (^) and (||) show that at leading order the transverse photon propagátor 
reduces to an instantaneous potential, which simply corresponds to approxi- 
mating the transverse photon propagátor l/(fcQ — k^) by — Due to the 
w-suppression of the couphngs of a transverse giuon to heavy quarks, the po- 
tentials from the transverse gluons are of order with respect to the Coulomb 
potential. The prescription of Labelle allows the construction of a Schrodinger 
theory where the leading order binding of the QQ pair is deseribed by the well 
known non-relativistic Schrodinger equation and where higher order correc- 
tions such as w-suppressed potentials and retardation effects can be computed 
separately as perturbations. 

All analytic bound statě calculations that were based on NRQED or 
NRQCD have essentially used the prescription proposed by Labelle. No ex- 
plicit discussion of QED calculations shall be given at this pláce. For a number 
of originál results in the NRQED framework on the muonium hyperfine split- 
ting and on positronium, see e.g. Refs. ]l^ , pO| , p]] | . A major appHcation of this 
formalism in QCD was the determination of the total cross section of QQ pairs 
close to threshold at next-to-next-to-leading order (NNLO) in the fixed order 
non-relativistic expansion, which is discussed in Secs. [7|-pX|. 

Although a cutoff regularization is in principle feasible, it has a number of 
disadvantages, particularly in QCD. In analytic QCD bound statě calculations 
the cutoff breaks gauge invariance at intermediate steps. This problém does not 
occur on the lattice, because there NRQCD can be defined in a manifest gauge- 
invariant way. In addition, in analytic QCD bound statě calculations a cutoff 
introduces also subtle scheme-dependences at intermediate steps, because the 
results for diagrams in the effective theory depend on the routing of loop 
momenta (see e.g. Ref. l22| ). 

A problém with both analytic and lattice computations is that the power 
counting is nevěr manifest with a cutoff prescription, regardless whether all 
relevant momentům regions are separated or not. (In the latter case, however, 
the effects of power counting breaking are more severe.) The reason is that 
the cutoff A ~ m arises as an additional scale in the diagrams of the effective 
theory. This introduces a sensitivity to the hard scale m, when there are linear, 
quadratic, etc. UV divergences. This means that an operátor that contributes 
for the first time, leťs say, at NNLO according to the non-relativistic power 
counting can lead to lower order contributions, if A ~ m. Such terms are 
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Figuře 1: The two-loop NRQCD vector current correlator with the exchange of one Coulomb 
gluon without relativistic corrections (a) and with the kinetic energy correction indicated by 

the crosses (b). 



compensated by NLO matching corrections of coefficients in the effective La- 
grangian. In practice this means that whenever additional higher dimensional 
operators of the effective theory are included to increase the order of approxi- 
mation, the matching conditions of all coefficients, i.e. also of the coefficients of 
lower dimensional operators, obtain additional lower order corrections. This 
is a feature well known in NRQCD lattice computations, see also Ref. [ p3| . 
To illustrate this issue consider the two-loop NRQCD vector current correla- 
tor at LO in Fig. |l|a, which contains the exchange of a longitudinal Coulomb 
gluon. The power counting in Labelle's scheme indicates that the dominant 
contributions to the absorptive part of the diagram is of order a^w^, 



Im[Fig. la] 



3 Nc Cf as 



TT 

T 



4^ 
A 



O 



(6) 



where ±p is the three-momentum of the produced quarks, and Nc = 3, Cp = 
4/3. Now consider the kinetic energy corrections in the same two-loop diagram 
in Fig. |l|b caused by the terms in the NRQCD Lagrangian. According to 
Labelle's scheme the dominant contribution to the absorptive part is of order 
a,v^, but the actual result reads 



Im[Fig. Ib] = 



A|P| 



2m2 



O 



(7) 



and contributes at order asV for A ~ m because of a linear UV divergence. 



2.3 Dimensional Regularization 

It has become common practice in analytic QCD computations to use dimen- 
sional regularization to regulate UV and IR divergences. Dimensional regu- 
larization maintains gauge invariance and provides a comfortable framework 
to determine the anomalous dimensions of the operators in effective theories. 
In addition, power counting breaking effects as mentioned at the end of the 
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previous section do not arise, because power divergences are automatically set 
to zero. However, it turns out that the NRQCD Lagrangian cannot be de- 
fined in dimensional regularization, if the non-relativistic power counting rule 
Dq ~ ^ is maintained. Manohar0 pointed out that the resulting efFective 
theory contains spurious particle poles at the hard scale, which lead to the 
breakdown of the power counting and make a consistent matching to QCD or 
even the determination of anomalous dimensions impossible. As an example, 
consider the one-loop selfenergy of a quark due to a gluon. Taking the external 
quark in its rest frame the selfenergy has the generic forni 



d^k 



(8) 



(fc2_k2+í<5)(fco-|l+zá) 

Carrying out the fcp integration by contours one finds a term proportional to 



2m) 



(9) 



which is non-zero and dominated by hard momenta fc ~ m. The same feature 
arises at higher loop orders and for QQ scattering diagrams. Thus NRQCD 
does not properly describe non-relativistic degrees of freedom, if it is defined 
in dimensional regularization and if the non-relativistic power counting Dq ^ 
^ ^ mv^ is imposed. 

Manohar pointed out that the problém can be resolved for the single quark 
sector (i.e. all operators bilinear in the quark field and with arbitrary number 
of light fields) by treating the kinetic energy term as a perturbation. In 
this case the selfenergy has the form 



d^k r 1 



{kl - + i6) \ (fco + i5) 2rn(fco -f- iSf 



(10) 



and vanishes at any order in the expansion. A sensitivity to the hard scale 
does not arise. It is important to realize that the integrals in Eqs. m\ and 
( |lo| ) are not equivalent, even if the terms in the l/m expansion in Eq. (^) are 
summed up. This is because in dimensional regularization, expansion and in- 
tegration do not commute, in contrast to a cutoff scheme. Therefore the single 
quark sector of NRQCD can onlv be consistently defined if the 1 /m counting 
of heavy quark effective theorytil (HQET) is adopted. A comparable prescrip- 
tion for the four-quark sector of NRQCD was not found. However, Grinstein 
and RothsteinĚa pointed out that Labelle's multi-pole expansion can be under- 
stood as a generalization of Manohar's prescription (see also Ref. [|6[). Thus, 
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Manohar's prescription foUows from the fact that in the multipole expansion 
the combination of propagators in Eq. (j^) does simply not exist. When keep- 
ing the form l/(fcQ — k^) for the gluon propagátor one has ^ li m and 
the quark propagátor has to be expanded in l/m. On the other hand, when 
keeping the form l/(fco — ^) for the quark propagátor one has fco ^ k <C m 
and one has to expand the gluon propagátor in fco. The latter contribution is 
a scaleless integrál and vanishes to all orders in the expansion. 



3 Method of Regions and Threshold Expansion 

In bound statě calculations one frequently needs a small-velocity expansion of 
Feynman diagrams involving a heavy quark pair close to threshold. The thresh- 
old expansion is a prescription developed by Beneke and SmirnovO to carry 
out such an expansion for loop diagrams in dimensional regularization. The 
prescription of the threshold expansion is based on the more generál "method 
of regions" and fuUy appreciated for the first time the meaning of the multipole 
expansion in dimensional regularization. 



3.1 Basic Idea 

The problems in constructing an asymptotic expansion and the basic idea of 
the method of regions can be illustrated with a simple example. Consider the 
one-dimensional integrál 

-t-oo 

f |arctan(fc)| 

— oo 

which one might take as a simplified version of a diagram with propagator- 
like factors of 1/(A;^ + a^) and a non-trivial numcrator structure. The task is 
to expand /(a) in a ^ 1. Here it is impossible to naively expand / before 
integration, because each term is "IR divergent" at fc = 0. Altogether there 
are two relevant integration regimes one has to consider, the "hard" regimes 
where fc ~ 1 and the "soft" regime where fc ~ a. In a cutoff scheme one deals 
with this situation by introducing a cutoff A with a ^ A ^ 1. The cutoff 
separates the fuU integration range into a "soft" regime |fc| < A and a "hard" 
regime |fc| > A. In the "soft" regime we have fc ~ a <C 1, and we can expand 
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the numerator. In the "hard" regime we have a and we can expand in a, 



Carrying oiit the Taylor cxpansions bcforc the intcgration considcrably sim- 
plified the computation. but it is not rcaUy necessary for the method to work 
because expansion and integration commute in the presence of the cutofF A. 
Adding together the "hard" and the "soft" contributions we readily find the 
correct asymptotic expansion of the integrál /, 



The A-dependent divergcnt tcrms cancel in the sum. The logarithmic divergent 
terms are speciál because their existence is still visible in the finál sum through 
the non-analytical In a term. Thus the "large" logarithmic In o terms originate 
from logarithmic divergences at the borders of the regions. If we would consider 
a similar situation for the computation of a complicated multi-loop Feynman 
diagram with a non-trivial structure of relevant regions in QCD, the method 
describcd above would be in principlc fcasiblc, but quite uncomfortable and 
cumbersome due to the existence of additional cutoff scales. In QCD this 
also leads to the breakdown of gauge invariance in intermediate steps of the 
computation. 

A more economic (but also less physical) regularization method is to use 
an analytic regularization method. In QCD computations the most common 
ehoice is to \mc dimcnsional regularization in the MS scheme, where the four- 
dimensional integration measure is continued to D = 4 — 2e dimensions, and e 
is an arbitrarily small complex parameter. In our simple example we continue 
the integration measure to .Ď = 1 — 2e dimensions, 






where 



n{D) = 



r(f) 



is the f-dimensional angular integrál and 
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7b = 0.577216 . . .. This is just onc spccific regularization mcthod of an infinite 
number of possibilities to regulate the integrál analytically. The separation of 
the "hard" and the "soft" regions is not implemented through a restriction to 
integration regions, but by strictly Taylor expanding out the hierarchies in the 
regions. For the contributions from the "hard" and the "soft" regions one now 
finds 



The sum of the terms gives again the asymptotic expansion of / and the fi- 
dependent terms cancel. Linearly divergent contributions do not arise at all, 
whereas the large logarithmic terms In a are associated with the 1 /e singu- 
larities and In^ terms that arise from divergences at the boundaries of the 
regions. The method also works if the originál integrál is divergent itself for 
e 0. In this oase the outcome of the method agrees with the result of the fuU 
calculation in the limit e — í- O supplemented by an expansion in a. Choices for 
the analytic regularization schcmc that are different from the one used in the 
simple example above lead to different intermediate results for the "hard" and 
the "soft" contributions, but to the same finál result, if the originál integrál is 
finite. The method breaks down if an expansion it not carried out in a strict 
way. In particular, nothing is gained keeping some contributions unexpanded 
with the intention to sum up certain terms. 

It is an important feature of this method that each term in the expansion 
only contributes to a single power of a and that the order of a to which each 
term contributes can be easily determined before integration by considering 
the a-scaling of k in the two regions. In the "hard" regime we have k ^ 1 
and dfc ~ 1, and we just expand in powers of a up to the required order in a. 
In the "soft" regime we have k ^ a and dk ~ a, and 1/(A:^ + a?') ^ for 
the propagator-like term. Thus the leading term in the expansion in the soft 
region is of order a~^, etc. 

3.2 Threshold Expansion and Regions 

The prescription for the threshold expansion goes along the lineš of the previ- 
ous simple example. In Ref. fl^ ] four different loop momentům regions were 
identified to be relevant by an analysis of Feynman diagrams in fuU QCD 





,3 



■soft' 
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involving non-relativistic QQ pairs close to threshold, 
hard : (fc'^,k) ^ {m,m) , 
soft : ífc*^, k) ^ (mw, TOw) , 

(11) 

potential : (fc^^jk) ^ {mv^,mv) , 

ultrasoft : (fc°,k) ^ {mv'^ , mv'^) . 

The energy and thrcc-momcntum componcnts of thc loop momentům in the 
potential region scale differently with v, because Lorentz- invariance is broken 
for QQ pairs close to threshold. The structure of the regions is most transpar- 
ent if the QQ pair is in the center-of-mass frame. 

The small velocity expansion of a fuU QCD diagram, which involves a non- 
relativistic QQ pair, is obtained by writing it as a sum of terms that arise from 
dividing each loop into all possible regions. The separation of the regions is 
achieved by strictly expanding loop momenta according to the hierarchy in the 
various regions. All terms from all regions that arise from this proceduře are 
integrated over the fuU D-dimensional space using dimensional regularization 
in leťs say the MS scheme. All scaleless integrals have to be set to zero. The 
explicit expansion of energies and three-momenta that are small is essential to 
make the method work, because the expansion is the (only) instrument that 
separates the regions. 

It is one of the amazing features of dimensional regularization that this 
heuristic prescription, which involves subtle cancellations of IR and UV diver- 
gences in neighboring regions, appears to work. This "method of regions" is 
quite generál and has been used earlier for other problems, such as the large 
mass or the large energy expansioncZI Also extended and modified versions of 
the threshold expansion for different kinematic situations and particle content 
have been devisedt30 A mathematical and more rigorously defined formu- 
lation of the method in terms of so called R-operations has been developed 
earlier by Smirno\E£l (see also Ref. [0). Some recent formal considerations 
pointing out potential problems can be found in Ref. Up to now the 

method has not been proven mathematically to work in generál for arbitrary 
diagrams with on-shell external particlespi but no counter example has been 
found in cases where the expansion of the complete result was available. This 
should be kept in mind, since for many multi-loop calculations an asymptotic 
expansion using the method of regions seems to be the only way to tackle the 
problém in the firat pláce. A mathematical proof, however, exists for off-shell 
external particlesEj It should also be noted that it is in principle not excluded 
that regions other than in Eq. (^) could become relevant in some cases. In 
such a case the new regions would simply have to be included without affecting 
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the method itself. Interesting examples (but not dealing with non-rclativistic 
QQ pairs), where the number of regions incrcascs indefinitely at higher loop 
orders, have been reported in Ref. [|32|| . 

In actual calculations, the threshold expansion has been apphed correctly, 
if each integrál only contributes to a single power in ti. The order of v to which 
a certain term in the threshold expansion contributes, can be determined by 
power counting rules that are obtained from the w-scaling of the regions in 
Eq. (|ll|) similar to the simple example discussed in the previous section. The 
feature of the power counting makes the threshold expansion method quite 
economic, because all terms needed for a certain order in v can be idcntified 
unambiguously before doing any integrations. For example, the integration 
measure ď^k in the hard, soft, potential and ultrasoft regions scale u**, 
and v^, respectively. 

The threshold expansion had a significant influence on the development 
of non-relativistic effective theories. Its features, such as the existence of the 
diffcrent separable momentům regions, the existence of power counting rules 
and the requirement to "multipole-expand" in small momenta in dimensional 
regularization also play an important role in effective theories. The present 
use of language, for example calling a gluon propagating in the soft region 
a "soft gluon" , or a quark propagating in the potential region a "potential 
quark", etc, was finally established with the threshold expansion. Since the 
time it has been devised, a number of calculations, which were important 
for matching calculations in effective theories for non-relativistic QQ pairs, 
have been carried out using the threshold expansion, see e.g. Refs. [ |3^ , p^ , ^ . 
However, the threshold expansion is a prescription to obtain the asymptotic 
expansion of Fcynman diagrams and not an effective theory itself, because 
it includes on-shell as well as ofF-shell modes. (An attempt to formulate the 
threshold expansion in form of an effective Lagrangian containing fields for 
the off-shell soft heavy quarks has been given in Ref. [^6|.) In addition, the 
concept of renormalization is not implemented into the threshold expansion. 
Therefore the summation of large logarithmic terms is impossible. It should 
also be noted that the definition of the regions in Eq. (^l]) is not necessarily 
cquivalcnt to the corrcsponding regions in effective theories duc to the freedom 
in the choicc of the operátor basis. For example, the contributions obtained 
from the hard momentům regime in the threshold expansion are in generál 
not equivalent ta the hard matching conditions of operators obtained in an 
effective theory.ĚZl However, an effective theory has to be able to reproduce the 
sum of contributions from the different regions of Eq. ([ll|). 
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3.3 Example: One-Loop Vertex Diagram 

As a simple example for the application of the threshold expansion let us 
consider the one-loop vertex diagram in Fig. H, which describes production 



Figuře 2: One-loop scalar vertex diagram with external massive lines close to threshold. 
Solid lines are massive, and the dashed line is massless. AU lines represent scalars. 

of two massive scalars (solid lines) with momenta pi/2 = (po,±p) and the 
subsequent exchange of a massless scalar (dashed line). The presentation given 
here is slightly different to the one of Ref . ||lj] with respect to the interpretation 
of the soft and ultrasoft contributions. Close to threshold we have mE = 
Pq — ~ p2 ^ (mvY . We assume that the produced massive scalars are on- 
shell, i.e. we have the relation mE — p^. The on-shell relation is, however, not 
essential for the threshold expansion as long as the w-counting of the external 
energies and three-momenta is well-defined. The fuU loop integrál reads {D = 
4- 2e) 

d°k 1 



{2t:)D [(fc +pi)2 - „,2]p _p2)2 - m2][fc2] 
d^k 1 



(27r)^ [(fco + Po)2 - (k+p)2 - m2] [(fco - po? - (k+p)^ - m^] [kl - k^] ' 

(12) 

The renormalization scale ^ is set to 1 for simplicity and an conventional over- 
all factor (47r)^' is undcrstood. The usual +i5 prescription in all propagators 
is implied. In the foUowing we determinc the respective leading order contri- 
butions from the different regions. 

In the hard regime we have to expand the propagators in p, p' and mE. 
To leading order we obtain, 

d^k 1 



(27r)^ [^2 + 2koPo + iS] [A:2 - 2knpo + iS] [k"^ + iS] 

i r(e) 



32 7r2(p2)i+e i + 2e 



(13) 
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The result is of order m^^, which can be obtained before integration from the 
scaling ď^k ~ to'* and fco ~ k ~ po ~ "i. 

In the potential regime we have to expand thc propagators in fcg. To 
leading order we obtain, 

d^k 1 



(27r)^ [+2fcoPo-(k+p)2+TO£; + z(5][-2fcoPo-(k+p)2+mí; + í(5][-k2] ' 

(14) 



For on-shell external particles the result reads 

(15) 



r(í + é; 



32po{-mE ~iS)2+'^ 2 7r2e 

It is of order m~'^v~^, which can be obtained before integration from the scahng 
ď^k ~ TO^w^, fco ~ mv'^ and k p ~ mv. The massless propagátor does not 
have an iS prescription since it does not depend on fco and cannot develop a 
particle pole. 

In the soft regime the fcopo terms in the massive propagators are dominant. 
The leading term is of order (mv)^ (rn?v)~'^ {rn^v'^)~^ = and reads 

ď^k 1 

{2Tr)D [+2fcopo + iS] [-2fcopo + i5] [fcg - + iS] ' ^ ' 

The integrál is scaleless and proportional to {-^ ~ ^^'^ there is no contri- 
bution from the soft region. However, the integrál contains a pinch-singularity 
at fco = O and is in fact ill-defined. However, the poles at fco = ±iS in the com- 
plex fco-plane are unphysical because the massive lineš cannot become on-shell 
in the soft region. Thus, the integrál can be defined by removing the i(5's from 
the massive propagators with the prescription that the fco is carried out first 
by contours. 

Finally, in the ultrasoft regime we have to expand the massive propagators 
in fc2 and k • p. The leading term is of order {mv'^)'^{m?v'^)~'^{m?v'^)~^ 
and reads 

d^k 1 



9 
TO 



{2n)D [+2koPQ + mE-p^+ iS] [-2fcoPo + mE - p2 + iS] [k^ -k^+ iS] ' 

(17) 



Interestingly, the ultrasoft contribution is the same as the soft one in Eq. (16) 



for on-shell external momenta. However, soft and ultrasoft contributions are 
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not equivalent, because for off-shcU cxtcrnal momcnta thc soft contribution 
in Eq. (|l^) has the pinch-singularity, the ultrasoft contribution, on the other 
hand, does not. In the ofF-sheU case Eq. gives an UV divergence oc 
In the on-sheU case there is an additional IR divergence and the overall contri- 
bution is zero. The complete one-loop resuh for Eq. ( p^ for on-shell momenta 
can be determined analyticaUy and reads 

* í T? -A^-i-e IP f ^ ^ , ^ Po ^ r(l + e) 
;[-mE-iS) 2F1 -,1 + e, -, 



167r2' ' \2' ' 2' mE + i6 J 2e 

where 2F1 is the hypergeometric function. An cxpansion in E yields the con- 
tributions determined with the threshold cxpansion. 



4 EfFective Theories II: Potential NRQCD 

The effective theory pNRQCD, first proposed by Pineda and Soto^ has been 
devised for the description of non-rclativistic QQ systems where m S> mv 3> 
mv'^ ^ Aqcd as well as for cases where mv'^ is of order or small than Aqcd- 
In particular, pNRQCD is proposed to describe dynamic QQ pairs, which 
are relevant for physical apphcation, as well as static QQ systems, which, for 
example, can be studied non-perturbatively on the lattice. In pNRQCD the 
separation of the momentům regions is achieved in two steps and dimensional 
regularization is used. In this presentation I will mainly discuss the case m 3> 
mv 3> mv'^ 3> Aqcd , because it has the most solid level of understanding. 

Basic Idea 

The basic idea behind the construction of pNRQCD is that the quark and gluon 
field components that are off-shell in the hard, soft, potential and ultrasoft 
momentům regions can be integrated out in two steps. To be more precise, 
pNRQCD is the theory that results from this two-step proceduře. Starting from 
fuU QCD, first the quark and gluonic off-shell degrees of freedom in the hard 
region are integrated out at the scale ^ — m. The resulting theory for < m is 
equal to NRQCD (Sec. ||) for bilinear quark terms and, for simplicity, also called 
NRQCD in this context. This theory is then scaled down to the soft scale /i = 
mv, where in a second step the quark and gluonic degrees of freedom in the soft 
regime and the gluonic degrees of freedom in the potential regime are integrated 
out. The theory that results for /i < mv is called "potential" NRQCD, because 
the four-quark interactions that are generated through the matching proceduře 
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are the potentials in thc Schrodinger perturbation theory that can be derived 
from pNRQCD. The theory pNRQCD is then scaled down to fi = mv^ , where 
matrix elements are calculated. At the scale /x = mv^ matrix elements should 
be free of any large logarithmic terms, and all logarithms should be resummed 
into Wilson coefíicients. GraphicaUy the scheme is as foUows: 

-Cqcd (m > m) 

-Cnrqcd {m> fi> mv) (18) 

-CpNRQCD {mv > ^ > toí;2) 



4-.2 Matching, Running and Power Counting in NRQCD 

The intermediate theory, NRQCD, is constructed in close analogy to the ef- 
fective Jjheory proposed by Caswell and LepageH and Bodwin, Braaten and 
LepageÉI discussed in Sec. ^ This nieans in particular that the gluon field 
describes all fluctuations of frequencies below to. To avoid a sensitivity of 



NRQCD loop diagrams to the hard scale (Sec. 2.3), the non-relativistic power 
counting is abandoned and a strict 1/to counting is carried out everywhere. 
Thus the theory is defined in a statijC-jexpansion, and the bilinear quark sector 
of NRQCD is equivalent to HQETI13 The 1/to counting is not only applied 
to the bilinear quark sector but, in particular, also to interactions between 
quarks and antiquarks such as four quark operators. Therefore, the inter- 
mediate theory in the scheme of Eq. ( p^ ) is notj^quivalent to the NRQCD 
theory originally proposecLby Caswell and Lepaga3 (even when supplemented 
by Manohar's prescriptionll3 for the single quark sector) but rather a specific 
extension of it. Nevertheless, for simplicity we will call the intermediate theory 
of Eq. (|l|) also "NRQCD". Up to order I/to^ the La grangian has the form 

D2 o- B 



t r D2 

-Cnrqcd = W \ + ck- 1- C4- — - + CFg 

[ 2m 8ni-^ 



2to 



(D-E-E-D) CT -(DxE-ExD) , 
+cng'- '- + .csg '- + ... + x) 

TO^ m"^ 

+%V^T-VX^T-X + ^^^TV^x^TVx 
TO^ 
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^G^^^G^, , (19) 



1 

where ijj is the Pauli spinor that annihilates the massive quark, x is the Pauli 
spinor that creates the massive antiquark, and 



iDo = ida - cjAq , ÍD = + gA , igG'"' ^ [Z?^, D''] . 

The convention for the four quark operators has been taken from Refs. [ ^8|j39| 
and 771 is the pole mass of the heavy quark. Operators involving hght quarks 
are not displayed. In the foUowing I take the convention that, when I mention 
gluons, I refer to hght quarks as weU. It is important that the overah 1/777- 
counting for QQ scattering treats the kinetic energy term ^ as a perturbation 
and, at the same time, removes the potential region from the theory, because 
the hierarchy 

777 Z777 

is enforced. This can be seen from the fact that in dimensional regularization 
an expansion of the quark propagátor in l/7ř7, 

(20) 



kO- }sL+iS kO + iS (fco + iSy 2777 



sets ^ |k| of order or smaller than vfiv. It is claimed in Refs. [|3|j3||^] that 
the removal of the potential region does not affect the matching and running, 
or any physical prediction for dynamical quarks made with pNRQCD. In par- 
ticular, it is stated that the potential region is still contained in the theory, but 
hidden in a subtle way (see Sec. ^^). For the bilinear quark sector the matching 
proceduře is carried out using the prescription of Manoharllj For amplitudes 
describing quark-antiquark interactions and also for currents describing QQ 
production the matching to QCD is carried out on-shell and at threshold, i.e. 
the four-momenta of the external quarks in the full QCD amplitudes need to 
be set to fc = (tti, 0). The idea is that in this kinematic situation all NRQCD 
loop computations give exactly zero, because the integrals are scaleless. For 
the loop computations in full QCD only the hard contributions should be left 
over. At the one-loop level this can be seen from the vertex computation dis- 
cussed in Sec. ^ using the threshold expansion. Setting the external momenta 



for the massive lineš in the hard loop, Eq. (|13[), to pi,2 = (777, 0) we obtain 
dPk 



{2t:)D + 2fco777 + i5][B - 2fco777 + + iS\ 

r(e) 

32 7r2 (7772)1+' l + 2e ' 



(21) 
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On the other hand, thc potcntial, soft and ultrasoft loop contributions are 
scaleless integrals that vanish. The 1 /m expansion is not carricd out for those 
integrals since they are part of the complete computation. This prescription 
should hold for any number of loops, because any connected potential, soft and 
ultrasoft loop integrál in the threshold expansion for QQ scattering diagrams is 
zero in this particular kinematic situation. To obtain the matching conditions 
for operators that contain covariant derivatives a suitable number of derivatives 
with respect to the external momentům should be applied before setting the 
external quark momenta to (m, 0). Matching appears to be impossible for any 
other choice of external momenta due to the l/m expansion in NRQCD. 

The matching conditions sá ii — m for the two-quark operators are equiva- 
lent to the results in HQET, and for the four-quark operators they were derived 
in Ref. Il^] (see also Ref. [|3^]). The RG equations for the coefHcients of the 
two-quark sector are again equivalent to those of HQETull The RG equations 
at one-loop for the four-quark operators and their solution for thc cocfficients 
can be found in Ref. | ]39[ |. 

4-3 Matching, Running and Power Counting in pNRQCD 

The finál theory, pNRQCD, is formulated in terms of a quark-antiquark pair 
and a single gluon field that-dcacribes ultrasoft fluctuations. The Lagrangian 
in configuration space readsEJM 

/:pNRQCD = Tr <^ SM í9o - Cfc5_ + C4-P^ _ - ^ _ li^ + . . . S 

-^gV^Tr {O^r • E S + S^r • E 0} + S-yTr {O^r • E O + O^Or • E} 

-^G^.G^% (22) 

where "S"' describes a QQ color singlet and "O" a QQ color octet, and where 
the Ci and Ví are Wilson coefficients. The terms displayed in Eq. ( p2| ) are 
sufficient for a next-to-next-to-leading logarithmic (NNLL) description of the 
QQ singlet dynamics. The operators involving two QQ fields and the ultrasoft 
gluon E-field describe the leading order interaction of an ultrasoft gluon with 
a QQ pair. In the matching proceduře this interaction is the gauge-invariant 
combination of NRQCD diagrams with (ultrasoft) gluon interactions with each 
of the quarks and the (potential) gluon that is exchanged between the quark 
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msmsL> 



Figuře 3: NRQCD diagrams that contribute to the r ■ E interactions of ultrasoft gluons 

with the QQ pair in pNRQCD. 



pair (see Fig. The bilinear QQ operators involving the V^s describe in- 
teractions between QQ pairs, the potentials. This is the origin of the "p" in 
pNRQCD. Because pNRQCD distinguishes between ultrasoft and potential en- 
ergy and momenta, the multipole expansion is carried out, i.e. the dependence 
of the gluon field on the relative distance between the quarks is expanded out 
and all gluon ficlds in Eq. (|2|) are only functions of the center-of-mass coordi- 
nate and time. The explicit r-dcpendence of the operators describing ultrasoft 
gluon radiation arises from the multipole expansion. Due to the multipole ex- 
pansion the potentials only depend on the relative distance between the quarks 
and not on time, i.e. they are instant aneous, but non-local in space. Up to 
order 1/m^ the singlet potentials read {Cf — 4/3, Ca — 3) 



yiO) = 


— Of , 




r 




CfCaDÝ^ 


m 






cfd^'1 r 1 


m 


2m'^ \ r 











y(0) ^ ['^„c 



2 r 



(23) 



2m2 



where S'i2(r) = 3r • cri r • (T2/r^ — ai ■ ct^-, S = cri/2 + (T2/2, and the í?'s are 
Wilson coefficients. The term vfi^^ is the static singlet/octet potential. The 
pNRQCD Lagrangian can in principle also be written in momentům space 
and/or in terms of quark and antiquark fields. However, such a formulation 
has not yet been given explicitly. Nevertheless, in cases where the transfer 
to the momentům picture is straightforward, I will use the momentům space 
representation at some instances in the foUowing presentation. 

Because NRQCD is treated in a l/m expansion also the pNRQCD ampli- 
tudes are power-counted in l/m during the matching computationst^cJ This 
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is necessary because the kinetic energy term provides an IR cut-ofF for the 
Coulomb singularity (in the non-relativistic power counting), and because IR 
regularization in fuU and efFective theory has to be equivalent to obtain the 
correct matching conditions. Thus NRQCD and pNRQCD amphtudes are 
matched order by order in í/m and in the muhipole expansion. In Ref. |40 
Brambilla et al. carried out the matching calculation with Wilson loops that 
describe a static quark-antiquark pair with fixed relative distance r interacting 
for a very large time T ^ -. In terms of quarks and gluons and in momentům 
space this matching proceduře corresponds to comparing QQ scattering amph- 
tudes, where the incoming and outgoing quarks have off-shell non-relativistic 
momenta (O, ±p) and (O, ±p'), respectively. In this kinematic situation all pN- 
RQCD loop diagrams (in the l/m expansion) are scaleless integrals because 
the quark propagátor behaves like a static propagátor, 



l/m expansion 



(k- 



(k+p)" 



-I- iS 



2m 



(24) 



and because radiation of gluons as well as potential interactions depend poly- 
nomially on the external momenta due to the multipole expansion. Thus the 
NRQCD scattering amplitudes give directly the pNRQCD matching coefS- 
cients. As an example, consider the matching calculation for the singlet static 
potential in the limit m — oo. The pNRQCD scattering amplitudě to all orders 
in as at the scale ^ = reads 



AniCpCív^ (Ms 
(P-P')' 



(25) 



The corresponding NRQCD amplitudes are built from static quark propaga- 
tors, the quark gluon coupling in the Dq covariant derivative and the gluon 
self- interactions. QuatLpinch-singularities are removed by the exponentiation 
of the static potentialca The calculation was carried out up to two loops and 
the result reads (k = (p ^ p')) 



'i^TiCFasi^J's) 



Po^^'{K)-{'2Poa^+í3,) In(^) 



a2 



(26) 



where /3o — ll — in£ and (3i = 102 — ^ ng are the one- and two-loop beta func- 



tions, and ai = 



10 



ne (see Refs. M]), as = 456.749 - 66.354 - 



1.235 ?i2 



(see Refs. M]). The constant ne denotes the number of massless quark species. 
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Finite light quark mass effects at two loops are also known0 but are not dis- 
cussed here. A numerical estimate of the three-loop corrections for massless 
light quarks based on Padé approximation techniques has been carried out in 



Ref. |46|. The two-loop matching condition for ay^ is obtained by demand- 
ing equahty of Eqs. (|25| ) and (|26|). The matching conditions for the other 
pNRQCD coefficients shown in Eqs. d22) and which are needed for a 



NNLL description of the QQ dynamics, have been given in Ref. |39|. Because 
the matching proceduře uses off-sheh scattering ampHtudes, the pNRQCD po- 
tential operators basis is different from vNRQCD, where on-sheh matching is 
carried out. 

Appelquist, Dine and Muzinich (ADM)cZl pointed out that at three loops 
the potential of a static QQ pair with relative distance r has an IR divergence 
from graphs of the form in Fig. 13. ADM showed that this IR divergence could 
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Figuře 4: Graphs contributing to the three-loop IR divergence of the QCD static potential 

in perturbation theory. 

be avoided by summing diagrams with an arbitrary number of gluon rungs. 
This summation gives a factor exp{[Vs{r) — Vo{r)]T) for the configuration space 
propagation of the intermediate color-octet QQ pair, which regulates the IR 
divergence. In the pNRQCD framework the diagrams considered by ADM 
are obtained in the intermediate theory NRQCD, where the heavy quarks 
are treated in the static approximation. In Ref. ||4^] Brambilla et al. made 
the important observation that, by construction, the ADM IR divergence is 
equivalent to IR divergences of (unresummed) pNRQCD graphs that describe 
the selfenergy of a quark-antiquark pair due to an ultrasoft gluon (see Fig. 
Thus the static potential defined as the pNRQCD matching coefficient ay^ is 



Figuře 5: Selfenergy of a QQ pair due to ultrasoft gluon radiation in pNRQCD. 
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an infrared-safe quantity. As pointed out in Rcf. , this also nicans that the 
static potential dcfincd via the Wilson loop 

V;tat(r) - ^lirn^ ^ In <J^TrPexp [ - ig£ A^da;^] ^ , (27) 

where P refers to pathordering, is not equivalent to Vs^\r) at short distances, 
because the former contains the fuU ultrasoft selfenergy contributions just men- 
tioned before. 

The anomalous dimensions of the pNRQCD coefíicients are determined 
from UV divergences in QQ scattering diagrams, where again the l/m expan- 
sion with off-sheh static QQ pairs at fixed distance r (or with external jnomenta 
(O, ±p) and (O, ±p') for the incoming and outgoing quark) is usedcJO At this 
stage, loop diagrams with UV divergences that occur in potential loop inte- 
gration, i.e. when potential interactions are iterated, are neglected and only 
UV divergences from ultrasoft gluon exchange are consideredcj In particular, 
the iteration of static potential exchanges between the quarks is treated effec- 
tively in D — A dimensions by using the resummed singlet/octet static QQ 
propagátor ^ {~ido + vj^J {r))~^ , see Ref. [|40|]. 

At leading order in the multipole expansion the entire ultrasoft renor- 
malization of the potential coefficients is obtained from the efFective one-loop 
selfenergy shown in Fig. ||, where the double line represents the static QQ prop- 
agátor. For example, the leading logarithmic (LL) running of the potentials 
cx 1 /m^ is induced by the "one-loop" diagrams where the double line contains 
only one static potential exchange, see Fig. ^ for typical graphs. At this level 




Figuře 6: Diagrams in pNRQCD that contribute to the anomalous dimension of the spin- 
independent 1/m^ potentials. The solid lineš are static quark propagators. The dashed line 

represents the static potential. 



the LL expression for the coupling of the static potential is sufhcient. Because 
the coupling of the static potential to the quarks does not run at LL order in 
pNRQCD, it is just as{fis)- Due to gauge cancellations only the p.A/m cou- 
pling of the ultrasoft gluons to quarks (which only runs trivially with as{fJ-) 
at LL order) contributes and generically leads to the following LL anomalous 
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Figuře 7: Diagrams in pNRQCD that contribute to the three-loop anomalous dimension of 
the static potontial. The solid lineš are static quark propagators. The static potential does 
not run at the one- and two-loop level. 



dimension of the spin-independent í/rri^ coeííicients, 

J^[D^{'}{^i,^is),D2{^^,^^s)} ^a4^i)a4^is). (28) 

At this order, the spin-dependent potentials do not run in pNRQCD. 

A more involved example is the NNLL running of the static singlet po- 
tential av,, which was considered in Refs. [ ^ , ^ . An ultrasoft anomalous 
dimension of ay, is generated for the first time by UV divergences in the 
"three-loop" diagrams in Fig. ^ which are also contained in the one-loop self- 
energy of Fig. |^. The ultrasoft loop with the gluon propagátor is UV-divergent. 
The two potential insertions are obtained by expanding the static QQ propa- 
gátor in terms of the static potential. At this order its coupling can again be 
approximated by as{fis)- In momentům space the computation consists of an 
ultrasoft gluon loop in Z? = 4 — 2e dimensions, where the coupling of ultra- 
soft gluons to the static potential is Q!s(/í), and static quark loops in Z) = 4 
dimensions, where the couplings of the potentials to the quarks are as{^J-s)■ 
The incoming and outgoing quark momenta are (0,±p) and (O, ±p'), respec- 
tively. The ultrasoft loop is UV-divergent, and the static quark loops lead to 
the overall factor l/(p — p')^. For a color singlet QQ pair the amplitudě reads 



Figs.7 = -i 



6 (P-P')' 



1 1 



(29) 



The IR divergence oc 1/ eiR is equivalent to the ADM infrared divergence men- 
tioned above. The anomalous dimension for the static singlet potential ob- 
tained from the UV divergence reads 



aVs(^J■, ^J■s) ^ -T— [avjps)] as{^J.) . (30) 
dm /i Iítt 

Using the initial condition of Eq. (EO) the NNLL result for the static singlet 
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Figuře 8: Diagrams in pNRQCD that contribute to the two-loop anomalous dimension of 
the ^ potential (a) and the QQ production current (b). 



potential at the scale ^ = fiu reads 



two-loop 



+ 



(31) 



where [av^ (/^s)]two-ioop is the two-loop matching result obtained from Eqs. (25) 
and dl). 



^.^ Potential Loops in pNRQCD 

If pNRQCD is apphed to dynamical quarks, the counting Dq ~ ~ ^ ~ it^t 
is reinstated and the potential region is reintroduced into the theory. As men- 
tioned before, it is claimed in Refs. [ p^ , p8||4Č| ] that the potential region is always 
contained in the pNRQCD framework, but hidden in a subtle way as long as 
the l/m power counting is imposed. Potential loops, which are generated by 
iterations of potentials, can lead to UV divergences. Consider for example 
the QQ scattering diagram in Fig. ^a with two D^f l potentials and one static 

(2) 

potential. The D^^ ^ potential is a á-function in configuration space, and the 
graph is UV divergent inducing a two-loop next-to-leading logarithmic (NLL) 
anomalous dimension to D^fl- The diagram in Fig. ||b describes QQ produc- 
tion and induces a two-loop NLL anomalous dimension to the coefficients of the 
QQ production currcnts such as the '^'^^^Lj = ^5*1 current O^^ — cr{ia2)x- 

(2) 

Thus the LL running of ^ and ay^ obtained in NRQCD and from ultrasoft 

(2) 

gluon radiation in pNRQCD mixes into ^ and c^^ at NLL order through 
UV divergences in potential loop diagrams. In order to quantify this mixing a 
correlation between the soft scale /x^, the ultrasoft scale fiu and the scale for 
potential loops, Hp, needs to be imposed. Such a correlation does, by itself, 
not exist in the pNRQCD framework, since the quarks are treated as static. 
So the correlation needs to be imposed by hand, once potential loops and the 
counting Dq ^ ^ 2_ ^ M_ are reintroduced into the theory. 

From the physical point of view this correlation should reflect the energy- 
momentum relation of the heavy quark: /Xp ~ /is ~ ^/J^jň. This concept has 
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been first realized in the framework of vNRQCD,liÍ where the correlation of 
scales arises naturally from the thcory itsclf. In Rcf. Pineda proposed the 
foUowing prescription for pNRQCD: after NRQCD is evolved to ^ = and 
the uhrasoft pNRQCD evolution is carried out down to the scale /x = /i„ as 
described before, the relations 



fJ-s 



fJ-p , l-íu 



m 



(32) 



are imposed. This scale correlation has to be used for the determination and 
the integration of the RGE's obtained from potential UV divergences as well 
as for the determination of matrix elements, i.e. in particular in the solution 
of the Schrodinger equation. For the determination of all matrix elements the 
1 /m expansion is abandoned and the non-relativistic power counting from the 
regions in Eq. ( [Til ) is reestablished. For the integration of the renormalization 
group equations fip has to be run from the hard scale m down to the soft 
scale of order mv. This proceduře means in particular that, for the description 
of dynamical quarks, the end points of the running in NRQCD (/i^) and the 
ultrasoft running in pNRQCD (/i„) are generally not the soft and the ultrasoft 



scales of order 



and 



respectively. It is claimed in Ref. |50| that this 



prescription sums properly all large logarithmic terms. 

The NLL anomalous dimensions induced by Figs. ^ and b have the generic 
form 



lnM;^'í'Am'^^ 



D 



(2) 



\ m 



InyCtp 



2 2 

Ci,xif^p) ~ c^xi^^p)DdU — '^^p) "vj — • (33) 



At this order, only the LL results need to be known on the RHS. The fuU NLL 
running of ^ has actually not been calculated so far; it would be relevant for 
the description of the QQ dynamics only at N'^LL order. The NLL runnning 
of in the pNRQCD framework has been determined in Ref. 1 50 1 . 



^.5 O-pen Issues 

There are issues in the pNRQCD framework that might require some investi- 
gations. One obvious point is that the matching procedures in the scheme of 
Eq. ( p^ appear to strongly rely on specific choices for the kinematic situation 
where the matching conditions can be calculated. It would be desirable that 
the matching proceduře could be carried out for any kinematic situation, for 
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which thc efTective theory is applicablc. For cxample, the NRQCD matching 
is carried out for on-shell scattering amplitudes with external quarks exactly 
at the threshold point, where the quark velocity is zero. If on-shell scatter- 
ing amplitudes with some sniall but non-zero quark velocity are employed the 
diagrams in fuU QCD also contain contributions from, leťs say, the poten- 
tial momentům region, which cannot be reproduced in the l/m expansion 
employed for the NRQCD computations. 

Another issue concerns the prescription of Ref. for the treatment of 
UV divergences from potential loops after the endpoints of the NRQCD and 
ultrasoft pNRQCD running, yUs and are correlated to the potential scale, 
tJ'P = ^J's — ^/JJmŤŤí. The potentials that arise order-by-order from matching pN- 
RQCD to NRQCD are independent of the matching scale fis up to highcr order 
terms, which are small corrections in the matching conditions for Hs ~ (p — p')i 
see for example Eq. ( p6| ) for the static potential. On the other hand, in the 
anomalous dimensions for D^^l and c^x Eqs. ( ^3| ) the residual dependence 
of for example av^ on /x^ leads to the summation of unphysical logarithms 
unless the corresponding logarithmic terms in the matching computation for 
the potential are summed too. This means in principle that the higher order 
matching corrections to the potentials cannot be considered as small when 
correlated running according to Rcf. is employed. 

It is a prediction of the pNRQCD framework that the Coulomb potential 
in thc Schrodinger equation that describes dynamical QQ pairs is_equal to the 
static potential determined in the l/m expansion in pNRQCDc3 The ultra- 
soft anomalous dimension of the static potential is determined from Eq. (p9|), 
which does not contain any scale dependence from the static ladder diagrams 
that iterate the potentials. On the other hand, after the scale correlation 

= = ^/Jhjň and the non-relativistic power counting are imposed, the 
corresponding UV divergent contribution from the diagram in the N'^LL ma- 
trix element describing leťs say on-shell QQ scattering in analogy to Fig. |^ 
reads 

. CfC\ [as{iis)An''[as{tiu)pl'] 
6 (p-pO" 

where the dependence on and /i„ has been made manifest. Here, the scale 
dependence of the potential loops is included, because the quarks are not static. 
The 1/e term is canceled by the counterterm of the Coulomb potential deter- 
mined from Eq. (p9|). However, only a part of the In/ip terms is canceled by 
the order ultrasoft logarithm summed into a\/ by Eq. (^^. The remain- 
ing In /ip terms seem to remain uncanceled in the hnal N'^LL result for QQ 
scattering. 



(34) 
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4-6 Aqcd ~ mv'^ or larger 

So far pNRQCD has been discussed for the case m ^ mv 3> mv^ ^ Aqcd, 
where the counting as{m) ~ as{rav) ^ as{mv^) <^ 1 was used. In this situa- 
tion the description of the QQ dynamics is predominantly perturbative and the 
effective theories serve mainly as a tool to handle the non-trivial perturbative 
calculations and summations in a manageable and systematic fashion. Non- 
perturbative efTects appear as local gluon and Ught quark condensates in an 
expansion in Aqcd/^t-í^^, see Refs. One motivation in the construction 

of pNRQCD was that its multi-layer structure should make it also apphcable 
to Systems where Aqcd is of order mv^ or larger. In this case also the non- 
relativistic dynamics of the charmonium or of higher bottomonium excitations 
could be studied in a systematic and quantitative manner. Qualitative discus- 
sions of the application of pNRQCD in cases where Aqcd is of order mv^ or 
larger have been given in Ref. |^^. For mv > Aqcd ^ 'mv'^ pNRQCD can still 
be matched perturbatively to NRQCD at the soft scale /Zs, but Aqcd screens 
the ultrasoft scale (or any other lower scale) and as{mv^) is of order 1. How- 
ever, the coupling of ultrasoft Aqcd) gluons to the quark-antiquark pair 
c>c r.E is still small because of the multipole expansion. The effects of ultrasoft 
gluons enter as non-perturbative non- local condensates, that are small correc- 
tions. For Aqcd ^ mv, the matching of NRQCD to QCD at the hard scale can 
still be carried out perturbatively, but the full non-relativistic QQ dynamics is 
screened by Aqcd- In this case NRQCD (Sec. ||) should describe QQ produc- 
tion and decay, but all non-relativistic QQ dynamics is non-perturbative. 



5 Effective Theories III: Velocity NRQCD 

The effective theory vNRQCD was first proposed by Luke, Manohar and 
RothsteinOJ It is devised for systems where the relation m 3> mv ^ mv"^ ^ 
Aqcd is valid. Its main fields of application are top pair production close to 
threshold, sum rules for T mesons and maybe low radiál excitations of bot- 
tomonia. The theory vNRQCD provides the separation of all different degrees 
of freedom in a single step at the hard scale and is defined in dimensional 
regularization. 
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5.1 Basic Idea 

The basic idea behind the construction of vNRQCD is that all resonant degrees 
of freedom (dof 's) (i.e. dof 's that can develop a particle pole in the momentům 
regions of Eq. ([Tl|)) are separated and that all ofF-shell fluctuations (i.e. fluc- 
tuations that cannot come close to their mass-shell in the regions of Eq. (pl])) 
are integrated out at the hard scale. This means that the theory contains dif- 
ferent gluon field operators for the gluonic fluctuations with soft and ultrasoft 
energies and momenta. The theory vNRQCD is matched to QCD at the scale 
m and the renormalization group equations are constructed such that large 
logarithms of the scale ratios m/(p), m/ (E) and {p)/{E) are summed simulta- 
neously into the coefhcients of the operators when the theory is scaled down. 
This is achieved by having different renormalization scales for potential/soft 
and ultrasoft loop integrals, which are correlated in accordance with the quark 
equations of motion. The scale correlation is fixed by the theory itsclf and not 
imposed by hand. The correlation of the scales and the renormalization group 
evolution are expressed in terms of the "velocity" scaling parameter v. The 
hard scale corresponds to = 1 and all large logarithmic terms are summed 
for v = v. 



5.2 Effective Lagrangian 

The effective Lagrangian is formulated in terms of quark and gluonic fields de- 
scribing modes that can become on-shell in the momentům regions of Eq. (pŤ|). 
For convenience, massless quarks will be referred to as gluons in the foUowing 
presentation. Heavy quarks are simply referred to as quarks. On-shell fluctua- 
tions for quarks exist only for the potential regime where (/cojk) ~ [mv^jinv) 
and on-shell gluonic fluctuations exist in the soft and the ultrasoft regime where 
(fco,k) ~ {mv,mv) and (fco,k) ~ {mv'^ , mv"^) , respectively. These three kinds 
of modes are referred to as "potential quarks", "soft gluons" and "ultrasoft 
gluons". Although soft gluons cannot be produced in a non-relativistic QQ 
systém with total energy of the order mv"^, they are kept as relevant dofs in 
vNRQCD. This is necessary because interactions with soft gluons are involved 
in the renormalization of vNRQCD operators. In fact, one can consider vN- 
RQCD as a theory that can describe interactions of a quark with soft gluons 
as well, such as in Compton scattering. 

In the effective Lagrangian only ultrasoft energies and momenta are treated 
as continuous variables. Soft energies and momenta are treated as discrete 
indices for potential quarks and soft gluons. For a heavy quark this is achieved 
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by writing its energy E and three-momcntum P as 

E ^ k\ P==p + k (35) 

where the three-vector p is of order of the soft scale mv and the four-vector 
k — (fc°,k) is of order of the ultrasoft scale mv'^. This can be understood 
as subdividing the fuU non-relativistic momentům space for the quark three- 
momentum (with side length of order mv) into small cubes with ultrasoft side 
length of order mu^, such that the number of ultrasoft cubes is (^^)^. The 
center location of each ultrasoft cube is given by the index p and a point 
within an ultrasoft cube is labeled by k, see Fig. Only the ultrasoft four- 

mv ■- 



mv" 

Figuře 9: Momentům space of size mv is divided into boxes of size mv^ . A point in 
momentům space is labeled by p and k. 

momenta k are treated as continuous variables, the soft three-momentum p is 
a discrete index. The number soft three-momenta is equal to the number of 
ultrasoft cubes. The quark and antiquark fields in vNRQCD are written as 
two-componcnt spinor fields 

V'p(a^), Xp{x), (36) 

where p is the soft index and x the continuous variable that describes ultrasoft 
fluctuations. This proceduře is similar to HQETJ^ where the four-momentum 
of a heavy quark is split into mv^^ of order m and a residual four-momentum 
k of order Aqcd, = mv^ + /c^. The index v^^ is a discrete variable one 
has to sum, whereas k^ is a continuous variable one has to integrate. In 
vNRQCD one has to sum over p and integrate over x (or k in momentům 
space representation) . 
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The energy E and three-momentum P of a soft giuon is writtcn in an 
analogous way as 

E ^ p^^ + k^, P = p + k, (37) 

where p = {po,p) is a four-vector of the order of the soft scale mv and 
k = (/cojk) is a four-vector of the order of the ultrasoft scale mv'^. The cor- 
responding field in vNRQCD that annihilates and creates a soft gluon with 
index p in vNRQCD is Al^{x). One has to sum over p and integrate over x (or 
k). The energy E and the three-momentum P of an ultrasoft gluon is just 

E , P =k, (38) 

where k = (fco,k) is a four-vector of the order of the ultrasoft scale mv^ . The 
corresponding field in vNRQCD that annihilates and creates an ultrasoft gluon 
vNRQCD is A^'{x). 

The decomposition into a soft and an ultrasoft momentům component of 
leťs say the quark is not unique. Taking an ultrasoft three-momentum q one 
can redefine 

k^k + q, P^p-q. (39) 

This reparametrization invariance is similar to HQET0 but it does not affect 
the spin axis of the quarks, because the latter is fixed by the choice of the 
center of mass frame, which is not affected by the transformation in Eq. (^9|). 
So the consequences of reparametrization invariance in vNRQCD are the same 
as for HQET for spin-zero particles. The basic outcome of reparametrization 
invariance in vNRQCD is that derivatives of il^p{x) or Xp(a^) are always of the 
form (ip -|- V). 

The effective Lagrangian is written in terms of quark fields, i/^pj antiquark 
fields, Xp, soft gluon fields, A^, and an ultrasoft gluon field, A^^. The covariant 
derivative is 1?^ = (-D°,D) with D° = d° + igA°, D = V-z.gA, and it involves 
only the ultrasoft gluon field. The effective Lagrangian is manifestly gauge- 
invariant with respect to ultrasoft gauge transformations, i.e. with slowly vary- 



ing gauge phases involving frequencies of order mv^ . It is argued in Ref. |14 
that the fuU gauge invariance including also higher frequencies of order mv is 
recovered by combination of reparametrization invariance and ultrasoft gauge 
invariance. A mathematical proof of this conjecture does, however, not yet 
exist. 

In the center of mass frame the Lagrangian has the formEj'Ěil'Ě3 



4 

V'p + (V' ^ x) 
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(P'-P 

■^(p^pOV-^-^pxlp-x-p 



E 



where (k = (p — p')) 



y(p,p') = (T-^^Ť-^) 



v, 
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U, 



(0) 



(P - P')' 



(40) 



V. 



(T) 



(41) 



(42) 



The matrices and Ť'^ are the color matrices for the 3 and 3 representations, 
igG^'^ = [D^, D"] and m is the heavy quark pole mass. Interactions involving 
ghosts and light quarks are not displayed. The effective Lagrangian (second 
line) contains interactions betwcen quarks and soft gluons. Due to energy 
conservation at least two soft gluons participate in interactions with quarks. 
Graphically the interaction of a quark with two soft gluons is represented by 
Fig. nOfi. The explicit form of the leading order soft interactions is displayed 

(0) 



in Eq. ([42|). The í/q term that is visible for example in Uqq arises from the 



00 

propagation of an off-shell "soft quark" , see Fig. ^ The term does not have 
aií iS prescription and does not lead to poles or pinch-singularities when the 



Heavy Quarkonium Dynamics 



36 



(a) (b) (c) 




Figuře 10: The Compton scattering diagrams (a,b,c) in full QCD generate the soft gluon 
coupling (d) in vNRQCD. Soft gluon modes are denoted by a zigzag line 



loop integrál over is carried out by contours. Terms in the expansion of full 
QCD such as 

1 1 



fco + i5 ko — i6 

lead to additipnal four-quark interactions with soft gluons in the matching 
procedureEjl£j There are terms in the Lagrangian where ultrasoft gluons cou- 
ple to potential operators. The term in the third line of Eq. (^) for example is 
generated by a QCD diagram similar to the third diagram in Fig. ^. There are 
potential-type four-quark interactions (fourth line) that depend non-locally on 
the soft indices, but are local interactions with respect to the dynamical ultra- 
soft fluctuations. They are graphically represented by Fig. [ll|b. In Eq. (41) 



Figuře 11: The scattering Born diagram (a) generates the potential four-quark interaction 

(b) in vNRQCD. 

the leading order Coulomb potential interaction and the v- and w^-suppressed 
potentials are displayed. The basis in terms of (1 ® 1) and (T'^ ® Ť^) can be 
converted to the color singlet and octet basis by the transformation 



^inglct 




Voctct 





-Cf 



Vt^t 



(43) 



The multipole expansion has to be applied strictly for all terms in the effec- 
tive Lagrangian in order to achieve separation of the modes in the diíferent 
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regions. For examplc thc term (p — íD)^ in the bilinear quark sector has to 
be considered as an expansion in D because íD <^ p. Interactions between 
soft and ultrasoft gluons without (heavy) quarks do not exist because ultrasoft 
modes cannot resolve the "high frequency" soft gauge transformations (that 
are reaUzed through reparametrization invariance). On the other hand, the 
soft gauge fielda^feel uftrasoft gauge transformations just as a globál phase 
transformationÉJ 

The velocity power counting of the fields can be derived from demanding 
that the action for the kinetic terms is of order and reads 



p 


^P, Xp 






D 




v 


„3/2 


V 









(44) 



Soft and ultrasoft massless quark fields count as u'^/^ and v^, respectively. 

External currents describing quark production, annihilation or decay are 
also built from the quark fields in Eq. (|36|). For example, consider QQ produc- 
tion at NNLL order in e^e~ annihilation. One needs the ^Si (vector) current 
with dimension three and five, 

Op,i = ipp^ (7{ia2)x-p, 

Op.2 = ^ V^p^ P^t(íct2) X-p , (45) 
and the ^Pi (axial- vector) current with dimension four 

Op.3 ^ ^'4'pHo-,fr ■ p]{ia2)x-p- (46) 

The currents describe production of a QQ pair with momentům ±p. The 
corresponding annihilation currents Oj, are obtained from complex conju- 
gation. 

5.3 Loops and Correlation of Scales 

A generál multi-loop graph in vNRQCD is divided into soft, potential and 
ultrasoft loop subgraphs. The diagram in Fig. |l2| a shows a typical ultrasoft 
loop. Because an ultrasoft gluon cannot change the soft index of a quark the 
loop does not involve any sum over soft indices of the quark. In dimensional 
regularization the integration measure of the ultrasoft loop reads (í? = 4 — 2e) 



(47) 
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Figuře 12: Typical one-loop diagrams in vNRQCD containing ultrasoft (a), potential (b) 

and soft (c) loop integrals. 

where fiu is the common renormalization scale for the ultrasoft dynamical 
fluctuations introduced to keep couplings dimension-less. In vNRQCD nu is 
rewritten as 

fiu = miy^ , (48) 

where v is called the "velocity scahng parameter" and v = 1 corresponds 
to fiu = m. The reason for introducing v will become clear below. The 
natural choice for the vefocity scahng parameter in matrix elements v ^ v. 
The integrál for the ultrasoft loop is carried out over the fuU D-dimensional 
non-relativistic space. The information that the loop is dominated by the 
ultrasoft region is implemented by the pole structure of the propagators and 
the multipole expansion. 

In Fig. ^2|b a typical potential loop diagram is displayed. Its integration 
measure reads 

p 

and involves an integrál over an ultrasoft momentům and a sum over all pos- 
sible soft indices of intermediate quark-antiquark pairs. In practical analytical 
calculations it is convenient to rewrite the combination of the sum and the 
ultrasoft integration as an integrál over the fuU D-dimensional non-relativistic 
momentům space. The number of terms in the sum is {j^Y, where 

US = mv . (50) 

Thus Eq. (^ ) can be written as 




and we find that vNRQCD itself "recognizes" that the proper renormaliza- 
tion scale for a potential loop is the soft scale fis- The pole structure of the 
propagators and the multipole expansion determine that the loop is dominated 
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by the potcntial region. For a typical soft loop such as in Fig. the same 
renormalization scale iis arises, when the combination of sum and ultrasoft in- 
tegration is written as an integrál over the fuU _D-dimensional non-relativistic 
momentům space, 




The pole structure of the propagators and the multipole expansion determine 
that the loop is dominated by the soft region. For more complicated multi-loop 
diagrams the same rules apply. Soft and potential loops have the renormaliza- 
tion scale fis and purely ultrasoft loops the renormalization scale iijj- Whether 
a loop is either soft, potential or ultrasoft is determined by the pole structure 
of the propagators and the multipole expansion. 

In vNRQCD there are two renormalization scales, which are correlated 
such that a "subtraction velocity" v can be defined. The parameter v serveš 
as the natural scaling parameter for all types of loop integrals in vNRQCD. 
The scale correlation is not imposed by hand, but an intrinsic property of the 

2 

theory itself. The physical origin of the correlation fj,u = ^ is the equa- 
tion of motion of the quark, which relates the soft and ultrasoft energy scales. 
The renormalization scales can be naturally incorporated in the effective_La- 
grangian in Z? = 4 — 2e dimensions to keep the couplings dimensionlessĚá In 
generál, a vertex with 2 + i quarks, j soft gluons and k ultrasoft gluons receives 

a factor [ms"^^-'^^']- 

Without the correlation of scales vNRQCD becomes inconsistent. Con- 
sider, for example, the three-loop diagram in Fig. ^3|a, which contributes to the 
renormalization of the dimension-three '^Si QQ current Op,! = ip^cr{i(T2)x*-p- 
The operátor Op,! describes for example the leading order production of QQ 




Cl Vc Vc Cl SV2,r Vc 



Figuře 13: Examples of graphs that contribute to the three-loop anomalous dimension for 
Cl- SV2,r are one loop counterterms. The coupling of the ultrasoft gluon to the quarks is 
gp.A/m, where the scale of g is mi/^. 

pairs in e+e" annihilation. After subtraction of all UV subdivergences for 
e ^ O through the diagrams in Fig. [l^b the remaining overall UV divergence 
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leads to thc foUowing countcr term contribution for thc coefficient ci of the 
operátor Op,i, 



<5ci ~ a,(mí.2)[vf)(i.)]2 



Í + -( \n ^iu -2ln fis + # 



e- 



(53) 



where # symbolizes real numbers. The anomalous dimension for ci can be 
determined|-U|nambiguously only if the correlation between fj.u and /ig is ac- 
counted forĚII 

The computation of matching conditions and anomalous dimensions pro- 
ceeds in the canonical way taking the velocity scaling parameter as the fun- 
damental renormahzation scaUng variable in dimensional regularization. This 
is thc origin of the "v" in vNRQCD. 



5.4 Power Counting 

The velocity power counting can be carried out using the velocity scaling rules 
of energies and momenta in the soft, potential and ultrasoft regions of Eq. (pŤ|). 
At this point the power counting is similar to the one in the threshold expan- 
sion and relies on assigning powers of v to all loop measures, propagators 
and vertices in vNRQCD diagrams. An equivalent but more convenient power 
counting prescription has been derived in Ref. [p^ . It relies only on í;-counting 
for the vertices of the effective Lagrangian and the topology of the soft momen- 
tům component in loop diagrams. Let vl^^\ vj.^^ and Vjf^^ denote the number 
of soft, ultrasoft and potential vertices of order in a given graph, where the 
velocity counting of the fields in Eq. (^^ is included. Here, ultrasoft vertices 
involve only ultrasoft fields, potential vertices involve at least one quark and 
no soft fields, and soft vertices at least one soft field. For example, the vertices 

have Vg , and , respectively. The ultrasoft gauge kinetic term 

C^Gpjy has V^^\ For soft and potential vertices vertices k > and for 
ultrasoft vertices fc > 8. A given diagram is then of order , wherellá 

= 5 + E [(^ - 8)^^"" + - ^^^k^^ + - 4)yf ^] - Ns , (54) 



Ns being the number of connected soft components of the diagram. Ultrasoft 
and soft vertices can only give positive contributions to the sum. Potential 
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vertices give positive contributions except for thc Coulomb potcntial intcrac- 

í p) 

tion, which has V^; . An insertion of n Coulomb potentials lowers 5 by n, 
but leaves the w-counting unchanged if one also counts the Coulomb coupling 
as Vc ~ Qís ^ v. Thus each additional insertion of a Coulomb potential gives 
a factor asV~^. It is therefore convenient to count the Coulomb potential 
oc Qís/k^ simply as of order asV~^, the l/m|k| potentials as of order a^w" and 
the 1/m^ potentials as of order agV^, ete.. 

For illustration let us consider the graph displayed in Fig. [ij, which has 
V^^'^ =2andV:i^^^ =2- Thus one obtains (5 = 5 + (6-5) x 2 + (4-5) x 2 = 5, so 
the graph is of order v^. Dividing this result by the velocity counting of the 




Figuře 14: Examples of graph or order a^v^. The coupling of the ultrasoft gluon to the 

quarks is gp.A/m. 



four external quark fields (w^), we find that the amputated graph is of order 
a^v~^. The same result is found from the u-counting of the loop measures, 
propagators and vertices based on the velocity scaling of the ultrasoft and 
potential regions in Eq. (11). The diagram in Fig. |lj is suppressed by agv'^ 
with respect to the same graph without the ultrasoft gluon, which is of order 
a^v^. Each additional ultrasoft gluon gives an additional factor a^w^. If we 
apply the same counting rules to QED we find that the Lamb-shift, which is 
caused by ultrasoft photons, is suppressed by a factor av^ ~ with respect 
to the leading Coulomb interaction. 



5. 5 Matching 

The determination of the matching conditions of the coefRcients in the effective 
Lagrangian at thc hard scale, v — 1, is carricd out with amplitudcs for on-shell 
quarks and gluons. This avoids operators that vanish by the equation of motion 
and gauge-dependent matching conditions. Any kinematic situation described 
by the effective theory can be used for the matching calculation. For example, 
at the Born level the matching-conditions for the potentials obtained from the 
on-shell graphs in Fig. ^ areEJ 

Vf )(1) = ina.im) , V^íl) = 47ra.(m) , Vf )(1) = , 
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Vr(l) = -6.a.M, Vr(l) = -^. (55) 

Matching conditions that vanish and annihilation contributions are not shown 
in Eq. ( |55| ) . The on-shell matching conditions for the potentials at the one-loop 
level were computed in Ref. ||5^ . For the order 1 / (rnk) potentials the results 
read 

Vr(l) = alim){^ - ^) , Vi'\l) = alirn)^ , (56) 

where Cd = SCp - 3Ca and Ci = Cf{^Ca ~ Cp)- For SU(iV) we have 

iV2 - 1 -A -l 

Cf = -^^TT- , Ca = N, Cd^ — — , and Ci = 



2N ' ' iV ' 4iV2 ■ 

For the description of the QQ dynamics at NNLL order the results in Eqs. 
( ^ ) and ( ^6|) are sufficient. The coefScients of the Coulomb potential do not 
receive any higher order matching corrections up to order , see Ref. . In 
generál, a non-zero matching corrections appears, when there is a an off-shell 
region such as the hard one that contributes in the matching condition. 

The Coulomb potential ^ is equivalent to the potential oc that de- 
scribes the QQ binding effects in the Schrodinger equation only at the lead- 
ing logarithmic level. The known one- and two-loop corrections obtained in 
Refs. 1^,0 are not contained in In vNRQCD these corrections arise in 
form of time-ordered products of the soft vertices in Eq. (^) , such as the one- 
loop diagram in Fig. ^2|c, i.e. they are contributions from matrix elements and 
not part of the matching coefHcient of the Coulomb potential. 

The matching conditions for the two-quark soft vertices are determined by 
Compton scattering diagrams in fuU QCD. At Born level the diagrams in Fig. 
p^,b,c are matched onto the local soft operators in the second line of Eq. (pO|). 
At leading order in v (í/^ií""'') the results are displayed in Eq. (|4^). The results 
up to ncxt-to-next-to-leading order in v {U^1~^'^'^) have been determined in 
Ref. . Interestingly, the higher order corrections to the two-quark soft ver- 
tices can be determined by computing the corrections in the HQET Lagrangian 
(see e.g. Ref. [^l|), ^uad then matching their timc-ordcred products to the two- 
quark soft verticesO This works because only the soft momentům region is 
relevant for the renormalization of the two-quark soft vertices. One outcome 
is that there are no higher order perturbative corrections to the leading order 
(cr = 0) soft vertices in Eq. (^). 

The matching conditions for external currents are determined in the same 
way as the matching conditions for the potentials. For example, for the de- 
scription of QQ production at NNLL order in e^e^ annihilation one needs 
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Cl Vfc 



Cl . Vc Vc 



Cl V2,s,r Vc 



Cl Vc V2,s,r C2 Vc Vc 



Figuře 15: Difference of fuU QCD and order a'^v^ vNRQCD graphs which gives the two 
loop matching for ci(l). The X denotes an insertion of the p'^/{8m^) operátor, and graphs 
with this insertion on a different propagátor are understood. 



the matching conditions for the ^Si vector current up to dimension five, 
Jp = ciOp_i + C20p.2 and the ■^Pi axial-vector current up to dimension four, 
Jp = C30p_3, see Eqs. (^Šj) and (|4^). The matching condition for ci needs to 
be known at order a^. For C2.3 Born matching is sufficient and is obtained by 
expanding the fuU QCD currents Q^^Q and Qj^jsQ at tree level, 



C2(i^ = l) = -l/6. 



(57) 



The one-loop matching condition for ci is well known and scheme-independent 
for mass-independent regularization schemes. The two-loop computation in- 
volves computing the difference between graphs in fuU QCD and in the effective 
theory as shown in Fig. 15. The two-loop vertex corrections in the fuU theory 
were computcd in Rcfs. [33,56 . The matching can be carried out at the level 
of amplitudes as shown in Fig. 15, which involves a cancellation of IR divergent 
Coulomb phases, or at the level of the total cross section, where the Couloiah 
phases are absenL The latter type of matching is called "direct matching" E3 
The result readsĚIl 



ci(l) = 1 as[m)^ 



:(m) 



Ci 



'In 2 25 2 \ / \ K 

^ ^^^^ 

whereEJ k = -2.556 + 0.08256 7i£, and rií is the number of light quark species. 
At order the matching condition is dependent on the subtraction scheme 
and on the definition of the operators in the effective theory. The result in Eq. 
( ^ ) is obtained in the MS scheme and with the form of all operators given in 
D — A dimensions. The result is different from Ref. [Q, where the potentials 
were defined with an explicit dependence on _D = 4— 2e, such that the matching 
conditions agrees with the contribution of the hard region in Eq. ( pŤ|) in the 
threshold expansion. 
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5.6 Soft IR Divergences 

One of the most important (and interesting) conceptual aspects of vNRQCD 
is the existence of the ultrasoft gluons at the hard scale through the relation 

2 

of scales nu — ^ = rav^ . This relation is an intrinsic property of the effective 
theory. One consequence is for example that ultrasoft gluons contribute in 
the determination of the matching conditions at z/ = 1, or the running of the 
coefRcients for any value for v smaller than 1. The necessity for having the 
ultrasoft degrees of freedom separated at the hard scale can be seen explicitly 
in multi-loop diagrams, where potential (or soft) and ultrasoft loops occur at 
the same time. An unambiguous determination of the anomalous dimcnsions 
is only possible, if the scale correlation is taken into account, see e.g. Eq. (p3|). 

This seems to be in contradiction to the fact that the ultrasoft gluons 
clearly arise from splitting the originál gluon ficld into soft (~ mv) gluons and 
ultrasoft (~ mv^) gluons. The latter clearly describe fluctuations at scales 
much smaller than mv. For the threshold expansion the presence of the ultra- 
soft momentům regions and the necessity for a proper expansion (separation) 
of the ultrasoft momentům region for any choice of renormalization scale is 
unproblematic conceptually - the threshold expansion is a technical prescrip- 
tion to obtain an expansion of Feynman diagrams and not an effective theory, 
where renormalization issues are relevant. 

In vNRQCD the problém is resolved by the property of IR divergences in 
soft loop integrations. The mechanism to account for the scale correlation = 

2 

^ = mv^ and the fact that soft and ultrasoft degrees of freedom fluctuate 
at different Jfingth scales is to treat the IR divergences in soft loops as UV 
divergences E3Ě3 This means that soft IR and soft UV divergences contribute 
to the anomalous dimensions of the operators in the effective Lagrangian. For 
illustration let us discuss the following simplified situation where only diagrams 
with single \ poles and either purely soft or ultrasoft gluons are discussed 
and real Coulombic IR divergences are droppedEj A generál diagram with a 
divergent soft loop then gives an amplitudě with the divergence structure 

lA^ ^ + — + C = + (C - S) 

euv eiR V^uv eiR/ ^uv \euv eiR 

(59) 

Pure dimensional regularization does not distinguish between UV and IR di- 
vergences, but an Identification of leťs say UV divergences can be made by 
performing the calculation with an additional IR regulátor. The term C repre- 
sents diagrams involving scaleless integrals such as tadpole graphs, and A and 
B represent all other graphs. For example, for the hydrogen Lamb shift (see 
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Sec. only C is non-zero, while for positronium or quarkonium also A and 
B are non-zero. At the same order in the power counting as Eq. ( p9[ ) there is 
an amplitudě with a divergent ultrasoft gluon loop, which has the form 

= + (60) 

In generál, D is independent of (C — B) since ultrasoft loops are in generál 

not proportional to (— —). Examples are the two- and three-loop dia- 

grams contributing to the renormalization of Vk (Ref. and Vc (Rcf. [p5| ), 
respectively. The IR divergence matches with an IR divergence in fuU QCD. 



The non-trivial issue is that the coefRcient of the UV divergence in Eq. (60) 
matches with the IR divergence in the soft amplitudě in Eq. (|59|). This state- 
ment can be understood intuitively from the fact that the originál gluon field 
is split into soft and ultrasoft modes in the effective theory. Thus any soft 
IR divergence cannot be an IR divergence of the fuU theory, because it is not 
associated with scales below mv'^. In other words, in Eq. ( |59| ) the euv's cor- 
respond to the scale m and the eiR's correspond to the scale mv, while in Eq. 
( |60| ) the euv's correspond to the scale mv and the eiR 's correspond to the scale 
TOu^. The validity of this argument for complicated multiloop diagrams, where 
soft, potential or ultrasoft loops exist simultaneously, is conjectured, and not 
yet proven mathematically. However, examples such as in Eq. (p3|), seem to 
support its validity. 

The running of the operators is determined by UV divergences. However, 
examining iA^ + iA^ one finds that the term (C — B) in the soft amplitudě 
acts like a tadpole contribution that puUs the in the ultrasoft amplitudě 
up to the hard scale. It is argued in Ref. |5^] that the scale-dependence of 
the coefficient (C — B) in Eq. (^0|) does not affect this argument since the 
scale-dependence of the (C — i?) term is Eq. (|59| ) can be chosen arbitrary. 
The hnal outcome is that interpreting all - terms in the soft amplitudě as UV 
divergences, running the ultrasoft modes from m to mv^ with an anomalous 
dimension oc (C — B) and running the soft modes from m to mv with an 
anomalous dimension oc {A + B) correctly performs the running between the 
scales. This is the basis of computations of the evolution of the couplings and 
the summation of large logarithmic terms in vNRQCD. 

The correspondence in Eqs. ( ^9|) and ( |60|) provides a useful tool. If the 
UV divergences (A + C) in the soft diagrams are known, and the combination 
{B — C) is determined from the UV divergences in the ultrasoft diagrams, one 
arrives at {A + C) + {B — C) ^ A + B, which is the combination needed to 
determine the soft anomalous dimension. This means that one can recycle 
results for the running of coefficients in HQET for the renormalization of two- 
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Figuro 16: Order ct^/v diagrams with potontial itorations. The X donotos an insertion of 
the p'*/8m^ relativistic correction to the kinetic term, Vjg) stands for a potential and 

V(i) for a 1/m^ potential. 

quark soft vertices in vNRQCdEÍ An explicit demonstration for the mechanism 
described above is given in Sec. |^. 

5. 7 Renormalization Group Equations and Running 

The renormahzation proceduře is carried out in the canonical way. For soft 
loops UV and IR divergences are not distinguished and both contribute to the 
anomalous dimensions. The renormahzation group equations are formulated 
in terms of the velocity scahng parameter v. For illustration let us consider 

the NNT j , order running of the color singlet Coulomb potential -^^^ in some 
detailĚj As mentioned previously, the hard matching conditions only contain 
an order Born contribution, Vc(T) = 47ras(?n), Vc(l) — 0. In the foUowing 
discussion it is understood that as runs with the three-loop /3-function. As for 
the matching calculation, the running is determined using on-shell four-quark 
amplitudes. 

2 

At LL order the only effective theory graphs cx -p- are the soft diagrams in 
Fig. Q2c, where the soft gluon vertices are given in Eq. (B3). The result reads 




k2 



(61) 



where (3q — 11 — | and ai — ~ ^ ítí the MS scheme, and is the 
number of massless quarks. The divergence is canceled by a counterterm for 
the operátor , which causes the coefScient Vc to run with the anomalous 

dimension ~2j3Qas{mv), i.e. Vc = '^■nas{mv) at LL order. The remaining 
terms in the soft graphs^axe identical to the one-loop static potential calculation 
of Fischler and BilloireCj This can be understood because the l/m expansion 



of the static calculation in Ref. |43| just separates out the soft region. 

At NLL order the effective theory diagrams cx ^ have iterations of poten- 
tials as shown in Fig. ^ and purely soft diagrams such as shown in Fig. 
The potential diagrams are UV-finite and reproduce the Coulomb singularities 
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Figuře 17: Examples of order af/f diagrams with soft vertices. The vertex with a cross 
dcnotcs an inscrtion of a one-loop counterterm. 




Figuro 18: Graphs with an ultrasoft gluon, which contribute to the three-loop running of 

thc Coulomb potential. 



in full QCD. The soft|-diagrams are identical to the static QCD calculation 
of Peter and SchroderEŽl up to the pinch-singularities caused by the id pre- 
scription of the static quark propagators. In the static QCD calculation the 
pinch-singularities are removed by the exponentiation of the static potentialeS 
(Sec. ^). The effective theory computation does not have pinch-singularities 
because there is no iS prescription in thc soft Feynman rules. After subtraction 
of subdivergences the remaining divergence is canceled by a two-loop counter 



term for 



Up to order aj the counterm reads 



1 



í2 



(62) 



where /3i = 102 — ^ng is the two-loop /3-function, so that vi"^' — A'Kas{mv) 
also at NLL order. At NNLL order, diagrams with ultrasoft gluons have to be 
considered for the first time. In Coulomb gauge we need to consider graphs 
with p.A/to vertices as well as the coupling of ultrasoft gluons to the Coulomb 
potential. After subtraction of subdivergences through counterterms of the po- 
tentials cx l/(m|k|) and 1/m^, the diagrams with ultrasoft gluons that have 
UV divergences not completely canceled are shown in Fig. |l8|a and b. Fig- 
uře shows a counterterm diagram needed to subtract a UV subdivergence 
associated with the l/(TO|k|) potential. Because the energy in the potential 
loop is of order mv^ , the potential and the ultrasoft loops are not separable. 
The sum of the diagrams in the singlet channel reads 



Fig. 18 = 



Ai CfC\ 



a, 



(47r)3 k2 
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AíCfC\ [as{miy)f as{mv'^) fif 



7 -'"(i) 



(63) 



The total anomalous dimension in the singlet channel from thc ultrasoft dia- 
grams with respect to the velocity scahng parameter v reads 



lu 



é£^A [as{mi^)f asimu"^) . 



(64) 



After subtraction of subdivergences the three-loop soft graphs contain a UV 
divergence associated with the three-loop MS /3-function and a UV divergence 
that is induced by the UV divergence in the ultrasoft graphs discussed above. 
So the sum of all soft graphs gives the foUowing contribution to the anomalous 
dimension 



(47r)2 



(65) 



where (52 = - -^fp ni + ^nj . The result in Eq. (|§) has been determined 
in Ref. without an explicit calculation. The íirst term on the RHS of 
Eq. (|65| ) can be obtained directly from Eq. (|63|), because the overall dependence 
on the scales and on as of the corresponding divergence in the sum of the soft 
diagrams is proportional to [as (TOi^)/i|ř]'*/e. Taking into account the matching 



condition for the singlet cc 
the fuU NNLL result reads 



icient at thc hard scale, Vc^\l) — — 47rCi?as(m) 



+ 



íttCfCI 



A „3 



3/3o 



a,(m) 



11 

T 



2z- 



z 

T 



- 41n(í«) 



(66) 



where z = as{rav) j as{m) and w = as{mv^) / as{m) . This result is needed for 
a NNLL order description of the non-relativistic QQ dynamics. At the same 
order the coefHcients of the l/m? potentials are needed at LL order, and the 
coefficients of the l/(m|k|) potentials are needed at NLL. The corresponding 
computations were carried out in Refs. ||5l| ] and [^^ , respectively. 

To describe the production of QQ pairs at NNLL order in e+e~ annihila- 
tion one needs the LL evolution of the coefficients C2 and C3 (sec Eq. (pŤ])) and 
the NNLL order evolution of ci. The complete NNLL running of ci is not yet 
known. The NLL running was determined in Refs. |^,^ and arises from UV 
divergent one- and two-loop graphs with insertions of the l/(m|k|) and Ijm? 
potentials. At LL order ci{v) — 1, i.e. it does not run. The LL evolution of 
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Figuře 19: Graphs with an ultrasoft gluon with p.A couplings which contributc to the 

running of 02(1') ■ 



the axial-vector coeSicient C3 would be determincd by UV divergent diagrams 
of order agV, but such diagrams do not exist and cs^v) = 23(1) = 1. The LL 
running of C2 was determined in Ref. ]37| ] and arises from UV divergences in 
the ultrasoft graphs of order asU^ shown in Fig. The puU-up to the hard 
scale is provided by tadpole-Hke soft one-loop diagrams that vanish in dimen- 
sional regularization. The soft graphs arise from those in Fig. by replacing 
the ultrasoft gluons by soft gluons and contracting the internal quark line to 
a point. The resulting anomalous dimension reads 



d 8Cf , 2\ / \ 

v-^C2{v) = — — as{mu ) ci{v) , 
au Stt 



(67) 



with the solution 



1 8(7f 
6 3/io 



as{m) 



(68) 



where Eq. (p7|) is used as the initial condition. 



6 vNRQCD versus pNRQCD 

The effective theories vNRQCD and pNRQCD are not equivalent. From the 
phenomenological point of view pNRQCD is the by fare more ambitious theory, 
because it was devised with the aim to treat systems with the hierarchy m 3> 
mv 3> mv'^ 3> Aqcd as well as systems, where Aqcd is larger. In addition, it is 
supposed to treat, at the same time_systems with dynamical quark pairs and 
systems where the quarks are staticEj It is fair to say that some conceptual and 
technical aspects are not yet fuUy worked out, as the authors admit at sorofi 
occasions, but it is argued that this does not affect the concepts of pNRQCDĚEl 
On the other hand, vNRQCD has been designed only for the oase m 3> mv 3> 
mv'^ ^ Aqcd , with the primary aim to provide a theory that is fully worked out 
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Figuře 20: One-loop vNRQCD diagrams contributing to the LL anomalous dimension of 

V2. 

technically and conceptually for this speciál case. The treatment of systems 
where Aqcd is larger is by construction difficult in vNRQCD. The theories 
pNRQCD and vNRQCD are not equivalent, because they lead to different 
predictions in the case m ^ mv ^ mv'^ ^ Aqcd ■ 

It is very instructive to have a closer look at the mechanisms working in 
vNRQCD and pNRQCD by means of two examples. One, the order m^a^ lna 
Lamb shift in hydrogen in QED, where both effective theories agree, and the 
other, the Coulomb potential in the QCD Schrodinger equation for a color 
singlet QQ pair, where both theories disagree. For simplicity, I will only con- 
centrate on the mechanism of renormalization and the structure of logarithmic 
terms, which is sufíicient to pinpoint the conceptual and technical differences. 
We note that both theories also disagree on their results for the l/m potentials 
and the QQ production currents. Those cases will not be discussed here, but 
the origin of the disagreement is the same as for the Schrodinger Coulomb 
potential. 



6.1 Hydrogen Lamb Shift 

The Lamb shift is the 2Si/2^'2'Pi/2 level splitting in the hydrogen spectrum. 
For infinite proton mass the order m^a^ In a term reads 

Aí; = -!!^lna, (69) 
3 vr 

where a = 1/137 is the fine structure constant and nie the electron mass. The 
result in jEjq. (|69| ) can be derived in vNRQED from the LL running of the V2 
potentialE^I (Eq. (^)) and in the framework of pNRQED from the LL running 
of the D^^l potential (Eq. (p3|)) (sec Ref. for a fixed order calculation). 

In vNRQED the LL anomalous dimension of V2 conies from soft scattering 
loop diagrams, such as in Fig. po| a and b and from ultrasoft graphs, such as in 
Fig. pOp and d. In QED the leading soft two-fermion vertex is of order l/m and 



there is no contribution from the soft diagram in Fig. 20 a,, if the proton mass 



is infinite. The soft diagram in Fig. |20p gives a massless tapole and vanishes in 
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dimension regularization. Distinguishing between soft UV and IR divergences 
the result reads 



Fig. 20h = -z—[af,ir 



1 1 



(70) 



In Coulomb gauge the ultrasoft contributions are generated by ultrasoft pho- 
tons with p.A/m couplings drcssing the Coulomb potential. Typical diagrams 
are displayed in Figs. EOp and d. For infinite proton niass only the diagram in 
Fig. pO| c contributes. For on-sheU external electrons {E = ^) and including 
wave function renormahzation constants the resuh reads 



Fig. 20c = -i—^[afil'] [a^if] 

óTfí„ 



1 1 



(71) 



The uhrasoft IR divergence in Eq. ( |7l| ) is present also in fuU QED and not 
relevant for the renormahzation proceduře. This can be seen from the fact that 

2 

the IR divergence is regularized for an off-shell electron with E The 
ultrasoft UV divergence is canceled by a counter term in the unrenormalized 
V2 potential. The soft IR divergence is unphysical and does not exist in fuU 
QED. It is associated with the soft scale rrieV and matches with the ultrasoft 
UV divergence, because in vNRQED the photon from fuU QED has been split 
into soft photon fields and an ultrasoft photon field. The full soft tadpole 
contribution puUs the ultrasoft running up to the hard scale. This feature exists 
for any ultrasoft running in vNRQCD. The finál result for the LL anomalous 



dimension of V2 is determined from the ultrasoft UV divergence in Eq. (71) 
and reads 

^|-V2(^)-|a^ (72) 

In practical calculations the UV and IR divergences in soft loops do not need 
to be distinguished, because the puU-up is an automatic mechanism. This 
means that in the case of the Lamb shift the soft diagram would not have to 
be calculated at all. However, the IR divergences in ultrasoft loops have to be 
identified, because they are IR divergences of full QED and do not take part 
in the renormalization of vNRQED. The anomalous dimension in Eq. ( [Ť2| ) can 
be trivially integrated giving 

V2{f) = -^ + ^a'\n,., (73) 

>> In Ref. Q the coefficient V2 is actually called U2- The coefficient is defined with an 
additional minus sign with respect to the convention for Vj^' in vNRQCD of Eq. (^). 
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where V2(l) = — ^ is obtaincd from tree level matching. Interestingly, 
Eqs. (^2|) and ( [Ť3| ) show that there is no infinite series of logarithmic terms 
oc a{ahih')"' . The series terminates because a does not run and because there 
is no non-trivial mixing from other running operatorsĚJ Since the V2 potential 
is the delta function in configuration space the contribution to the 26*1/2^2^1/2 
level splitting can be readily calculated, 



25* 



V2(z^ — a) 



2S 



i*«=2(o)r 



(74) 



where |í'„(0)p = m^a^/(Trn^), which gives the known order mga^lna term 
in Eq. (|69|). A complete calculation of the order niea^lna Lamb shift for 
arbitrary masses in vNRQED can be found in Ref. |59|. (The same work 
also contains the determination of the order mea^ a Lamb shift, the order 
nieO^ In^ a. hyperfine splitting for arbitrary masses, and the order c? lna and 
c? In^ OL corrections to ortho- and para-positronium decay in vNRQED using 
the renormalization group equations.) 

(21 

In pNRQED the determination of D\ ^ is divided into two steps. In 
NRQED (supplemented by the l/m expansion) the running of the Darwin 
coefíicient is carried out from fi = m to fi = iis- The anomalous dimension 
of Cjj is obtained from one-loop vertex corrections and wave function diagrams 
such as displayed in Fig. and has been carried out before in the framework 



Figuře 21: Onc-loop HQET vertex and wave function diagrams that contribute to the 
anomalous dimension of Cjj. The photon in the loop couples with p.A/rríe and the external 

photon with A". 

of HQET0 In these diagrams the photon in the loop couples with p.A/we 
to the electron. For on-shell electrons with momenta p and p' the sum of 
diagrams relevant for the renormalization of cd reads 



Figs. 21 



(P-P') 



M2 



6tt 



1 



1 

eiR 



(75) 



when IR and UV divergences are distinguished. Interestingly, the result seems 
to be equivalent to the soft vNRQCD amplitudě in Eq. (yW, when the photon 
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P 



Figuře 22: 



(2) 

NRQCD diagram contributing to the Born matching condition of thc ^ 
potential in pNRQCD. 



exchange to the proton is included. However, here the IR divergence is the IR 
divergence from fuU QED, because NRQED is supposcd to describe the whole 
momentům space below the scale nie- At this point this observation seems 
to be only a curious fact, but it is at the heart of thc conccptual difFerences 



between vNRQCD and pNRQCD. The UV divergence in Eq. (|75|) is canceled 
by a counter term in the unrenormahzed cd interaction and the anomalous 
dimension of cd reads 



with the solution 



CD(A^)-l-|^ln(^), (77) 

ÓTT Vme/ 

where coime) = 1 is obtained from tree level matching. Thc second step is 

(2) 

carried out in pNRQED. Thc matching condition for the I?^ ^ potential at 
/i = /is is obtained by calculating the NRQED Born scattering diagram in 
Fig. ^ for off-sheU electron momenta (0,p) and (0,p') for the incoming and 
outgoing electron, respectively. The little square indicates the c^j-coupling. 
This gives 

í)g(. = ..) = |cM..) = |-í^.n(£). (7«) 

The running of D^^l in pNRQCD is determined as described in Sec. ^ from 
the ultrasoft one-loop diagrams in Fig. ^ in the l/m expansion. For electrons 
with external momenta (0;p) and (0,p') the result reads 



2a 



3to2 



1 1 



Fig- 6b = -z—, [a M^^] — — ■ (79) 



The IR divergence is again the fuU QED IR divergence. The UV divergence is 
canceled by a counterterm in the unrenormahzed D\ ^ potential and leads to 
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(2) 

thc foUowing LL anomalous dimcnsion for ' 



ÍÍ,S' 



This is again integrated trivially giving 



= fcz5(Mn) (81) 

at the scale /i = For /i„ = mea^ this agrees with V2 up to a convention- 
dependent factor and also leads to Eq. (|6^ ) for the order mga^ In a Lanib shift. 
Comparing the computations we see that the difFerence between vNRQED and 
pNRQED appears to be simply that the running is distributed in a different 
way. In vNRQED running happens in one step and in pNRQED it happens in 
two steps. Agreement is in fact found for the LL running of all l/m? potentials 
in vNRQCD (Refs. |l|,|o|) and pNRQCD (Ref. [||). 

One of the main differences between QED and QCD is that the coupHng 
constant runs itself . If we consider the Lamb shift of muonic hydrogen and treat 
the electron as massless we can "simulate" , in a simphfied way, one difFerence 
of the calculation presented above to the QCD oase. If massless electrons are 
included, the fine structure constant becomes scale-dependent using leťs say 
the MS renomialization scheme. This makes the couplings in all potentials 
scale-dependent and makes the 2S'i/2^2Pi/2 level splitting more complicated. 



As shown in Sec. 6.2, the a in the Coulomb-potential has to be set at the scale 
/Í3 or mv and, as the consequence, also the a in the wavefunction at the origin 
in Eq. ( [7^) . In addition, also the a in the Spin-Orbit potential (the D^^^g ^ term 



in Eq. (23) and the Va terms in Eq. (M)), has to be set at the scale /is or mv 



and now contributes logarithmic terms to the 2<S'i/2-2Pi/2 level splitting. For 

(2) 

the purpose of this discussion we only consider the coefRcients ^ and V2 in 
some detail, because they contain a non-trivial ultrasoft running involving the 
scales /i„ or mv^ . In the pNRQCD computation for cd thc coupling in the 
anomalous dimcnsion in Eq. ( [76| ) is a (/i) and the result for c^j is 

e.(M)-l + :^ln(^), (82) 

where j3o = — ni = 1. The a in the matching condition of Eq. ( [ŤS] ) is 
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(2) 

evaluated sX fi — fis and thc LL anomalous dimcnsion for Z?^ ^ reads 

d {2), . 4:a{fis) , . , . 

= (83) 

The coefHcient of the LL static potential is set at the scale fis and does not 
run below fis- The solution just reads 

In vNRQED with massless electrons the anomalous dimension for V2 also gets 
a contribution from loops in Fig. ^ a and b where the soft photons are replaced 
by soft (massless) electrons. The modified version of Eq. (fŤ^ ) then reads 

^^-^V2(!^) = ^ a{m^v) a{m^v'^) + ^ [a{mf_,v)Y Ci^,2e{v) , 

Cifi2e{v) = coim^v'^) , (84) 



where C4^2e is the coefficicnt of the 6-fermion operátor contained in Fig. |2C 
(This type of operators, involving four heavy quarks and two soft gluons or 
masslfiss quarks, has been missed in Ref. jsi] ]. Their numerical contribution is 
smalltHl and does not affect conclusions drawn based on the results of Ref. [pT| .) 

It turns out that this coefficient is equal to the co obtained in NRQED at 
the scale /i = rrif^v^. Integrating Eq. ( |8^ ) with the initial condition V2(l) = 
— 7ra(m)/2 gives 

V2{ry) = CD{my) , (85) 

which agrees again with the pNRQED result up to a convention-dependent fac- 
tor. The vNRQCD and pNRQCD computations agree, because only one-loop 
diagrams contribute to the LL anomalous dimensions of the 1/m^ potentials. 
So at this level it seems to be merely a choice of convention whether one uses 
vNRQCD or pNRQCD, because the running is just distributed in a different 
way. The actual reason for agreement, however, it that at LL order the (one- 
loop) diagrams cannot contain non-separable ultrasoft and potential loops at 
the same time, which connects in a subtle way the evohition with respect to 
soft and ultrasoft scales. 



6.2 Schródinger Coulomb Potential at NNLL Order 



The subtle difference in the renormalization of vNRQCD and pNRQCD be- 
comes apparent beyond the LL approximation, when non-trivial multi-loop 



Heavy Quarkonium Dynamics 



56 



diagrams with UV divergences contain loops dominated by the ultrasoft and 
the potential region. Examples are the Coulomb potential oc l/k^ in the 



NNLL order Schrodinger equation for a color singlet QQ pair (see Ref. [55 
for vNRQCD and Ref. [|§ for pNRQCD) or the l/(™|k|) potential at NLL 
order (see Ref. ||] for vNRQCD and Refs. ||@ for pNRQCD). Because 
the l/(TO|k|) potentials contribute to the NLL running of the QQ produc- 
tion current, vNRQCD and pNRQCD predictions also differ for the latter case 
(see Refs. [||,|o| for vNRQCD and Ref. @ for pNRQCD). Let us discuss 
the NNLL order Coufomb potential that appears in the Schrodinger equation. 
Its determination in vNRQCD and pNRQCD has already been deseribed in 
previous sections and will not be repeated in detail. 

In vNRQCD the NNLL l/(p — p')^ potential that appears in the Schrodin- 
ger equation, called in the foUowing, consists of the sum of the Coulomb 
potential Vc(s)/(p — p')^ at NNLL order and the time-ordered products of soft 
vertices proportional to l/(p — p')^ at one (Figs. [T^ c and Eq. (|6l|)) and two 
loops (see Figs. that determine the running of Vc at LL and NLL order. 
Im momentům space for Z) = 4 the sum is (k = (P ^ p')) 



K"(p,q) = 

'as(Aí, 



+ 



k2 



47r 



/3o In 



Ol 



47r 



Eq. (26) 



^■kCfC\ 3 



3/3ok2 



af(m) 



2z h 41n(w) 

4 2 ^ ^ ' 



(86) 



where z = as(mv) j as{ra) and w = as{mv'^)/as{mv). In principle it is in- 
appropriate to consider the soft time-ordered products together with the Vc 
term as a single quantity for D = A because they are matrix elements depend- 
ing non-trivially on D — A — 2e dimensions. Therefore, they should not be 
considered for e ^ O independently of the other matrix elements relevant for a 
given process. However, as long as insertions of the soft time-ordered products 
do not generate UV divergences relevant for the process, Eq. (|8^) can be used. 
This is also the case for the NNLL order description of QQ pairs, when the 
quarks are stable. 

In the pNRQCD framework the NNLL l/(p - p')^ potential that appears 
in the Schrodinger equation, called in the foUowing, is obtained directly 
from the two-step matching and running through NRQCD and pNRQCD. The 
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NNLL result reads 



K^(p,q) 



k2 

The first term is the two-loop pNRQCD matching condition and the second 
term arises from the pNRQCD running from to AU computations are 
carried out in the static hmit abandoning the non-rclativistic power counting 
for QQ pairs because otherwise the NRQCD computations cannot bc properly 
dcfincd. It is one of the most important assumptions in the pNRQCD frame- 
work that the static calculation of the running coefficients can be apphed also 
for moving dynamical quark pairs. 

Comparing Eqs. and (^Ťj) one finds disagreement in the contributions 



associated with the ultrasoft gluons. Expanding the second term in Eq. (86) 
(times a factor k^) in as{m), we find 

- i Cpď^ atim) Hu) + g Cpď^ aUm) ln'{,.) + .... (88) 

The analogous resuh for the static calculation in Eq. (^) for = miy and 
Mu — miy^ reads 

- l CfCI atim) Hu) + ^ CpCl al{m) \n\v) + .... (89) 

The single In u terms agree, but the higher order logarithms differ. From the 
phenomenological point of view the difference is quite small for most applica- 
tions, but both theories claim to sum the correct higher order logarithms. 

In vNRQCD the result arises from the UV divergences in the three-loop 
ultrasoft-potential diagrams in Fig. |l^ and the associated UV divergences in 
three-loop soft diagrams. The individual quark propagators i/{kQ~'k'^ /m-\-i5) 
in the ultrasoft-potential diagrams describe dynamical quarks. The resulting 
anomalous dimension in Eq. (|6J) describes running from v — Ito v = v and in- 
cludes the scale-dependence from all three loops, i.e. the factors /z^ř for the 
ultrasoft-potential diagrams and fi^g for the soft as well as the scale-dependence 
of all the coupling constants that occur in the diagrams. In pNRQCD the re- 
sult arises from the UV divergence of the one-loop diagram Fig. ^ with one 
ultrasoft gluon and a static QQ propagátor i/{idQ + {^Ca — CF)oíVa/f)- This 
propagátor sums an infinitc number of static potential rungs in Z? = 4 dimcn- 
sions. The resulting anomalous dimension in Eq. (|30|) describes running from 
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Figuře 23: Loop graph with an insertion of the l/(m|k) and l/k^ potentials. 



ji = jis mv to /Lí = /Lítí ~ mv'^ and includes only the scale dependence of 
the ultrasoft loop, i.e. the factor /x^*^ and the ultrasoft coupUng as{fi)- The 
quarks in the loops with the Coulomb rungs are not dynamical. Clearly, if 
the calculation would involve only true one-loop diagrams, as was the case for 
the Lamb shift discussed before, we would find again agreement, because loops 
involving multiple iterations of potentials wouldn't exist. 

At the conceptual level the difference can be understood as the difference 
between a calculation for static and dynamical quarks. In the dynamical vN- 
RQCD calculation the correlation between soft and ultrasoft scales through the 
quark equation of motion is maintained at all times. In the static pNRQCD 
calculation the scale correlation is broken at intermediate steps through the 
l/m expansion, i.e. energy and three- momentům are treated independently. 
The discussion above shows that dynamical and static QQ systems are differ- 
ent. In particular, it can be problematic to use results obtained in static quark 
systems for problems involving dynamical quarks. 

AjL-jthis point it is, once more, instructive to see with a very simple ex- 
amplcĚa how the l/m expansion breaks the non-relativistic power counting. 
In the static calculation the time-ordered product of a l/(TO|k|) and a 1/k^ 
potcntial is 



ď^k 1 1 111 



(27r)4 m|k| (fco + iS) (fco - i6) k^ m 



(90) 



and hence suppressed by l/m. In the dynamical calculation (see Fig. |2^) the 
time-ordered product of a l/(m|k|) and a 1/k^ potential is 

d^k 1 1 11 



i2nr m\k\ {ko-^+^S) (fco + ^ - .<5) 

d^k 1 TO 1 , , 

^ 1 (91) 



(27r)3 TO|k| (k2 - iS) k2 



and of order 1. 
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7 Heavy Quark Mass and Renormalons 



In the discussions of the previous sections on non-relativistic effective theories 
we have assumed the pole mass definition m for the heavy quarks. From 
the effective Lagrangians one can derive the one-particlc-irreducible two-point 
function for a quark with energy E = ^/š — SíTip^^j^, being the total center- 
of-mass energy, and three-momentum p, 

E~-^ + -^ + ..., (92) 

which vanishes for an on-sheU quark in the centcr-of-mass frame (_B = p = 0). 
The on-shell condition is useful for the determination of matching conditions 
and anomalous dimensions in the effective theories. So at this stage the use of 
the pole mass definition is appropriate and useful. For a generál mass definition 
M = JTipoic + SAI the effective Lagrangians contain additional bilinear quark 
operators SMip^ip, SMip^ -^ip, etc, and the one-particle-irreducible two-point 
function for a quark reads 

E,, + 2SM ~ + + ^ + ... , (93) 

where \/š - 2M . 



7.1 Pole Mass and Static Potential 

Since quark masses, just like coupling constants such as cts, are parameters 
of the Lagrangian and not physical observables, in principle any quark mass 
definition can be employed as long as its perturbative definition is known and 
well defined. A suitable mass definition should also have reasonable properties 
such as being gaugejjivariant and IR-safe. The pole mass definition certainly 
belongs to this classĚS However, there are good reasons not to use the pole mass 
definition for the analysis of experimental data, because it induces artiíicially 
large perturbative corrections. In the framework of the thcory of asymptotic 
series these perturbative corrections can, with certain assumptions on their 
large order behavior, be interpreted as an intrinsic ambiguity of the pole mass 
definition of the order Aqcd ~ 300 MeV caused by a strong linear sensitivity 
to small momenta of order Aqcd- (The exact amount of the uncertainty is 
diíHcult to estimate, and I have just given a personál (and probably biased) 
number to be definite. It is also not intended to rely on a specific scheme when 
the notion of "Aqcd" is used.) This unpleasant property is, at least from the 
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perturbative point of vicw, not characteristic for all quark mass definitions. 
For the pole mass this property wasJicst discussed in the analysis of B meson 
decays in the framework of HQETEÍ'L3 In the context of non-relativistic QQ 
Systems the same problém was first realized by Aglietti and LigetiE3 in an 
analysis of the large order behavior of the perturbative static QQ potential. 
Initially, however, the connection to the properties of the pole mass parameters 
was not noticed. The problems associated with the pole mass definition are 
not academie. For heavy quarkonium systems it is important to consider the 
issue seriously, since the parametric perturbative uncertainty of present next- 
to-next-to-leading order computations is already at a level where an ambiguity 
of a mass definition of prder Aqcd is relevant. 

Aglietti and LigetiE3 considered the color singlet static potential in the 
"large-/3o" approximation, 



''static 



OO 

E 

n=0 



/3o In 



(94) 



where c ~ —5/3 in the MS scheme. The large-/?o approximation assumes 
that the highest power of /3o in each orden-of perturbation theory dominates 
numerically the rest of the contributionsĚII At this point the distinction be- 
tween the static pNRQCD potential and the Schrodingcr Coulomb potential 
in vNRQCD is irrelevant, because the large-/3o approximation is the same in 
both cases. The large order behavior can then be analyzed through the Borel 



transform of Vj,^"^.;^. with respect to as/3o/(47r), 



OO 



/3ok2 
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/kV\ 
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The inverse Borel transform reads 



(95) 



Po 



(k) = r duB[Vtt]{u) exp ( - 



The impact of the large order behavior of Eq. (94) in bound statě calculations 
can be visualized by considering the static potential in configuration space, 
describing the static energy of a QQ pair with distance r <C in the pole 

mass scheme, 



^tóticK 
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stotic] (^) 
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f3or 
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r(^ + ^)r(^- 
r(2M + 1) 



(96) 



The RHS of Eq. ( P6[ ) has poles on the positive real axis at w = (2fc + l)/2, 
fc = 0, 1, 2, 3, . . ., with residues oc {firy^'^^ /r. These poles cause ambiguities 



in V£l^^{r) proportional to 



(/ir) 



2fc+l 



■ exp 



47r 2fc + 1 
/3oas(Aí) 2 



(Aqcd ^ 



\2fe+l 



since there is no a priori prescription how the poles on the positive real axis 
should be treated in the inverse Borel transformation. The pole closest to the 
origin at u = 1/2 is dominant and causes an ambiguity of order Aqcd- K 
corresponds to an asymptotic large order behavior of V!,f ÍJ^. (r) proportional to 



2ti 



One might wonder how such an r-independent behavior can arise in the static 
potential, which has an overall factor 1/r. It was shown in Ref. |68| that it 
arises from an exponential structure of the ln(/ir) terms in the order a^^^ term 



in Vgfatí'^) the form 
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In" "^(^r) + ... + - In^(^r) + In(^r) + 1 



r \ 2tt ^ \ 2tt 



(97) 



It should be noted that the large-/3o approximation is just a heuristic concept. 
However, explicit multi-loop calculatinná for example for the perturbative re- 
lation between MS and pole masg£j'L2rO have shown a remarkable consistency 
with th£j|irge-/3o approximation. On the other hand, it has also been re- 
portedOO that the large-/3o approximation gives only order of magnitude es- 
timates for higher order corrections in mass relations that do not involve the 
pole mass definition. Numerical evidence for the bad behavior of the static po- 
tential in configuration space was in fact noticed quite early, see e.g. Refs. [||, 
but the observations were not interpreted as a consequence of the pole mass 
definition. 

If the previously described large order behavior would be real, the deter- 
mination of perturbative quark mass parameters would be forever limited to 
an uncertainty of order Aqcd- However, it was found by Hoang et al.E3 and 
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Beneketj (see also Refs. |7^]) that the ambiguity of order Aqcd is unphysical 
and merely an artifact of using the pole mass definition. The pole mass coun- 
terterm subtracts, apart from the usual UV divergent term, also a finite piece 
that correspond to the static quark selfenergy, 

-Jj^y^M^)- (98) 

This expression is just the leading term in the non-relativistic expansion of the 
finite piece in the pole mass counterterm. Its asymptotic large order behavior 
is exactly 1/2 times the asymptotic large order behavior of the static potential 
displayed in Eq. (|9^) . If a mass definition is chosen that does not subtract the 
static selfenergy, the order Aqcd ambiguity and the corresponding bad large 
order behavior of the peptiirbative seriea^described above do not arise in the 
first pláce. Hoang et alEÉ and BenekeLj demonstrated that the total static 
energy, 

Estat = 2mpolc + Kítatír) , (99) 

is free of the ambiguity of order Aqcd, at least from the point of view of 
perturbation theory. 

It was shown by Smith and WillenbrockEa that even if the quark is un- 
stable and decays with a finite lifetime smaller than 1/Aqcd, the static quark 
selfenergy in Eq. ( ^Šl ) produces the same large corrections described above. 
Thus the problematic behavior of the pole mass definition remains also for 
short-lived quarks such as the top quark with Tt ~ 1.5 GeV in the Standard 
Model. 

The previous discussion disqualifies the pole mass a priori for the use in 
analyses of experimental data, since it induces artificially large corrections that 
are compensated e.g. by shifts in the pole mass value itself, when higher order 
corrections are included. In principle, the pole mass can still be employed 
as an order- and smle-dependent correlated quantity, similar to the value of 
a matrix elementJIj In such an approach the large perturbative corrections 
associated with the pole mass are containcd in its numerical value. Since these 
corrections depend on the order and on the choicc for renormalization scales 
and couplings, the numerical value of the pole mass needs to be treated as 
a function of these parameters in order to achieve the proper cancellation of 
the large corrections, if the pole mass value is used in computations for other 
mass-dependent quantities. It should be noted, however, that this stratégy can 
become increasingly unreliable for orders n where the corrections induced by 
the pole mass are large enough that their numerical cancellation is incomplete. 



Heavy Quarkonium Dynamics 



63 



Results for analyses in the pole mass scheme for the tt production close to 
threshold in e~^e^ annihilation and for bottom mass extraction from bb sum 
rules are discussed in Secs. ^ and |o|. 

The subdominant anibiguity in the static potential is order AQ(^j-,r ~ 
Aq^^/ (mas) and quite smaU parametrically. The associated large order be- 
havior of the perturbative series was analyzed by Brambilla et al. in Ref. , 
and it was found that it cancels with the large order behavior of the leading 
ultrasoft corrections shown in Fig. ^ coming from two insertions of the r.E 
interaction. The cancellation is independent of the choice for the quark mass 
definition. The ambiguity of order Aq^qT can be ignored for any practical 
analysis of the non-relativistic QQ dynamics. 

7.2 MS Mass and Threshold Masses 

To avoid the problems of the pole mass definition it is advantageous to use 
quark mass definitions that are less sensitive to low momenta and do not 
contain the static quark selfenergy in the mass counter term. Such quark mass 
definitions are called "short-distance" masses and have a parametric ambiguity 
or order Aqq^/M or smaller. In principle, the most natural choice of a short- 
distance mass would be the MS (or MS) mass, because its mass counter term 
subtracts only the UV poles oc 1 /e" and does not contain any infrared-sensitive 
subtraction. This means that the MS mass is only sensitive to scales of order 
or larger than m and that it is from the conceptual point of view not allowed 
to lower the scale of the MS mass below m. The relation between the pole 
mass and the MS mass m(m) reads (a^ = ai"*^^'' (m(m))) 

Wpoic , , 4as _ (a 



. . , . , , 13.4434 - 1.0414n,l 

m[m) óTT V TT 

^ ^190.595 - 26.655n<> + 0.6527n^ ) + . . . (100) 



for n£ massless light quark species. The two-loop correction was determined 
in Ref. and the three-loop term in Ref. (see Ref. |^l[] for a numerical 
computation) . For b quarks with ág = 0.22 and ng = 4 the series on the RHS of 



Eq. (lOC) reads 1 — 0.093 + 0.045 + 0.032. Its rather poor convergence properties 
at already low orders illustrate the infrared sensitivity of the pole mass. The use 
of the MS mass for non-relativistic QQ systems would in principle be welcome, 
because it is also the preferred mass definition for high energy processes or 
electroweak precision observables. Unfortunately, the MS mass breaks the 
non-relativistic power counting. The parametric size of the terms of the one- 
particle-irreducible two-point function of Eq. (621) is of order mv^ ~ mal . In 
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the MS scheme, however, Sm-^ is of order mas and exceeds all the other terms 
in Eq. (p^), which would in principle require that Smy^ is not tijeated as a 
correction. This leads essentially back to the pole mass definitionEj Explicit 
examples for the use of the MS mass for tt production close to threshold in 
e~^e~ annihilation and for bottom mass extractions from bb sum rules have 
been given in Ref. From the physical point of view, the incompatibility 

of the MS mass with non-relativistic QQ systems is a consequence of the fact 
that it is a mass definition adapted to the situation where the heavy quark has 
a virtuality |q^ — of order or larger than m^, q being the heavy quark four- 
momentum. Thus the MS mass is the preferred mass definition to describe 
virtual effects of a heavy particle in low energy processes with characteristic 
energies E <^ m ov mass effects in high energy processes with characteristic 
energies E ^ m. 

An appropriate mass definition to use for non-relativistic QQ systems is a 
short-distance mass that also complies with the non-relativistic power counting, 
i.e. 5M is parametrically of order Mv'^ ~ Mol. Such masses are called "thresh- 
old masses" , since they are adapted to the situation where the quark virtuality 
is smalltil Obviously an arbitrary number of such threshold masses can be in- 
vented, because the number of subtractions possible to achieve 5M ~ Ma^ is 
infinite. Meanwhile, a considerable number of threshold mass definitions can 
be found in the literatuře. In the foUowing I briefly comment on the threshold 
masses in the current literatuře in historie order. 



Kinetic Mass 

The kinetic mass has been proposed by Bigi et al. in Ref. |^^. Originally it 
was devised as a mass definition for the description of B mesons. It is defined 
as 



pert 



kin-' 



277lpolc 



pcrt 



_ 16 as 



(101) 



where [^(Aíkin)] port [^^(^kin)] pert ^""^ pcrturbativc evaluations of HQET 
matrix elements that describe the difference between the pole and the B meson 
mass. The two-loop contributions to the kinetic mass have been calculated in 
Ref. 1^^. The three-loop contributions, which correspond to NNLO in the 
non-relativistic power counting, are known in the large-/3o approximation. In 
the first line of Eq. (101) the ellipses indicate matrix elements of operators 
with higher dimension, which have not been included in any analysis so far. 
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The kinctic mass is dependcnt on thc scale which is thc momentům eut- 
off for the evaluation of the matrix elements. The cutofF has to be chosen of 
order mag to comply with the non-relativistic power counting. Thus the order 
Mkiii'^s term is LO in the non-relativistic power counting, the order A<kin'^s term 
is NLO, etc. If the kinetic mass is used in the perturbative series for quanti- 
ties that do not have non-relativistic power counting, the explicit counting in 
powers of as is used, i.e. /ij^j^^ is formally counted of order m. This is the case 
for example in the relation between the MS and the kinetic mass or when the 
kinetic mass is used in computations for B mesons. 

Potential Subtracted Mass 

The potential subtracted (PS) mass mps has been proposed by BenekeEl and 
is defined by 

|q|<Mps 
1 f d h ~ 
AÍps(mps) = mpoic +2 J J2ň)^ V;tat(k) 

4 a 

= TOpoic - o ~ '^PS + ■ • ■ . (102) 
where V^tat is the singlet static potential in momentům space representation 



(see Eqs. ( |26| ) and (pT\))- In the last hne of Eq. ( |102| ) the first term in an 
expansion in as is displayed. The additional subtraction of the PS mass is 
equal to the non-relativistic selfenergy mentioned before and can be regarded 
as a kind of minimal way to cancel the large higher order corrections in the 
static potential in Eq. (99). This subtraction has in fact been proposed before 
by Bigi et al. in Ref. |64], but no explicit mass definition resulted from the 
subtraction prescription. The PS mass is dependcnt on the scale /ipg, which 
is the cutoff for the selfenergy integration. The cutoff has to be chosen of or- 
der mas to comply with the non-relativistic power counting. Thus the order 
/ipgOís term is LO in the non-relativistic power counting, the order /ipga^ term 
is NLO, etc. At NNLL order in perturbation theory thc PS mass receives an 
additional dependence on the pNRQCD renormalization scale /i„ due to the 
ultrasoft anomalous dimcnsion of the static potential caused by the selfenergy 
diagram shown in Fig. |^. This scale is correlated to the renormalization scale 
/ip used in the computation of the pNRQCD matrix elements by Eq. (^2|). 
If the PS mass is used in perturbative series for quantities that do not have 
non-relativistic power counting, the explicit counting in powers of is used, 
i.e. /ipg is formally counted of order of m. A modified versiťpn of the PS mass, 
called the PS mass was proposed by Yakovlev and Grootelá In the PS scheme 
additional 1/m-suppressed corrections from the non-relativistic expansion of 
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the quark selfenergy are included in Eq. (102) 



IS Mass 

The IS mass has been proposed by Hoang et aL and is defined as half of 
the perturbative series for the mass of the n — 1, "^^^^Lj ^ ^Si quarkonium 
bound State, 



1 



pcrt 



(103) 



In the last hne of Eq. (103) the LO term is displayed. The IS mass is scale- 
independent, and it is determined with the same methods that are used for 
the computations of the non-relativistic quantities for which it is used as the 
quark mass definition. By construction, the value of the IS mass very weakly 
correlated to other parameters such as Us or the renormalization scale in ex- 
tractions from the QQ spectrum. As such it will in generál have smaU pertur- 
bative unccrtainties. If the IS mass is used in perturbative series for quanti- 
ties that do not have non-relativistic power counting, the order term has 
to be counted as order a^, the order term has to be coujited as order 
a^, etc. This prescription is called the "upsilon expansion"E3 In the case 
m 3> mv 3> rav^ ^ Aqcd the physical l^Si quarkonium mass is equal to 
2Afis up to non-perturbative corrections, which are parametrically of order 
mKQQY)/ (ma)^ and which can be related to local condensates of the operators 
product expansion of Shifman et al.a (Sec. [ll| ). 



Renormalon Subtracted Mass 

The renormalon subtracted mass has been proposed by PinedaS and is defined 
as the perturbative series that results from subtracting all non-analytic pole 
terms from the Borel transform of the pole-MS mass relation at u = 1/2 with a 
fixed choice for the renormalization scale = A*j^g. The scale //j^g is then kept 
independent from the renormalization scale used for the computation of the 
non-relativistic quantities of interest. The terms in the relation between the 
pole mass and the RS mass are formally known to all orders, but the numerical 
value of the individual coefficients of the series are known only approximately 
due to yet undetermined subleading contributions to the residue at w = 1/2. 
To order the relation between RS mass and pole mass reads, 



(104) 
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where the constant c depends on the number of light quark species and has 
an uncertainty because the residue at m = 1/2 in the Borel transform of the 
pole-MS mass relation is known only approximately. The scale /j,j^g has to 
be chosen of order mas to comply with the non-relativistic power counting. 
If the RS mass is used in perturbative series for quantities that do not have 
non-relativistic power counting, /Xj^g is formally counted of order of m. Pineda 
also proposed a shghtly modified version of the RS scheme, which he called 
RS' scheme. 

At this point I refrain from a generál numerical comparison of the various 
threshold masses. Comparisons for some cases can be found in the references 
given above (see also Ref. fŠ^). To allow for a comparison between the results 
that have been obtained for the various threshold masses in the literatuře it 
has become practice to determine in a second step the MS mass nT{m). The 
perturbative uncertainty in the MS mass that is obtained from the conversion 
can be larger than the perturbative uncertainty in the threshold masses, be- 
cause the large order as term in the pole-MS mass relation (Eq. ( |100[ )) is quite 
sensitive to the uncertainty of ag ■ 

7.3 Threshold Masses and Heavy Quark Decay 

The application of threshold masses is not restricted to non-relativistic QQ 
systems. Threshold masses are in generál useful for systems, where the virtu- 
ality of the heavy quarks is small, i.e. if the heavy quarks are close to their 
mass-shell. For non-relativistic QQ systems the virtuality — rr? is of order 
rn}v^ <C rri} . A similar situation arises in B mesons, where the heavy quark 
virtuality is of order toAqcd ^ vr? . 

As an example consider the inclusive semileptonic B X^ev decay rate, 
which is relevant for the determination of the CKM matrix element \Vub\- The 
incliisive semileptonic partial rateJiito light hadrons has been analyzed in the 
1SE3 and the kinetic mass schemeO (see also Refs. [^,0). The semileptonic 
B — !■ X^ev partial rate has a particularly strong dependence on the bottom 
quark mass oc G\{m^)^ and its perturbative series illustrates the properties 
of the various mass definitions. In the foUowing I show the behavior of the 
perturbative series of the semileptonic B — > X^ev partial rate in the pole, MS, 
IS and kinetic mass schemes. For simplicity, as = 0.22 at the scale of the 
respective mass definition has been used as the expansion parameter. For a 
realistic phenomenological application of the results in the IS and kinetic mass 
schemes including the effects of non-perturbative contributions and I refer to 
Ref. H. 
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In the pole mass scheme the perturbative series for thc inclusive decay rate 
reads 



r (Š -> X, 



1927r3 



1 ~ 0.17 - 0.10 - 



, (105) 



where Gf is the Fermi constant, andJJae second and third term in the brackets 
correspond to the one- and two-loop£j corrections. The contributions of non- 
perturbative matrix elements are not displayed. In the large-/3o approximation 
the two-loop term is estimated as —0.13, which agrees reasonably well with the 
true result —0.10; the three-loop term is unknown and estimated as —0.12 in 
the large-/3o approximation. The convergence of the series is not very good. In 
particular, the estimate for the three-loop correction indicates that two- and 
three-loop corrections are of the same size. This is an artifact of the pole mass 
definition. If the MS scheme is employed the series reads | 



T(B Xuev) = 



Vub 



1927r3 



(m^(m'')y 



1 + 0.30 + 0.13 



(106) 



In the large-/9o approximation the two-loop term is estimated as -1-0.19 and the 
three-loop term is estimated as -1-0.05. While the use of the MS mass leads to 
a cancellation of the bad pole mass behavior at high orders, the corrections are 
still large at low orders. On the other hand, in the IS mass scheme the series 
for the inclusive semileptonic B X^ev decay rate reads 



Gl\Vub 
1927r3 



(mís)' 



1 - 0.115 - 0.031 - 



, (107) 



where the upsilon expansionE3 has to be used to obtain the terms in the series. 
In the large-/3o approximation the two-loop term is estimated as —0.035, which 
agrees well with the exact result —0.031; the three-loop term is estimated as 
-1-0.005 in the large-/3o approximation. So the series in the IS scheme has 
considerably better convergence properties than in the pole and theJVIS scheme. 
In the kinetic mass scheme for /ij^j^^ = 1 GeV the decay rate readsĚfl 



T B^X„,ei>) - 



Gl\Vub\ 
1927r3 



^"kin 



1 - 0.022 - 0.006 



(108) 



The series converges even better than in the IS scheme, which is a consequence 
of the specific choice ^j^j^^ = 1 GeV. In Ref. all inclusive semileptonic B 
decay modes were studied in the IS mass scheme including an analysis of 



The order contributions in Eqs. {tod) and (107) have been takcn from Rcf. ]9l[ . 
where the new results from Ref. [po| were included. 
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the semileptonic B decay forni factors. Othcr applications of thrcshold mass 
schemes in B decays can be found e.g. in Refs. p^. 



8 Heavy Quark Pair Production at Threshold 

The production of heavy quark-antiquark pairs in the kinematic regime close 
to threshold ^/š ~ 2m is currently the most important apphcation of non- 
relativistic effective theories for quark-antiquark systems in QCD. Moments 
of the total cross section cr(e+e^ bb + X) in the threshold region are an 
important tool for precise determinations of the bottom quark mass param- 
eter, and ti production close to threshold is a major part of the top physics 
program of the future Linear CoUider with the aim to provide high precision 
measurements of the top mass and other parameters. In this section I will 
review the calculational steps necessary for a computation of the total produc- 
tion cross section of a genuine color singlet QQ pair close to threshold assuming 
the hierarchy m ^ mv ^ mv^ 3> Aqcd and concentrating mainly on e~^e~ 
annihilation. For simplicity the quarks are assumed stable and electroweak 
effects are neglected except for the production mechanism. Non-perturbative 
effects are also neglected. Applications to tt and bb production including dis- 
cussions on non-perturbative effects (to the extent they are known) are given 
in Secs. |9| and F or an application on r pair production close to threshold 
in QED see Ref. |93|] . I will hrst review a recent computation of the total cross 
section at NNLL order in the framework of vNRQCD and then comment on 
computations in NRQCD at NNLO in the fixed order approach, i.c. without 
summation of QCD logarithms of v. 

8.1 QQ Production At NNLL Order 

Generically, the normalized total heavy QQ cross section in e+e^ annihilation 
at NNLL order in the non-relativistic expansion takés the form 

R=-i-^^-^QQ) = .y(^)Y(a.in.r 



x<{ l{LL);as, v{NLL);al, asv, (NNLL)^ , (109) 



where the indicated terms are LL, NLL and NNLL order. The summation 
of logarithms can be performed using the renormalization group equations in 
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the framework of the non-relativistic efFective theories discussed in Secs. |^ and 
^. In the foUowing I will review the computation of a{e^e^ QQ) in 1 1£ 
framework of vNRQCD using dimension regularization in the MS schemeEZI 
In the framework of pNRQCD such a computation has not yet been carried 
out, but foUows essentially the same lineš. For QQ production with a different 
initial State (77, etc.) or production of other colored massive particles such as 
squarks or for QED systems the same method can be apphed with straightfor- 
ward modifications. 



Schródinger Equation and Potentials 

At NNLL order ultrasoft corrections enter only through mixing into the co- 
efficients of potentials and currents (Sec. ^). Therefore the NNLL order non- 
relativistic QQ dynamics is described by a common time-independent two- 
body Schródinger equation. In momentům space and with a threshold mass 
definition M (Sec. ^ the NNLL order Schródinger equation takés the form 



SM 



I 1 + r::! ) _ _íl^^ -~(e + 2 SM) 



G{p,p',E) 



- J i?"k í/(p, k) G(k, p',E)= (2^)" (p - p') 



(110) 



where 



E = VŠ-2M = Mv^ 



3-2e , and L>"k = e"^''^' 



In 4tt' 



f^l^ d"k/(2^)" . 



The term 



SM = Af - m„ 



5Ah 



- 6M^ 



'pole — 'JJ^LL "^"NLL I ''-'"NNLL 

is the difference between pole and threshold mass. For the term pi^SM/KP only 
the LL expression needs to be taken into account. The Schródinger equation in 



Eq. (IIC) is obtained from the NNLL order equation of motion of a QQQQ four 



point function with center-of-mass momenta (-j zbpojip) and (-j ±pQ,±p') 
for the incoming and outgoing quarks, respectively, supplemented by a contour 
integration over po a-nd pf, (see e.g. Ref. [Q). The Green function G(p, p', E) 
describes production and annihilation of an off-shell QQ pair with relative 
momenta 2p and 2p', respectively, at total center-of-mass energy ^ = 2M+E. 
The potential V arises from the potential-type operators up to order 1 /m? (see 
Eq. (^) and Fig. [TTIb) and from time-ordered products of the two-quark soft 
interactions (see Eq. ( ^ and Figs. |l^, |l2|c). The fuU potential at NNLL order 
for a QQ pair in a S-wave statě reads 

V{p,p') = Vc{p,p') + Vs{p,p') + Vrip,p') + Vk{p,p') , (111) 
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where Vc has been given before in Eq. (p6|) and (q = p — p') 



^2 



V;(P,P') = ^^j^^V(^)(.). (112) 

The color singlet components of the potential are obtained from E q. (|4^) . The 
running of the potential coefficients was determined in Refs. [ ^l|j5^j55| , |60| . For 
v ^ v a\l QCD logarithms of v that arise from the non-relativistic QQ dynam- 
ics at NNLL order are summed into the coefhcients of the potentials. In e+e~ 
annihilation the QQ pair is predominantly produced in a statě where the 
spin functions in Eq. ( ^ ) are = 2 and A — T — 0. For a ^Sq statě (e.g. 
77 ~* QQ) the spin functions read = A = T = 0. The spin function are 
evaluatcd in the convention that the traces of the Pauh-niatrices are carried 
out in thrce dimensions. The difference between using three.and n dimensional 
Pauh matrices is a change in the renormahzatipij schemeE3 A similar scheme 
dependence arises in chiral perturbation theoryO 



Currents 

To describe QQ production at NNLL order in e+e" annihilation one needs 
the '^Si vector current with operators of dimension three and five and the '^Pi 
axial-vector current with dimension four, 

= Cl(z.)Op,l(i.) + C2(!/)Op,2(í^) , 

j; = c3HOp,3M, (113) 

where the currents Op^i_3 are given in Eqs. ( p5| ) and (^6|). In this op- 
erators basis there is an additional dimension five vector current, Op, 4 = 

^ í/jp^ (p((T.p) — crí|-) (i(T2) X-p- However it produces a Z?-wave quark- 
antiquark pair and thcrefore does not contribute at NNLL order. The cur- 
rent Op i is dominant and contributes at LL order. The currents Op, 2-3 are 
suppressed and only contribute at NNLL order. Thus ci has to be known at 
NNLL order, i.e. one needs to know the two-loop matching condition and the 
three-loop anomalous dimension. The two-loop matching condition has been 
determined in Ref . |]33| , ^7|| and is displayed íxl Ea . ( p8| ) . The anomalous dimen- 
sion is currently only known at NLL orderEju The coefficients C2 and C3 are 
needed at LL order and have been determined in Ref. [BTl, see Eq. (pq). 



Total Gross Section 
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In fuU QCD the expression for thc total cross section tr^óf (e+e 7*, Z* 
QQ) for quarks at center of mass energy ^/š is 



where apt = 47ra^/(3s). The vector and axial-vector i?-ratios are 

-t í d^xe''^--{0\Tj;{x)f^{0)\0) 



(114) 



R"is) = — Ini 
s 



47r 

R^is) = — Im 
s 



-i / d4xe''--(0 |Tj^(x)j-''(0)| 0) 



(115) 



where q — (y^, 0) and (j'^) is the vector (axial-vector) current that pro- 
duces a quark-antiquark pair. With both 7 and Z exchange the prefactors in 
Eq. (|114|) are 



Q 



rj -2Qf sin^ 9w 
2 sin6'w cos 9w 



s — m\ {s — 



a/ = 



^3 



[s — m^)2 ' 



2 sin 9w cos 6w 



(116) 



Here Q/ is the charge for fermion /, is the third component of weak isospin, 
and Oy/ is the weak mixing angle. 

In vNRQCD at NNLL order the current correlators are replaced by the 
correlators of the non-relativistic currents Op i, so that 

47r r 1 
R%s) ^ —IxYí cl{v)Ai{v,M,v) + 2ci{v)c2{,v)A2{v,M,v) , (117) 
s L J 



Áqr r 

i?''(s) = —Ini cI[v)A:í{v,M,u) 
s L 



(118) 



with 



rOp,i(a;)O;,,(0) 



'^x e*«-^ ( O 



T 



A^^iY^ j(řx e^« " ^ O 



TOpMx)OlJ0) 



0), 
(119) 
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Here q = {y/š — 2M,0). The correlators Ai are functions of thc quark mass 
M, the velocity scaling parameter v, and the velocity 



/i - 2M + 26M^^ 
M 



(120) 



where SM^^ is the LL contribution in the relation between pole mass and the 
threshold mass that is used for the computation. The correlator A2 can be 
related to A\ using the quark equation of motion (see also Ref. |^^) giving 

A2{v,M,v)^v'^Ai{v,M,v). (121) 

The correlators can be computed from the Green function of the Schrodinger 



equation (110) 



Ai{v,M,v) = 6N, J D^pD^p'G{p,p'), 

A3{v,M,iy) = / i?"pDV(p.p')G(p,p'), (122) 

n J 

where Nc = 3 is the number of colors. 

The computation can be divided into two parts. The higher order contri- 
butions coming from the Coulomb potential Vc do not lead to any divergences 
in thc absorptive part of Ai, so they can be determined in 71 = 3 dimensions. 
An extensive technology has been developed to carry out thc corresponding 
computations. For example, in Refs. p8| , p6|j9^ ] various analytic methods were 
used within time-ordered perturbation theory, 

5G{p, P') = - y i?"kiD"k2G,(p, ki)áy(ki, k2)G,(k2, p') + . . . , (123) 

where Gc is the Coulomb Green function that solves the LL Schrodinger equa- 
tion. In a different approach the Coulombic corrections were deter mined by 
solving the corresponding Schrodinger equation exactly. In Ref. [ p8|j99|, 10C] the 
computations were carried out in configuration space, and in Refs. [|57| , ^5[[3^ 
they were done in momentům space based on numerical routines developed in 



Ref. [101|. The corrections to Ai from the l/(m|k|) and the í/rri^ potentials, 
on the other hand, lead to UV divergences, which are associated with the NLL 
anomalous dimension of the leading order current Op i and the evolution of ci. 
Here, the computations need to be carried out in n = 3 — 2e dimensions in the 
MS scheme. In generál, the scheme has to be the same as for the determination 
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of the anomalous dimcnsions and the matching conditions. The rcsults for ^3 
and for the UV-divergent corrections to Ai in n = 3 — 2e dimensions in the 
MS scheme have been determined in Ref. ||3^] . It is also nccessary to treat the 
difFerent orders in SM in the same way as the corrections to the potentials to 
achieve a complete cancellation of the large unphysical corrections associated 
with the pole mass scheme. 



8.2 QQ Production In Fixed Order Perturbation Theory 

In renormaUzation-group-improved perturbation theory the quantity as In v is 
considered as a quantity of order 1. It is summed to all orders in as as indi- 



cated in Eq. (109). After scahng aU coefficients and couphngs of the effective 
theory down from 1/ = 1 to = t; in vNRQCD (or from fip — m to fip = mv 
in pNRQCD) QCD logarithms of v are summed into the coefficients and cou- 
phngs, and matrix elements of the operators are free of QCD logarithms of v. 
This is the approach that I have discussed up to this point in this review and 
that is used in most modern QCD computations involving scale hierarchies. 
Experience has shown that summing potentially large logarithms can lead to 
better controlled QCD predictions. 

It is also possible to consider \nv not as large. In this case the term as Inv 



^ ^ 'yje+e- ^ QQ) 



is counted of order as and the normalized cross section has the generic form 

my 

X |l (LO); as, v (NLO); , a^w, (NNLO)| . 



(124) 



I will call this counting scheme "fixed order perturbation theory" . The indi- 
cated terms are leading order (LO), next-to-leading order (NLG) and next-to- 
next-to-leading order (NNLO) in fixed order perturbation theory and each 
are a subset of the LL order, NLL order and NNLL order corrections in 
renormalization-group-improved perturbation theory. Of course, the count- 
ing asi v ~ 1 is maintained, since the summation of the Coulomb singuarity 
is at the heart of the non-relativistic expansion and is what makes up the 
Schrodinger equation. 

At present, most computations of the total QQ cross section have in fact 
been carried out in fixed order perturbation theory. In Refs. []l2| , p^j5^ ] cutoff 
schemes were employed using NRQCD as the conceptual basis for the com- 
putation. In Ref. [ ^Sf dimensional regularization and the threshold expansion 
were used. Since the cutoff schemes still employed the running QCD coupling 



Heavy Quarkonium Dynamics 



75 



as from fuU QCD, thc computations that used a cutoff were unavoidably de- 
pend cnt on scvcral scales, which are a priori independent. I will report in 
Secs. 9^ and |l0| on the fixed order results for tt production close to threshold 
and non-relativistic bb sum rules. In the fixed order computations only a par- 
tial summation of QCD logarithms is carried out, for example by setting the 
scale of as to a non-relativistic scale. From the point of view of the effective 
theories pNRQCD and vNRQCD these summations are inconsistent. For the 
rest of this section I review the basic technical steps in the computation of 
the total QQ production cross section at NNLO in fixed order perturbation 
theory and discuss the main differences to the renormahzation-group-improved 
calculation. 



Schródinger Equation and Potentials 

At NNLO in fixed order perturbation theory ultrasoft corrections do not enter 
anywhere in the computations. In addition, the coefficients of the operators 
do not run apart from the evolution of ag that is already know from fuU 
QCD. The NNLO non-relativistic QQ dynamics is also described by a common 



time-independent Schródinger equation, which is similar to Eq. (110). Here 



dM = SM^Q + áMj^LQ -I- SM^^j^Q . The potential V is |CDnstructed from the 

The fuU 



NRQCD interactions using Labelle's multipole expansiorJlZl (Sec. p.2[ 
fixed order potential at NNLO for a QQ pair in a ^Si statě reads 

f^(p,p') = Vjip,p') + v/{p,p') + V/Íp,p') + v/ip,p'), (125) 



where (og = as{n)) 

4 t: Cf fls 



V/ÍP,p') 



4 TT 



/?oln(H, 



(Ií)'[''»'"HS)-(^*-+*)'"(S 



0-2 



V/ÍP,P') = 



(126) 

slightly different conventions, 



In some publications (see e.g. Refs. |12| , |98| , p7 
employing off-shell potentials cx (p'^ — p'^) were used. This does, how- 
ever, not affect the finál results. The coefficients of the potentials in the 
fixed order approach do not run except for the evolution of as- To get 



ti'^Cf{'2.Ca-Cf) 
2M|q| 

2^Cf(p^+P^^) 
M2q2 
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Figuře 24: Evolution of renormalization-group-improved potential coefficients in vNRQCD 
(black lineš) and fixed order coefficients (gray Unes) for v <1 \n the tt systém. 



an impression of thc impact of the summation of QCD logarithms, it is in- 
structive to consider the numerical difference between the fixed order coef- 
ficients and the coefficients in renormalization-group-improved perturbation 
theory for scales below M. In Fig. ^ the ratios Vc'^\y) /v'f\l) (sohd hne), 
^k\^)l^k\^) (dashed hne), {V^^\iy) + 2v'^'\v)) / {v['\l) + 2v'^'\l)) (dot- 
ted hne) and VÍ''\v) /vl-'^\l) (dash-dotted hne) are displayed for < 1 for 
M = 175 GeV and as{M) = 0.1074 using four-loop running of the strong cou- 
phng. The resuhs show the evolution of the potential coefficients in vNRQCD 
for the tt systém. For the bh systém the results are similar. The results in 
pNRQCD agree for the \/m? potentials, but disagree for the Coulomb and 
the l/(m|k|) potential. On the horizontál scale of Fig. ^ there is no notable 
numerical difference between the Coulomb coefficient in vNRQCD (vi^'') and 
the corresponding expression in pNRQCD, see Eqs. (86) and (^). In the fixed 
order approach the ratios simply reduce to a1{Mv) / a1{M) for the Vk poten- 
tial (upper gray line) and to {Mv) jas (M) for all the other potentials (lower 
gray line). The curve for as{Mv)las{M') is on top of the sohd hne that dis- 
plays the running oivY^ (y) lv't^ (V) . This shows that the contributions in the 
three-loop anomalous dimension of Ve*"* (y) that are induced by the ultrasoft 
corrections are quite small for the ti systém. 

We find that the evolution in the fixed order and the renormalization- 
group-improved approach roughly agrees for the Coulomb coefficient and for 
the Vá potential. Up to small corrections, both evolve basically with the QCD 
/3-function. For the 14 and the Vy. potentials, however, the evolution of the 
fixed order and the renormalization-group-improved coefficients is profoundly 
different. Whereas the fixed order coefficients grow with as for smaller scales, 
the renormalization-group-improved coefficients become smaller. Expanding 
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for example the vNRQCD coefScient vj^^v) in terms of as{m), we findH 

In the expansion of the corresponding fixed order coefficient the order a^/So 
term agrees, but the third term is missing, because the anomalous dimension 
of the Vk potential is not accounted for properly. Numerically the coefficient of 
the third term is larger than the coefficient of the second term and has a differ- 
ent sign. This explains the diffcrent running for scales below M, and leads to 
a profoundly diffcrent behavior of higher order corrcctions in both approaches. 

Currents and Total Cross Section 

In the fixed order approach the QQ production currents can be defined in 
analogy to Eq. ( |113| ), see e.g. Ref s. [^Št , and the coniputation of the correla- 
tors Ai-3 foUows the lineš of Sec. ^.l| . Since a consistent summation of higher 
order QCD logarithms is not intended, the basic requirement is that the regu- 



larization scheme used in the solution of the Schrodinger equation in Eq. (110) 
is also employed in the matching proceduře. For the currents Op, 2-3 Born 
matching is sufficient. For the dominant current Op i two-loop matching has 
to be carried out. 

For example, in Ref. [p5|] a cutoff scheme was used, where the Schrodinger 



equation in Eq. (IIC) was solved in a three-dimensional momentům space 
sphere with rádius A M. The coefficient ci was then computed by de- 
manding that the total cross section in NRQCD equals the total cross section 
in fuU QCD at order in the limit <C w <C 1. The resulting expression for 
Cl is a complicated function of M and A and also contains the power counting 



breaking effects discussed at the end of Sec. |2.2| . This matching proceduře at 
the level of the total cross section is called "direct matching"E3 and also works 
for ad-hoc regularization schemes that do not have a straightforward definition 
in mathematical terms. For example in Refs. |12| a cutoff A was suggested. 



where the solution of Eq. ( 110 ) was expanded for A —^ 00. Logarithmic diver- 
gent terms were kept and linear, quadratic, etc. divergences dropped by hand. 
This scheme was later adopted by all fixed order computations at NNLO ex- 



cept for Refs. |5Ě,14C]. In this scheme also the power counting breaking effects 



just mentioned are removed by hand. The resulting expression for ci is sim- 
pler than in Ref. and has the same order as contribution as Eq. ( |5^ ) . The 
scale of as in ci obtained in this scheme is independent of the scale of as in 
the correlators. 
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9 Top Quark Pair Production at Threshold 



Top-antitop quark pair production close to the threshold wiU provide an in- 
tegrál part of the top quark physics program at the Linear CoUider, which is 
supposed to be the next major accelerator project after the LHC. It is fair 
to say that top threshold physics was the main motivation for the conceptual 
and technical progress in heavy quarkonium physics in the recent years. The 
theoretical interest in the top-antitop quark threshold arises from the fact that 
the large top quark mass ajid. width {Tt w 1.5 GeV) lead to a suppression of 



non-perturbative effectst22H123 Effectively, the top quark velocity is 



- 2M + 2(5 Aí + iTt 
M 



1/2 



(128) 



SO the hierarchy Mt » MíIí;^,^! ^ Mt|-ř;j,gp » Aqcd is satisfied for any energy 
in the threshold region. This makes top threshold physics an ideál application 
of the non-relativistic effective theories reviewed in the previous sections. In 
particular, perturbative methods are a reliable tool to describe the physics of 
non-relativistic ti pairs and allow for measurements of top quark properties 
directly at the parton level. 

Due to the large top width the total tt production cross section line shape 
is a smooth function of the energy, which rises rapidly at the point where the 
remnant of a toponium IS resonance can be formed. From the energy where 
this increase occurs, the top quark mass can be determined, whereas shape 
and height of the cross section near threshold can be used to determine Tt , the 
coupling strength of top quaxks,to giuons and, if the Higgs boson is not heavy, 
the top Yukawa couplingE2il'E£j For the total top quark width only very few 
other methods to determine it arp Irjinwn From differential quantities, such 

the forward-backward asymmetry or 



as the top momentům distribi| 
certain leptonic distributions 



tirm 



107 



101 



108 



106 



one can obtaitL measurements of Tt , the 



top quark spin and possible anomalous couplingst£2 



A considerable number of experimental studies_wei:e_c. 



IIC 


111 


112 


113 


114 



pd mít in the 

past to assess the feasibility of such measurementslHÍ]'tiil'Ei3lli3ilÍ3 The mea- 
surement of the top quark mass from a threshold lincLshape scan is particularly 
interesting, see Fig. p5| a for an early simulation studylllíJ In contrast to the stan- 
dard top mass determination method, which relies on the reconstruction of the 
invariant mass of jets originating from the decay of a single top quark, the line 
shape measurement has the advantage that only color-singlet tt events have to 
be counted. Therefore, the effects of hnal statě interactions, uncertainties from 
hadronization modeling, etc. are suppressed, and systematic uncertainties in 
the top mass determination are small. Recent simulation studies have shown 
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a) 




.150 CeVMachine+ISR+BeamíLralilung+0.3 íiLi 



b) 




Figuře 25: (a) Sensitivity of the íť line shape, convoluted by the TESLA e+e" beam lumi- 
nosity spectrum, to the top quark mass. The sohd line corresponds to the nominal mass, 
while the dotted lineš to values which difíer from the nominal value b y 4- 900 and ±400 MeV. 
The diagram is taken from an early DESY simulation study in Ref. [11C|, where the "exper- 
imental" points and error bars each correspond to total luminosity of 1 fb~^. Theoretical 
uncertainties are not taken into account. (b) Typical form of the luminosity spectrum at 
the Linear CoUider at nominal beam energy of 350 GeV due to ini tial s tatě radiation, beam- 
strahlung and beam energy spread from Ref. [112|. 



that, for a total luminosity of 100 fb"'^, statistical and systematical expeeka 
tal uncertainties in the top mass determination are well below 50 MeVtl3 
Such a precision would be relevant for example in supersymmetric scenarios, 
where the lightest CP-even Higgs mass caituhave a very strong dependence on 
the top quark mass due to top-stop loopstij 



9.1 Theoretical Achievements 

With the excellent experimental prospect in view it is obvious that a careful 
analysis and assessment of theoretical uncertainties in the prediction of the 
total cross section and various distributions is mandatory, in order to determine 



whether the theoretical ximcision can meet the experimenta 



number of leading ordeillij and next-order computations, 



nit 



i allv a 



LOS 



117 



were carried out. The latter relied basically on QCD-inspired potential modcls 
that used phenomcnological input from T and charmonium data. As such 
they did not represent true first-principles QCD calculations, and thcre was 
no systematic way how the computations could be consistently improved to 
include higher order radiative or relativistic corrections. 

After adoption of the concepts of effective theories in the framework of 
NRQCD (Sec. 0), c onsi ste nt f ixed order NNLO computations appeared in 
Refs. ||,|^j9|||,||,^00|j8|,|ll|,|||. The results were not just some new higher 
order corrections, but led to a number of surprising and important insights. 
The NNLO corrections to the location where the cross section rises and the 
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height of the cross section were found to be much larger than expectcd from 
the results at NLO. The large corrections to thc location of the rise were found 
to be an artifact of the on-sheU pole mass definition and it was reahzed that the 
top pole mass cannot be extracted with an uncertainty smaller than ©(Aqcd) 
from non-relativistic heavy quark-antiquark systems (Sec. 0). Subsequently, 
carefuUy designed threshold mass definitions were suggested to allow for a 
stable extraction of the top quark mass parameter (Sec. 0). The reinaining un- 
certainties in the normalization of the cross section were estimatedĚlI to about 
20% and seemed to jeopardize precise measurements of the top width, the top 
quark coupling to gluons and the Higgs boson. SomerDHLO logarithmic correc- 
tions to the cross section were found to be large toolli3il23 and at present the 
only way to possibly improve the situation in the fixed order approach seems 
to be the determination of N^^LO corrections in order to learn more about the 
structure of the perturbative series. 

An alternativě way is to reorganize the perturbative series by using renor- 
malization-group-improved perturbation theory (Sec. which counts as In v 
of order one and sums all QCD logarithms of w at a certain order into the coef- 
ficients of the effective theory. Such a treatment has become standard in many 
other areas of QCD where physics is governed by widely separated scales. The 
theoretical framework that allows to carry out such a summation consistently 
is more sophisticated than NRQCD (Sec. |2|). Currently, there are two alter- 
nativě formulations in the literatuře, called pNRQCD (Sec. ^ and vNRQCD 
(Sec. H). The two effective theories are not equivalent (Sec. ||). Recently, 
a renormalization-group-improved calculation in vNRQCD of the total cross 
section at NLL order and including most NNLL order corrections has been 
provided in Refs. It was found that the normalization uncertainties are 

an order of magnitude smaller and the remaining relative uncertainty in the 
normalization of the cross section was estimated at the level of 3%. 

In the foUowing sections I discuss how the top quark width and non- 
perturbative effects are implemented into computations of the non-relativistic 
tt dynamics, and I give a more detailed discussion of results for the total cross 
section in fixed order perturbation theory and the recent results obtained in 
vNRQCD. 



9.2 Top Width and N on- Perturbative Effects in the Total Cross Section 

For computations of the total cross section the top width can be understood as 
a modification of the coefHcients of the effective theory caused by electroweak 
corrections, since one is not interested in the dynamics of the top quark decay. 
(In generál, the top decay would be associated with external currents similarly 
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to the current that produces the tt pair.) Because the particles involvcd in 
these corrections are lighter than the top quark, they can lead to non-zero 
imaginary parts in the coefhcients of the effective theory. This is a known 
concept in quantum mechanics of inelastic processes. In the context of top 
quark threshold production the effects of the top quark decay are the most 
important part of a whole array of electroweak corrections. 

The dominant effect of the top quark width is obtained by including the 
electroweak top quark selfenergy in the matching conditions for the coefíicients 
of the bihnear top operators. Since the matching is determined for on-shell top 
quarks, the dominant effect of the top decay is that of an imaginary mass term 
proportional to the top quark width being added to the effective Lagrangian. 
For example, in vNRQCD the operators 

<5£ = 53 V^t ^r, Vp + E 4 xp . (129) 
p p 

have to be added to the Lagrangian, Tt being the on-shell top decay width. In 
the Standard Model the dominant decay channel is í bW^ and gives |a_width 
of Tt = 1.43 GeV for íTip^j^, = 175 GcV including one-loop electroweaklllll (see 
also Rcf. |]l22[| ) and two-loop qcdEÍ effects. In the effective theories, Tt is 
treated as a fixed parameters that can be determined from experiment. Since 
the typical energy of the top quark is E ^ M|t;^jjp 4 GeV, it is natural in 
the Standard Model to count Tt of order Mw^. Thus, the propagátor of a top 
quark with momentům (p'^ , p) is 

; • (130) 



pO^£^+SM + ^Tt+ iS 

In extensions of the Standard Model, where the total top width is different, the 
u-counting might hpzE to be modified. The use of this propagátor is equivalent 



to the replacementlii!3 E ^ E + iTt in tt Green functions obtained from the 



Schrodinger equation for stable quarks in Eq. (|llO[ ), i.e. it is sufRcient to replace 



v by the effective velocity v^g defined in Eq. (128) in the Schrodinger equation 



of Eq. (IIC). For the total tt cross scction the operators in Eq. (129)_have 
been shn wn to be sufRcient to account for all electroweak effects at LOE23 and 
NLOli£j in the fixed order approach. At NNLO the treatment of electroweak 
effects is more involved, since also first order electroweak corrections to top 
quark production and decay and background processes consistent with the top 
decay finál statě need to be taken into account. The complete set of elec- 
troweak effects at NNLO order is currently unknown—However, some one-loop 



electroweak corrections to the tt production current 



121 



and some non-resonant 
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backgroimd contributionslirJ have been computed. Also, the structure of some 
higher order efFects associated with the electroweak top quark selfenergy is un- 
derstood. The next order corrections from the selfenergy leads to a top wave 
function correction in the efFective theory, which is up to a difFerent sign equal 
to the corresponding vertex correction to the gluon-quark coupling contained 
in the Coulomb potential. In computations of energy levels both contributions 
cancel due to gauge-invariance. The remaining subdominant width effect arises 
from the non-relativistic expansion of the top quark spinors of fiill QCD. In 
vNRQCD this leads to the foUowing additional operators 

2 ■ 2 

sc = -J2^l '-r, ^v^p ~Y.4 2^xp . (131) 
p p 

Since the effective Lagrangian describes the non-relativistic tt dynamics in the 
center-of-mass frame, these operators can be interpreted as a time-dilatation 

2 

effect ~ rÉ(l — (j^p-) + . . .) that arises from the Fermi motion of the in- 
dividual top quarks in the tt systém at rest. The cancellation bctween the 
electroweak wave function corrections and the electroweak vertex corrections 
to the Coulomb potential was verified at the onc-loop level by explicit calcula- 
tion in Ref. [127| (see also Rcf. ]12^ ]). For muonic atoms it was noted before 
in Rcf. |12£] and later again in Ref. [130|. An analogous discussion for the Bc 
lifctime can be found in Ref. [131|. 

The large top quark width plays a crucial role in suppressing the size of 
non-perturbative hadronic contributions governed by the scale Aqcd- Even 
arbitrarily close to the threshold point, the tt systém has an effective energy 
set by the perturbativei^ale Tt, see Eq. (128). For stable quarks, an oper- 
átor product expansionElIj in powers of Aqcd /E can be used to incorporate 
non-perturbative contributions. The first non-pcrturbative correction arises 
from the radiation and absorption of a gluon with an energy of order Aqcd . 
In the framework of pNRQCD and in configuration space representation the 
correction can be understood from the inscrtion of two r.E interactions, where 
the gluon field carries no energy and momentům as shown in Fig. Due to 



Figuře 26: Leading non-perturbative correction in the operátor product expansion due to 
the vacuum average of two chromoelectric fields. 



Lorentz and gauge invariance the one-particle-irreducible two-point function 
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of thc two chromoclcctric ficlds is related to the vacuum average (or vacuum 
condensate) of two gluon ficld strength tensors, 

(o|ií;,^ií;f |o> = -lá.,á^^(o|G^,G^''^|o) . 

The resulting leading correction to the color singlet Green funrtion of the 
Schródinger equation in configuration space representation readsO 

(5G,(r,r')--^(0|a,G^,G^nO) J d^^ Jd^y^.yGs{r,^)Go{^,y)G,{y,r'), 

(132) 

where Go stands for the Green function of the Schródinger equation with the 
octet Coulomb potential. For tt production the relative size of the gluon con- 
densate correction is of order [Aqqb / {E + iT t)]'^ x v'^^ ~ 10~^. The first factor 
arises for dimensional reasons and the second term from thc p.A/M interac- 
tion that comes from the multipole expansion. The radiation of the gluon off 
the potential (Fig. |^) gives the same estimate. An explicit expression for the 



gluon condensate corrections to the correlator Ai was derived in Ref. [ |l32 | 

9 i Clcf 



5A,{v,s,M,y) = ^^^J^ (0|a.G^.G^^|0) 



X V r(n + 4)r2(n-A) 

(n + 2+|)r(n + l)r2(n + 5-A) 



where A = —CfOís/{'2,v^q) and F is the gamma function. Thc size of this correc- 
tion for the (conservative) literatuře value (0|as G^^G'^'' |0) = 0.05 ±0.03 GeV'' 
is consistent with the previous parametric estimate and can be safely neglected 
in view of the current perturbative uncertainties. It is a quite interesting ques- 
tion, whether this estimate of the size of non-perturbative corrections is indeed 
correct. 

9.3 Gross Section at NNLO in the Fixed Order Approach 

In the fixed order approach NNLO computations of the tt cross section clofifi 
to threshold in eie~ annihilaiion have been cajried out by Hoang et al.,LZfĚ3 
Melnikav et al.j£3 Yakavlevj£3 Beneke et al.,Ěj Nagano et al.,E2íj and Penin 
et al£Í3 Penin et al.lli3 have also provided a NNLO computation of the tt 
cross section in 77 collisions. Since the fixed order computations basically use 
NRQCD as the conceptual framework (Sec. 0), which does not fix a preferred 
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convention for the regularization of UV divergences, the results obtained by 
the various groups are, partly, quite difFerent. The differences of the resuhs 
in the various conventions and regularization schemes can be interpreted as 
differences in the treatment of corrections from beyond NNLO and serve as 
an important tool to estimate the inherent uncertainties in the fixed order 
approach. 

In Ref. \ ^Ť\ Hoang et aL solved the momentům space Schrodinger equation 
in Eq. (IIC) exactly for the Coulomb potential based on the numerical routines 
of Ref. [101|; the other l/m suppressed contributions whcrc calculatcd with 
time-ordered perturbation theory using a momentům cutoff keeping only terms 
that depend logarithmically on the cutoff and discarding terms that depend 
hnearly, quadraticaUy, etc. on the cutoff. In Ref. Hoang et al. solved 
the complete NNLO Schrodiner equation exactly in a momentům sphere with 
rádius A ~ M. In Ref. |9|] Melnikov et al. solved the NNLO Schrodinger 
equation exactly in configuration space representation for the Green function 
G{0, To), where the distance tq served as the cutoff parameter. They also kept 
only terms that a logarithmically dependent on tq and discarded power-like 
divergences. The method was also employed by Yakovlev in Ref. and 
Nagano et al. in Ref. [100|. A similar regularization method was employed by 
Penin et al in Refs. [118], but the entire Schrodinger equation was solved with 
time-ordered perturbation theory supplemented by partial resummations ta 
avoid multiple energy poles. AU methods used the direct matching proceduralj 
to determine the coefhcient ci of the dominant tt production current. The 
results obtained in these cutoff schemes had a residual dependence on three 
different scales, the renormalization scale ^ of as in the potentials, the cutoff 
and the renormalization scale of as in the coefficient ci . A different stratégy 
was used by Beneke et al. in Refs. [|5^. They used dimensional regularization 
identifying the contributions from the various momentům regions in Eq. (^l|), 
i.e. potentials, coefRcients, etc, using the threshold expansion and solving the 
Schrodinger equation in time-ordered perturbation theory. In this approach 
the potentials and cocfficients were determined without an explicit matching 
calculation. The result had a residual dependence on two different scales, the 
renormalization scale of tts in the potentials, and the factorization scale. In 
addition, a summation of some NLL order logarithms was included in ci . The 
prescription used for this summation, however, is not consistent with the known 
NLL order anomalous dimensions in pNRQCD or vNRQCD. In addition, ci 
was not included as a globál factor of the correlator Ai but expanded out. 

In Ref. |8^] the results for the normalized total photon-induced cross sec- 
tion Q'^R'' = cr(e+e- ^ 7 ^ tt)/a^t of all groups were compiled and com- 
pared numerically in detail using an equivalent set of parameters to analýze the 
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Figuře 27: The total normalized photon-induced tt cross section R at the Linear Collider 
versus the cm. energy in the threshold regime at LO (dotted curves), NLO (dashed) and 
NNLO (solid) in the pole mass scheme for mP°'° = 175.05 GeV, o^(Mz) = 0.119 , Tt = 
1.43 GeV and = 15, 30, 60 GeV. The plots were given in Ref. | |8l[ from results provided 
by Hoang-Teubner (HT), Melnikov-Yelkhovsky-Yakovlev-Nagano-Ota-Sumino (MYYNOS), 
Penin-Pivovarov (PP) and Beneke-Signer-Smirnov (BSS). 



scheme-dependence of the fixed order approach. Since the axial-vector current 
contributions are only a smaU correction at the few percent level it is justi- 
fied to consider only the photon-induced cross section. Figuře ^ shows the 
results obtained in Ref. in the pole mass scheme for m^^^^ — 175.05 GeV, 
asiMz) = 0.119, Tt = 1.43 GeV and ^l = 15, 30 and 60 GeV. The other scales 
where fixed at A/poio- The effects of the luminosity spectrum of the e^e~ 
beams, which leads to a smearing of the line shape, were not taken into ac- 
count. Since the method used by Melnikov et al., Yakovlev et al. and Nagano 
et al. were equivalent, a single figuře is displayed for their results. The result 
of Hoang and Teubner was based on the momentům sphere method of Ref. [85| 
and has the smallest variations. In particular, the LO result and the NLO and 
NNLO corrections are smaller than for the other groups. The results in Fig. 
demonstrate the instability of the "peak position" in the pole mass scheme 
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and the necessity of using thrcshold masscs. It was concluded in Ref. |^] that 
the perturbative uncertainty in an extraction of the pole mass parameter from 
the peak position (i.e. when the beam smearing effects are neglected) is around 
300 MeV. The results in Fig. ^ also show a large uncertainty in the normal- 
ization of the cross section. This uncertainty is particularly puzzhng, because 
there is no obvious physical reason for its existence. At present it does not 
seem to be related to renormalon-type higher order corrections, although it 
has also been sp^pulated that the large size of the corrections co uld h ave some 
infrared originiHj On the other hand, it has been shown in Ref. |120|1 that the 



dominant N^LO logarithmic corrections to the cross section of (relative) order 
Qíg In^ z; are at the level of 10% and not small. It was estimated in Ref. |81 



that the present theoretical relative uncertainty in the normalization of the 
cross section in the fixed order approach is around 20%. 

In Ref. [RlJ it was also studied how well the peak position can be stabilized, 



when threshold mass schemes are employed (Sec. 7^2). Figuře |2^ shows the 
results for the photon induced total cross section obtained in Ref. |8^] where 
the respective values of the threshold masses were obtained from the top MS 
mass value m(m) = 165 GeV as a reference point. Numerically, the series (in 
an expansion in powers of as for /i = m) that relate the pole punass and the 
threshold masses used in Fig. G§ to the MS mass are as followsM 



"^pole 



; 165 + 7.64 + 1.64 + 0.52 + 0.25(*' ] GeV 
Mis = [ 165 + 7.21 + 1.24 + 0.23 + 0.05^*^ ] GeV 
Mpg(20 GeV) = [ 165 + 6.72 + 1.21 + 0.29 + 0.08(*) ] GeV 
mki„(15 GeV) = [ 165 + 6.68 + 1.15 + 0.27^*^ ] GeV . (134) 



The numbers indicated by the superscript (*) are estimated in the large-/3o 
approximation and as{Mz) = 0.119 (which corresponds to q;s(165 GeV) — 
0.1091 using four-loop running) has been used. The order terms for the 
threshold masses are about a factor two smaller than for the pole mass. The 
results shown in Fig. ^ demonstrate that the peak position is significantly 
stabilized when threshold masses are employed. It was concluded in Ref. |81[ 
that the perturbative uncertainty in the determination of threshold masses 
from the peak position (i.e. when beam smearing effects are neglected), is 
between 50 and 80 MeV. It was also pointed out that the MS mass can be 
determined with a comparable perturbative precision, only if as{Mz) is known 
with an uncertainty of around 0.001. This restriction arises from the relatively 
large order as correction in the relation between threshold masses and MS 
mass. For example, for a given measurement of the IS mass, leťs say M^g = 
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•[qHGeV) ^(GeV) 

Figuře 28: The total normalized photon-induced tt cross section R at the Linear Col- 
lider versus the cm. energy in the threshold regime at LO (dotted curves), NLO (dashed) 
and NNLO (sohd) for a/{Mz) = 0.119 , Tt = 1.43 GeV and fi = 15, 30, 60 GeV. 
Hoang— Teubner used the IS mass M-^^ = 173.68 GeV, Melnikov— Yelkhovsky the kinetic 
^kin(15 GeV) = 173.10 GeV, and Beneke-Signer— Smirnov arid_ Yakovlev the PS mass 
Mps(20 GeV) = 173.30 GeV. The plots were given in Ref. ||lj from results provided 
by Hoang— Tcubncr (HT), Melnikov- Yelkhovsky (MY), Beneke— Signer-Smirnov (BSS) and 

Yakovlev. 



175 GeV±(5Mig, and UsiMz) = 0.118 ±x 0.001 the result for m(m) readřl 



m{m) 



175- 7.58 - 0.96 - 0.23 ± SM^^s ± a; 0.07 



GeV 



(135) 



where the first four numbers represent the perturbative series. 
The results in Fig. 



28 



also show that the threshold masses do not affect the 
normalization of the cross section; the uncertainty of around 20% remains in 
threshold mass schemes. This can be understood from the fact that the transfer 
from the pole mass to threshold masses corresponds predominantly to a shift 
in the threshold point. The large normalization uncertainties could have severe 
consequences for experimental measurements. For the measurement of the top 
quark mass it is conceivable that, after the luminosity spectrum of a speciíic 
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collider design is folded into thc calculations for the cross section, the resulting 
"peak-less" line-shape (see Fig. p5^) could lead to a transfer of normalization 
uncertainties into the top mass measurement. In generál, such beam effects 
are smaher for /í^/í^ colhders, because beamstrahlung is suppressed by the 



muon mass, see e.g. Ref. |134]. On the other hand, the large normaUzation 
uncertainties would make precise measurements of a^, Tt or the top quark 
Yukawa couphng to the Higgs, yt impossible. In the fixed order approach, 
effects of a Standard Model Higgs can be implemented by a Higgs ^iotential 
Vj.jjj = —yf/[2{q^ + mf^)] and a correction to the coefhcient cill23'E2£l For a 
light Higgs they are typically at the level of several percent. Figuře ^ shows 
the variation of the theoretical total cross section a{e^e^ Z ^ tt) for 

changes in the input parameters as, Tt and yt- 

The question how well the top mass and the strong coupling can be detei: 



e collisions with the TESLA-LC design 



135 



mined from a threshold scan ir 

was addressed by Peralta et al.liiíl in a simulation study based on the theoreti- 
cal NNLO calculations for the cross section from Hoang et al. in Ref. ||8^]. The 
theoretical computations from Ref. [|5| have the smallest theoretical variations 
of all NNLO computations, so the estimates for the theoretical uncertainties 



obtained in the experimental study of Ref. |112] are likely to be larger, if 



the NNLO computation of the other groups are used in a similar analysis. 



However, the results of Ref. |113] illustrate very well the dependence of the 
measurements on theoretical uncertainties. In the simulation Peralta et al. 
included the TESLA luminosity sprectrum and realistic assumptions on efR- 
ciencies, background and other systematics. In the study nine scan points were 
placed onto the tt line shape similar to Fig. p5| a and a tenth scan point was 
placed well below the threshold to evaluate the background. At each point an 
integrated luminosity of 10 fb~^ was spent. For a two-parameter fit of top mass 
and as{Mz) Peralta el al. found the correlations shown in Fig. |29|a. The error 
ellipses correspond to Ax^ = 1. Because the corrections to the '^S'i toponium 
ground statě mass are zero in the IS scheme, the correlation is smallest for the 
IS mass leading to the smallest experimental uncertainties. The experimental 
uncertainty is 40 MeV (50 MeV) for M^g (A/ps(20 GeV)) and 120 MeV for 
the pole mass. The correlation is larger in the pole mass scheme. In order 
to test the impact of the theoretical uncertainties in the normalization of the 
cross section the fits were repeated for "experimental" data generated from 
the LO, NLO and NNLO theoretical predictions and for three different values 
of the renormalization scale, /z = 15, 25, 60 GeV. The difference in the loca- 
tion of the ellipses are a measure for the theoretical uncertainties; see Fig. p9| b 
for the result in the IS scheme (Ax^ = 1). For M-^^ and Aípg(20 GeV) the 
theoretical uncertainties were found to be at the level of 10-40 MeV, whereas 
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Statistical error-- 




b) 



Systematic errors. 1 5 Míiss 




m, (GeV) rr, (GeV) 

Figuře 29: (a) Results of a two-parameter (M-Os(Mz)) fit (Ax^ = 1) for the pole, IS 
and PS (Aíp3 = 20 GeV) masses using the same NNLO theoretical cross section {fi = 25 GeV) 
for "experiment" and theory. The eUipses show the experimental errors. (b) Result for a 
two-parameter M-^g-as(Mz) X^ St {^x^ = 1) using difFerent theoretical cross sections for 
the "experiment" while keeping the theory fixed to NNLO with /i = 25 GeV. The shifts 
in the locations of the ellipses show the theoretical errors. The a nalv sis was carried out in 
Ref. [1131 based on the computations in Ref. 



the error in as{Mz) was 0.013, which is several times larger then the current 
world average. For the pole mass the theoretical uncertainty was found to be 
about 100 MeV, and the error in as{Mz) was 0.012. This result shows that 
in a two-parameter fit the large normalization uncertainties of the total cross 
section are fed almost entirely into as and essentially leave the measurement 
of threshold masfajmrameters unafFected. 

Peralta et al.ti3 also carried out a constraint fit taking ag{Mz) — 0.120± 
0.001 as an input. Here it was found that the experimental uncertainties for 
Mj^g, Mpg(20 GeV) and the pole mass are slightly smaller than in the two- 
parameter fit, but that the theoretical uncertainties for the masses are about 
3-4 time larger. This behavior can bc understood from the corrclation between 
masses and as shown in Fig. p9| a and from the fact that the largo normalization 
uncertainties of the total cross section are partly absorbed into the mass for a 
constraint fit. 



In a more recent simualtion study by Martineztlll the size of the experi- 
mental uncertainties for simultaneous measurements of the \S mass, as{Mz)^ 
Tt and yt in a top threshold run at an e^e~ Linear Collider was reexamined 
taking into account "data" on total cross section, top three-momentum dis- 
tribution and forward-backward asymmetry for a 9+1 point threshold scan 
using a total integrated luminosity of 300 fb~^. The intention was to deter- 
mine which experimental uncertainties could be eventually achieved, and not 
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to estimate theoretical uncertainties. Fixing yt to the Standard Model (SM) 
value Martinez et al. obtained AMis = 18 MeV, Aas{Mz) = 0.0012 and 
ATt = 30 MeV from a three-parameter fit. In a one parameter fit, fixing all 
other parameters, they obtained 5yt/yt between 12 and 24%. In a fit where Tt 
is fixed to the SM value and as is constrained to as{Mz) = 0.120 ±0.001 they 
obtained (5M^g = 25 MeV and Syt/yt = (+0.31,-0.49). In a fit where only 
as is constrained to as{Mz) = 0.120 ± 0.001 they obtained 5M^g = 30 MeV, 
STt = 33 MeV and 6yt/yt = (+0.33, -0.57) for a Higgs mass of = 120 GeV. 
For simulation studies for fj,^^~ colliders see for example Ref. [134|. 

At present it appears questionable whether the experimental precision for 
measurements of as, Fj and yt can be matched by the precision of the theo- 
retical computations in the fixed order approach. 



9.4 Cross Section at NNLL Order in vNRQCD 

At present there exists only one computation of the total tt production cross 
section closej-tp threshold in the renormalization-group-improved approach by 
Hoang et al.É3 in the framework of vNRQCD. Since ci does not run and the 
Coulomb potential evolves trivially with at LL order, the LL cross section 
in the renormalization-group-improved approach and the LO fixed order cross 
section are equivalent. The NLL cross section in vNRQCD and the NLO fixed 
order cross section differ because ci runs non-trivially at NLL order. Thus 
the NLL order cross section and the fixed order NLO cross section differ by a 
normalization factor. At NNLL order all relevant coefficients for potentials and 
currents run non-trivially, and the difference to the fixed order NNLO results 
can be significant with respectJ.o the normalization as well as to the shape of 
the cross section. Hoang et alO determined the cross section at LL, NLL and 
NNLL order. The NNLL order result is incomplete, since it does not include 
the three-loop evolution of ci, which is presently unknown. 

In Refs. \RT\i the calculations were carried out in _D = 4 — 2e dimensions in 



the MS scheme. The Schródinger equation in Eq. (IIC) has been solved exactly 
with numerical methods for all Coulombic contributions, and the corrections 
from the kinetic energy term p^/(4M'^) and the 1/m-suppressed potentials 
were determined in time-independent perturbation theory. In Fig. ^o| a the 
normalized photon-induced cross section obtained in Ref. ||3^ ] is displayed at LL 
(dotted lineš), NLL (dashed lineš) and NNLL (solid lineš) order in the IS mass 
scheme for M-^g = 175 GeV and four different choices of the velocity scaling 
parameter, 1/ = 0.1, 0.125, 0.2, 0.4. The peak position is stable due to the use of 
a threshold mass scheme. In contrast to the fixed order computations discussed 
in the previous subsection an excellent convergence of the normalization and 
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Figuře 30: (a) Normalized photon-induced cross section Q^R" at LL (dotted lineš), NLL 
(dashed Unes) and NNLL (solid lineš) order for M^g = 175 GeV, as{Mz) = 0.118, Tt = 
1.43 GeV and u = 0.1,0.125,0.2,0.4 at each order. (b) Total NNLL order cross section 
in units of pb for thc same parameters. The dotted, dashed, dot-dashed, and solid curves 
correspond to u = 0.1,0.15,0.2,0.4. The figures were n hta ined with renormalization-group- 
improved perturbation theory in Ref. [p7| using vNRQCD. 



a decreasing sensitivity to changes in v at higher orders is observed. It was 
shown that this behavior is a consequence of the evolution of the vNRQCD 
coefficients, which sum QCD logarithms oi v ií v is chosen of order v ag. In 
particular, the non-trivial evolution of the potential coefficients (Fig. ^) leads 
to a partial canceUation of NNLL order corrections because some coefficients 
develop a relative sign diffcrencc when v is lowered. Such a behavior does not 
exist at NNLO in the fixed order approach, where all coefficients scale trivially 
with ag and grow in the same relative direction when the scale /i is lowered. 
In Fig. pO| b thc total NNLL order ti cross section including 7 and Z exchange 
in units of pb determincd in Refs. is displayed for i/ = 0.1, 0.15, 0.2, 0.4. 

From the size of the NNLL order corrections, the í^-variation of the the cross 
section and an examination of some known higher order corrections, Hoang et 
al. estimated the remaining perturbative uncertainty of the normalization of 
the cross section as 3%, which is a reduction by almost an order of magnitude 
with respect to estimates in the fixed order approach. If the uncertainty of the 
renormalization-group-improved approach is established, the prospects for the 
size of theoretical uncertainties in measurements of top mass, as, Ft and the 
Yukawa coupling yt from a top threshold scan are significantly better than in 
the fixed order approach. In particular, for fits where the top couplings are 
fixed, the danger that the perturbative normalization uncertainties are fed into 
the top mass measurement is avoided. On the other hand, measurements of as, 
Tt and yt with acceptable theoretical uncertainties appear only to be possible 
in the renormalization-group-improved approach. In Figs. |^ the variation of 
the NNLL cross section cr^o^ obtained in Ref. is displayed as a function of 
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Figuře 31: Variation of the NNLL cross section in vNRQCD | |37[ without boam cíTects for a) 
the value of the strong coupling, b) the top quark width, and cj the inclusion of a Standard 
Model (SM) Higgs boson. Changes relative to the centrál value (solid lineš) are shown by 
dashed red lineš. In c) there are two solid lineš, the lower black line is the decoupling limit 
for the Higgs boson, and the upper blue line is for a SM Higgs with mass ttih = 115 GeV. 
The parameters Mis = 175 GeV, Tt = 1.43 GeV, as{Mz) = 0.118, j/j = O and ^ = 0.15 
have been used unless stated otherwise. 



for difFerent choices of the input parameters as{Mz), and yt- In all fig- 
ures the parameters M^^ = 175 GeV, Tt = 1.43 GeV, a,{Mz) = 0.118, ?/* = O 
and v = 0.15 are chosen, unless stated otherwise. The upper left panel shows 
the cross section for as{Mz) = 0.116 (lower dashed red line), 0.118 (solid black 
line), and 0.120 (upper dashed red line). At the peak one finds a ±2.7% varia- 
tion when varying as{Mz) by ±0.002. The upper right panel shows the cross 
section for Tt = 1.43 GeV (solid black line). The dashed red lineš correspond 
to variations of by ±10%, where for a smaller width the peak becomes more 
pronounced. At the peak one finds a (—2.3%, +2.6%) variation of the cross sec- 
tion when varying Tt by ±5%. Finally, the lower panel shows the cross section 
for zero Yukawa coupling (solid black line), the Standard Model value for the 
Yukawa coupling (solid blue line), and a ±20% variation of the coupling with 
respect to the SM value (upper/lower dashed red lineš). The Higgs mass was 
chosen to be mn = 115 GeV. At the peak one finds a (+3.3%, —2.6%) variation 
when varying yf^ by ±20%. For ttih = 130 and 150 GeV the corresponding 
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variations are (+2.9%, —2.3%) and (+2.5%, —2.0%) respectively. Thus, given 
that the remaining theoretical uncertainty in the normalization of the cross 
section is 3%, measurements with theoretical uncertainties das{Mz) ~ 0.002, 
STf /Tt 5% and Syt/yt ^ 20% appear feasible. These uncertainties match 
quite weh the estiraated experimental uncertainties from a recent simulation 
study by MartineziHj (see end of Sec. 9.3). 



10 Bottom Mass from T Sum rules 

The determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 
ments is one of the main goals of current B physics experiments. A sufficiently 
accurate determination of the size of the CKM matrix elements and their rel- 
ative phases will lead to a better understanding of the origin of CP violation, 
the structure of the weak interaction, and, possibly, to the estabhshment of 
physics beyond the Standard Model. For the extraction of the CKM matrix 
elements from inclusive B decay rates, such as of Vet from the semileptonic 
decays into D mesons or of Vub from semileptonic decays into light hadrons, 
an accurate and precise knowledge of the bottom quark mass is desirable due 
to the strong dependence of the total decay rate on the bottom quark mass 



parameter (Sec. 7.2) 



Non-relativistic sum rules for the masses and electronic decay widths of 
T mesons, bottom-antibottom quark bound states with '^'^^^Lj = '^Si, are an 
ideál tool to determine the bottom quark mass: using causality and globál 
duality arguments one can relate integrals over the total cross section for the 
production of hadrons containing a bottom-antibottom quark pair in e~^e~ 
coUisions to derivatives of the vacuum polarization function of bottom quark 
currents at zero-momentum transfer. The n-th derivative is referred to as 
the n-th moment of the total cross section. For a particular range of n the 
moments are dominated by the experimental data on the T mesons and, at 
the same time, can be calculated reliably using perturbative QCD in the non- 
relativistic expansion. Because the moments have, for dimensional reasons, a 
strong dependence on the bottom quark mass, P„ ~ these sum rules 

can be used to determine the bottom quark mass to high precision. Within the 
last few years there have been several analyses at NLO and NNLO in the non- 



relativistic expansiop 
mass was extracted 



n in the fivpd order approach. TnitiaHy the bottom pole 
l^'E£3aEa whereas later analysesollB^EfjQ determined 



E2a whereas later analyses 
threshold mass parameters and the MS mass m(m). For an application of 
non-relativistic sum rules to the V' family for an extraction of the charm quark 
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mass see Ref. [142]. 



10.1 Basic Theoretical Issues and Experimental Data 

The sum rules for the T mesons start from thc corrclator of two electromagnetic 
vector currents of bottom quarks at momentům transfer q, 

n{q'^)^-i J dxe"^-'' {0\Tf^{x)f''{0)\0) . (136) 

The n-th moment P„ of the vacuum polarization function is defined as 

(137) 



dq"^ I q^ 



q^=0 



where Qfc = — 1/3 is the electric charge of the bottom quark. Due to causaHty 
the n-th moment P„ can also be written as a dispersion integrál 

oo 

Pn = Ql J ^^"(S), (138) 

where 

QlR^is) = -i^^^^^-^^^bb+X) ^ 4^j^n(.) (139) 

CTpt S 

is the total normahzed photon-induced cross section of bottom quark- 
antiquark production in e+e^ annihilation. The lower Hmit of the integration 



in Eq. (138) is set by the mass of the lowest-lying 66 resonance. Assuming 
globál duality, the moments P„ can be either calculated from experimental 
data on R or thcorctically using perturbative QCD. 

The experimental moments are determined by using the data on the 
T meson masses, Mj^g, and e"*'e~ partial widths, Fj^s, for k = 1,...,6, see 



Ref. [143|. The formula that is used for the experimental moments has the 



typical form 

6 ^ 

Pn-" = ^ E -TTĚn + I 4^1 -cont(s) , (140) 
"cm ^1 ^^kS J * 

and is based on the narrow width approximation for the known T resonances; 
ácTíí is the electromagnetic coupling at the scale 10 GeV. Because the differ- 
ence in the electromagnetic coupling for the different T masses is negligible, 
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nS 


Mns/[GeV] 


rns/[keV] 


IS 


9.460 


1.32 ±0.04 ±0.03 


2S 


10.023 


0.52 ±0.03 ±0.01 


3S 


10.355 


0.48 ±0.03 ±0.03 


4S 


10.58 


0.25 ±0.03 ±0.01 


5S 


10.86 


0.31 ±0.05 ±0.07 


6S 


11.02 


0.13 ±0.03 ±0.03 



Table 1: The experimental numbers for the T masses and electronic decay widths used for 
the calculation of the experimental moments P^". For the widths, the first error is statistical 
and t he senori d one systematic. AU errors are estimated from the numbers presented in 
Rcfs. [143 144 1. The electromagnetic coupling at 10 GeV and the BB threshold point are 
= a;:J(WGeV) = 131.8(1 ± 0.005) and (v^)bb = 2 x 5.279 GeV. The errors in the T 
masses and the BB threshold (\/š)bb small and can be neglected. 



one can chose 10 GeV as the scale of the electromagnetic coupling for all res- 
onances. Since the continuum cross section above the BB threshold is quite 
poorly known experimentally a nuniber of different methods to estimate it has 
been used in the hterature. Usually it is approximated by some theoretically 
motivated model, based on the perturbative predictions for the continuum re- 
gion. For very low values of n < 4 the uncertainties in the knowledge of the 
continuum cross section are a significant contribution in the uncertainties of 
the bottom mass extractions. For n > 4 the continuum contribution is sup- 
pressed sufficiently so that only a rough estimate of the continuum region is 
needed. Here, it is sufficient to assume that Tc is a constant compatible with the 
perturbative predictions. In Tab. |^ a compilation is given of the experimental 
numbers for the T resonances needed in Eq. (140). 



A reliable computation of the theoretical moments based on pertur- 
bative QCD is only pos sible if the effective energy range contributing_tD the 
integration in Eq. (138) is sufRciently larger than Aqcd ~ 300 MeVEž3 For 
large values of n the theoretical moments are dominated by the non-relativistic 
kinematic regime close to the T resonances, and it can be shown that the size 
of the energy range is proportional M /n. This implies that n should be chosen 
sufRciently smaller than 15-20. A practical upper bound for n is 10, which 
has been used in Refs. pŠIPlliČI. ,|73| . L arger values for n up to 20 have 
been used in the analyses of Refs. [136 137, 96| , p7| . One argument in favor of 
using also values of n larger than 10 is that the leadiag xtcder gluon conden- 



sate contribution of the large-n theoretical momentslža-tžS (see Fig. and 
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Eq. (^)), 

Pn""^^ = -Ť|^e-°-4''""=^(0|a,G,.G^nO)Pr- + ■•■ , (141) 

where the P^'^' is the n's moment obtained from perturbation theory, is less 
then a percent even for n < 20. On the other hand, the estimate of non- 
perturbative efFects based on a local expansion in gluonic and hght quark 
condensates is only rehable if M/n ^ Aqcd, see Ref. In addition, the size 
of the non-perturbative corrections coming from the dimension-6 (and higher) 
condensates as well as all perturbative Wilson coefficients for the condensates 
are not well known. It has also been argued recently that tljp condensates 
might actually not be the dominant non-perturbative effectsE 



147 



In order to suppress uncertainties from non-perturbative corrections, n 
should be chosen as small as possible. If n is chosen very small, typically 
n < 4, the dynamics of the bb pair encoded in the moments is dominated by 
relativistic energies and momenta of order M , and the usual perturbative ex- 
pansion in the number of loops can be employed to determine the theoretical 
moments. At present the normalized photon-induced heavy quark pair pro- 
ductioH-jCross section including the fuU mass dependence is known at two-loop 



orderEda Very few analyses in terms of Iow-ti moments exist in the literatuře. 



A recent analysis in terms of Iow-ti moments and extracting the bottom MS 



mass as m(rn) = 4.21 ± 0.05 GeV was carried out in Ref. [149|. Since this re 



view concentrates on non-relativistic dynamics, no discussion on this method 



can given here. It should be noted that the result from Ref. [149| is consistent 
with the large-n sum rule analyses that used threshold masses (see Tab. |). 

The stratégy that has been employed most frequently in the literatuře is to 
use n as large as possible in order to suppress the contribution from the bb con- 
tinuum. Here, n can already been considered large for n > 4. The experimental 
uncertainties for large-n extractions of the bottom mass are around 15 MeV. In 
this case, theoretical uncertainties dominate and there is the chance to reduce 
the uncertainties in the bottom mass extractions by theoretical considerations. 
In large-n moments the bb dynamics is dominated by non-relativistic physics 
and it is necessary to sum the Coulomb singularities {oia/vY in the cross section 



(see Eqs. (10£) and (124)). Because the size of the energy range contributing 
to the n-th moment is of order M/n, the typical center-of-mass velocity of the 
bottom quarks in the n-th moment is of order d„ = 1/ ^/n. This counting rule 
allows for a quantitative formulation of the two requirements for the choice 
of n: theoretical reliability demands Mv^ 3> Aqcd, and dominance of the 
non-relativistic dynamics demands <C 1. If both requirements are met, the 
bb dynamics encoded in the moments is non-relativistic and perturbative, i.e. 
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thc hierarchy M ^ A/w„ ^ Mv^ ^ Aqcd is vahd and thc momcnts can be 
computed with the effective theory methods described in the previous section. 
Thus, from the conceptual point of view there is only a quite narrow window of 
n values between about 4 and 10 that should be used for large-n sum rules. If 
these considerations are apphed to cč systems, there is no window of n values 
left for which a sum rule analysis based on non-relativistic moments can be 
justified. 

At present only fixed order analyses have been carried out for the large-n 
T sum rules. Initially the bottom pole mass was determined. After the re- 
levance of the bad high order behavior of the pole mass for QQ systems was 
fuUy appreciated and suitable threshold mass schemes were proposed (Sec. 0), 
subsequent analyses extracted threshold masses, which were in a second step 
converted to the MS mass mijň). In the foUowing I give a historically ordered 
compilation of the various sum rule analyses including brief descriptions of the 
respective methods. The main results that have been given by the authors are 
collected in Tab. |2[ In Tab. ^ I have also coUected bottom quark mass results 
from hh sum rules not discussed in this review and from perturbative computa- 
tions of the T(1S) mass (Sec. 11.2). The results are compatible to recent lattice 



analyses (see Refs. [150Land references therein) and the measurements from 
three-jet rates a LEP.E^ It should be noted that the uncertainties obtained 
from recent lattice analyses have uncertainties that are comparable to the ones 
obtained from the perturbative analyses discussed in this revie^^f—tjlowever, the 
lattice numbers for mijn) generally tend towards larger values. 
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10.2 Pole Mass Results 



ExtractipiL^ of the bottom gt^arV pole mass were carried out at NI^ by 
Voloshin 



by Hoang 



liM^nd by Kiihn et al.JHZI and at NNLO by Penin and Pivovaro vN and 
r|i3^ In this section the methods and results are briefly described. The 
corre spon ding values for m(m) are computed from the formula for Tip^j^/r™ in 



Eq. (lOC). The two-loop relation is used for a NLO sum rule analysis and the 
three-loop relation for a NNLO sum rule analysis. This order-dependent con- 
version is necessary to achieve the numerical cancellation of the lar ge in frared- 
sensitive contributions contained in the pole mass definition (Sec. [7.21) . Since 
in all analyses the Information on the exact correlation of the numerical pole 
mass value with the set of theoretical parameters used for the fits was not 
quoted, I carry out the conversion treating the pole mass result as uncorre- 
lated. All uncertainties are added quadratically and perturbative conversion 
uncertainties are included. The values for the MS mass obtained from the pole 
mass in this way needs to be interpreted with some caution, because the fuU 
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author 


mi,(mi,) 


othcr mass 


comments, Rcf. 






Voloshin 95 




íTip^i^ = 4.83 ± 0.01 


NLO T sum rulcs, no 


ft^unccrt. |1 


se 




Kuhn 98 




TOp^i^ = 4.78 ± 0.04 


NLO T sum rulcs 






Ponin 98 




rrip^i^ = 4.78 ± 0.04 


NNLO T sum rulcíj 






Hoang 98 




mp„i^ = 4.88 ± 0.13 


NLO T sum rulcs 


1 


M- 




Hoang 98 


4.26 ± 0.09 * 


TOp^i^ = 4.88 ± 0.09 


NNLO T sum rulcs 


Li 




Mclnikov 98 


4.20 ± 0.10 


M^^^"" = 4.56 ± 0.06 


NNLO T sum rulcs 






Penin 98 


4.21 ± 0.11 * 


rrippi^ = 4.80 ± 0.06 


NNLO T sum rulcs | 


^ 








Jamin 98 


4.19 ± 0.06 




T sum rulcs; no cxar 




ifo 


L41 




Hoang 99 


4.20 ± 0.06 


Mjs = 4.71 ± 0.03 


NNLO T sum rulcs 


m 


i 






Bcnokc 99 


4.26 ± 0.09 


M'^'g'=^ = 4.60 ± 0.11 


NNLO T sum rulcs | 


m 


) 






Hoang 00 


4.17 ± 0.05 


Mjs = 4.69 ± 0.03 


NNLO T sum rulcs. charm ruas^ 






1 


Kuhn 01 


4.21 ± 0.05 




low-n T sum rulcs, O(a^) 






Pincda 97 




"pole = 5.00 + 07 


NNLO T[1S) mass & non-pcrt. cff. 


I54I 


Bonckc 99 


4.24 ± 0.09 


Af^^^v = 4.58 ± 0.08 


NNLO T(1S) mass & non-pcrt. cff. 




Hoang 99 


4.21 ± 0.07 


Mis = 4.73 ± 0.05 


NNLO T(1S) mass &. non-pcrt. cff. 




Pincda 01 


4.21 ± 0.09 


Aíj^^"^ = 4.39 ± 0.11 


NNLO T(1S) mass & non-pcrt. cff. 




Brambilla 01 


4.19 ± 0.03 




NNLO T(1S) mass, p.th. only | 





Table 2: CoUection in historical order of results of bottom quark mass determinations from 
T sum rules (upper part) and calculations of the T(1S) mass (lower part). AU results are 
given in units of GeV and rounded to units of 10 MeV. Only results where Os was taken as 
an input are shown. The uncertainties quoted in the respective references have been added 
quadratically. For analyses with several authors only the respective first author is quoted. 
The MS mass result that are indicated with a star have been determined with an inconsistent 

conversion formula. 



correlation should have been taken into account (Sec. |7.l| ). 

The NLO analysis by VoloshiiJHj was the first large-n sum rule anal- 
ysis consistent with non-relativistic power counting. Voloshin used time- 
independent perturbation theory to determine the NLO cor rections to the cor- 
relator Ai (Eq. (119)). The moment integration of Eq. (138) was carried out in 
a strict non-relativistic expansion. After deforming the contour into the com- 
plex energy plane away from the bound statě poles the integration can be solved 
using inverse Laplace transformations. This method has the feature that bound 
State and continuum contributions are treated at the same time and in exactly 
the same way with respect to the non-relativistic expansion, i.e. the corrections 
to bound statě energies in the energy denominators of the resonances are also 
expanded. For each moment Voloshin set the scale /i in the Coulomb potential 
such that t he N LO correction to the Green function vanishes. The scale of ag 
in Cl (Eqs. (117)) was set to the bottom mass. He carried out a two parameter 
fit for mpoic and Us using four moments P„ (n = 8, 12, 16, 20) and obtained 
TOpoic = 4.827 ±0.007 MeV and as{Mz) = 0.109 ±0.001. The result does not 
account for perturbative uncertainties. Using Voloshin's value for the cor- 



Heavy Quarkonium Dynamics 



99 



responding MS mass, and two-loop conversion, gives m{m) — 4.33 ±0.01 GeV. 
For asiMz) = 0.118 ± 0.003 one obtains m{m) = 4.27 ± 0.04 GeV. 

Voloshin's NLO analyses was repeat ed b y Kúhn, Penin and Pivovarov 
(KPP) in Ref. flŠ^] and Hoang in Ref. ||l3|]. KPP used time-independent 
perturbation theory for Ai and solved the moment integration exactly in the 
form of Eq. ( 138 ). This method treats bound statě and continuum contri- 
butions diffcrcntly, because the bound statě energy corrections in the energy 
denominators of the resonances are not expanded out. KPP also improved the 
convergence of the perturbative series by absorbing some part of the higher 
order corrections in the Coulomb potential into the LO Coulomb charge. KPP 
fitted individual moments for 10 < n < 20 taking as{Mz) — 0.118 as an in- 
put. The scale of in ci was set equal to ^, the scale in the LO Coulomb 
potential. They obtained TOpoic — 4.78 ± 0.04 GeV, where the uncertainty was 
estimated from the variation of the result with n, jj. ~ 1.2-4.8 GeV and when 
the NNLO corrections to the Coulomb potential in the Schrodinger equation 
are included. For the corresponding MS mass, using the two-loop relation, one 
obtains míro) = 4.18 ± 0.05 GeV for a^Mz) = 0.118 ± 0.003. 

Hoangll3 used Voloshin's method to determine the NLO moments, but 
did not fix /i. He carried out fits using moments for various scts of ?i's in 
the range between 4 and 10 for a given value of fi. Perturbative uncertainties 
were estimated by repeating the fit many times choosing random values for 
fi between 1.5 and 3.5 GeV and for the scale of in ci between 2.5 and 
10 GeV. In a two-parameter fit he obtained rripoic = 4.78 ± 0.14 GeV and 
as{Mz) = 0.109 ± 0.023. The result for the pole mass corresponds to m(m) — 
4.18±0.13 GeV at the two-loop level for asiMz) = 0.118±0.003 as an input. 
Taking as{Mz) as an input Hoang obtained nipoic = 4.88 ± 0.13 GeV for 
asiMz) = 0.118 ±0.004. This corresponds to m(m) = 4.28 ±0.13 GeV at the 
two-loop level for asiMz) = 0.118 ± 0.003. 

In the first NNIJD analyses by Penin and Pivovarov (PP)@ and the NNLO 
analysis by Hoangi^S time-independent perturbation theory was used to de- 
termine the NNLO corrections to the correlators. The UV-divergences were 
regularized with a cutoff and the tazo-loop corrections to ci were determined 
using the direct matching mcthodE3 The moment integration and the fitting 
proceduře were equivalent to the respective methods by P P an d Hoang in their 
NLO analyses. PP improved the convergence as in Ref. [137| and fitted indi- 
vidual moments for 10 < n < 20 taking asiMz) = 0.118 as an input. The scale 
of as in Cl and the cutoff scale were set equal to /i, the scale of as in the LO 
Coulomb potential. They obtained mpoic — 4.78 ± 0.04 GeV, where the uncer- 
tainty was estimated from the varying n and the scale /i = mpoic ± 1 GeV. In 
order to determine Wiijn) from this NNLO pole mass result one now has to use 
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thc threc-loop rclation for rí£ = 4 massless quarks. For as{Mz) = 0.118±0.003 
one obtains m(m) = 4.05 ± 0.07 GeV. 



In his NNLO analysis HoangE£3 carried out fits using moments for 
various sets of n's in the range between 4 and 10 for a given value of fj,. As 
in his NLO analysis, perturbative uncertainties were determined from varying 
independently /x between 1.5 and 3.5 GeV, and the scale in in ci and 
the eutoff between 2.5 and 10 GeV. In a two-parameter fit Hoang obtained 
TOpoie = 4.83 ± 0.07 GeV and as{Mz) = 0.110 ± 0.014. The pole mass value 
corresponds to ni{m) = 4.10 ± 0.09 GeV at the three-loop level for as{Mz) = 
0.118±0.003. Taking as{Mz) = 0.118±0.004 as an input he obtained mpoic = 
4.88±0.09 GeV, which corresponds to TO(m) = 4.14±0.10 GeV using the three- 
loop conversion for as{Mz) = 0.118 ± 0.003. For the determination of the MS 
mass Hoang used the two-loop relation and obtained mini) = 4.26 ± 0.09 GeV 
for the two-parameter fit and nT{m) — 4.24 ± 0.08 GeV, if as is taken as an 
input. In both cases Hoang accounted for the correlation between TOpoic and 
as{Mz) obtained from the fitting. 

In the second NNLO analysis by PPEa the methods from their first analysis 
were used, but was varied between 3.5 and 6.5 GeV and moments for n 
between 8 and 12 were employed. In this analysis PP obtained mpoio = 4.80 ± 
0.06 GeV, taking as (Mz) = 0.118 as an input. Fot as (M z) = 0.118±0.003 this 
result corresponds to mim) = 4.06 ±0.08 GeV if the proper three-loop relation 
is used. PP used the two-loop relation and obtained m{m) = 4.21 ± 0.11 GcV. 



10.3 Threshold Mass Results 



ExtractionSrAf threshnl 



-cnasses at NNLO were caxried out by Melnikov and 



Yelkhovsky P Hoan^^O and Beneke and SignerEžj In this section the meth- 
ods and results of the respective analyses are briefly descnbed. 

In the analysis by Melnikov and Yelkhovsky (MY)l£!l the kinetic mass 
Mj^jjj(l GeV) was determined. The theoretical moments were determined in 
time-independent perturbation theory in the pole mass scheme using the same 
methods as previous NNLO analyses. The moment integration was carried out 



in thc form Eq. (138) and the energy denominators of the bound statě reso- 
nances were not expanded out. MY fitted individual moments for n = 14, 16, 18 
taking as{Mz) — 0.118 as an input and extracted the bottom pole mass as 
a correlated parameter for a given choice of /i between 2 and 4.5 GeV. The 
scale of as in ci and the eutoff were fixed to 5 GeV. For each value of fi MY 
determined subsequently the kinetic mass using the numerical value for /i for 
which the pole mass value was determined. To achieve the numerical cancel- 
lation of the large infrared-sensitive corrections it was important to keep the 
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Scales 


Order of PT 


14 


16 


18 


^polc 


Mkin 


^polc 




^polc 


-Aíkin 


fi = 4.5 GeV 
Hd = 5 GcV 
A = 5 GcV 


LO 


4.69 


4.57 


4.7 


4.57 


4.7 


4.575 


NLG 


4.74 


4.49 


4.75 


4.50 


4.76 


4.51 


NNLO 


4.89 


4.51 


4.90 


4.52 


4.91 


4.53 


= 3.5 GcV 
/ici = 5 GeV 
A = 5 GcV 


LO 


4.73 


4.59 


4.73 


4.59 


4.73 


4.595 


NLG 


4.77 


4.50 


4.78 


4.511 


4.79 


4.52 


NNLO 


4.95 


4.54 


4.96 


4.55 


4.96 


4.55 


Ai = 2.5 GeV 
Hd = 5 GeV 
A = 5 GeV 


LO 


4.79 


4.63 


4.79 


4.63 


4.79 


4.63 


NLO 


4.82 


4.505 


4.83 


4.517 


4.84 


4.525 


NNLO 


5.08 


4.595 


5.08 


4.595 


5.09 


4.6 


/i = 2 GeV 
yUci = 5 GeV 
A = 5 GeV 


LO 


4.84 


4.67 


4.84 


4.666 


4.84 


4.66 


NLO 


4.84 


4.49 


4.85 


4.505 


4.86 


4.514 


NNLO 


5.21 


4.66 


5.21 


4.66 


5.21 


4.66 



Table 3: The bottom quark pole mass and the kinetic mass M]^i^{l GeV) in uuits of GeV 
for ^ between 2 and 4.5 GeV obtained by Melnikov and Yelkhovsky in Ref. ||97| for the 

moments n = 14, 16, 18. 

bound State energy denominators unexpanded. In Tab. || the results obtained 
by MY for between 2 and 4.5 GeV are displayed. MY found pole mass val- 
ues in the range between 4.5 and 5.2 GeV, where the NNLO shifts were larger 
than the NLO one. For the kinetic mass the results were more stable in the 
range between 4.49 and 4.67 GeV. However, the NNLO shifts were still quite 
large compared to the NLO ones. For fj, = 4.5 the convergence is quite good, 
but for /i = 2 GeV the NNLO and NLO shifts are of equal size. Since the 
resulting series for mkin(l GeV) is alternating, MY argued heuristically that 
the convergence could be improved by an Euler transformation on the series. 
In this way they obtained the result 

A4i„(l GeV) = 4.56 ± 0.06 GeV, 

which is somewhat between the NLO and NNLO values and has an uncertainty 
that is smaller than the range covered by the NLO and NNLO values obtained 
in the analysis. Since the kinetic mass is a short-distance mass there is no 
conceptual need to convert the result to the MS mass at a particular order. 
In generál, one should go for the best (i.e. highest order) conversion. MY 
determincd the MS mass at two-loop order and obtained m(m) = 4.20 ± 
0.10 GeV. 
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In the analysis by Hoangi^j the IS mass was determined. The theo- 



retical moments were computed directly in the IS mass scheme using time- 



independent perturbation theory and the methods of Ref. |138| ] (Sec. 10.2 ). 
All NLO and NNLO corrections in the Schrodinger equation in Eq. (110) 
were strictly treated perturbatively, i.e. the cancellation of the large infrared- 
sensitive corrections associated with the pole mass were eliminated analyti- 
cally. For the fitting proceduře ex actly the same method was used as for the 
earher pole mass analysis in Ref. 138|. Hoang carried out fits using var- 



ious sets of n's in the range between 4 and 10. For all moments in a given 
set the same value of /i was chosen. This fitting proceduře is different from 
fits of individual moments, because it put highest statistical weight on lin- 
ear combinations of moments that are the least sensitive to the experimental 
uncertainties in the electronic partial widths of the T mesons. Since the corre- 
lation of the experimental moments coming from the e"'"e~ partial widths (see 



Eq. (140)) is similar to the correlation of the theoretical moments coming from 
their dependence on the bb wave function at the origin, the resulting func- 
tion has a smaller sensitivity to variations of fi than the individual moments. 
In contrast to fits of individual moments, which exclusively use Information 
on the size of the moments, the fitting proceduře also uses Information 
on the relative size of the moments as a function of n. Hoang estimated the 
perturbative uncertainties by repeating the fit many times choosing indepen- 
dently random values for fi between 1.5 and 3.5 GcV, for the scale of in 
Cl and for the cutoff between 2.5 and 10 GeV. In a two-parameter fit Hoang 
obtained Aí^g = 4.67 ± 0.07 GeV and as(Afz) = 0.112 ± 0.023 at NLO and 
Mis = 4.70 ± 0.03 GeV and asiMz) = 0.116 ± 0.013 at NNLO. A graphical 
representation of the NLO and NNLO fit results is displayed in Fig. ^a and 
b. In a single parameter fit taking as{Mz) — 0.118 ± 0.004 as an input Hoang 
obtained Aí^g = 4.70 ± 0.05 GeV at NLO. At NNLO he obtained 

Mig = 4.71 ± 0.03 GeV, 

which was considcred the main result of the analysis. A graphical represen- 
tation of the NLO and NNLO fit results is displayed in Fig. |3^ and d. The 
results in Fig. ^ show that the IS mass has a very weak correlation to a^. 
Thus the uncertainty in as is practically irrelevant for th£ error in the IS mass. 
In contrast_tD the analysis by Melnikov and YelkovskyÉZi and later by Beneke 



and SignerEfj the results at NLO and NNLO by Hoang are consistent with 
each other. The range of the NNLO result is contained entirely in the range of 
the NLO result. This is a consequence of the fitting method used by Hoang. 
Hoang converted the result to the MS mass at two loops using the large-/3o 
approximation for the three-loop correction as an estimate for the perturbative 
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FÍOTír n .'^2 : Result for the allowed region in the M^^-Os plane obtained by Hoang in 
Ref. |139|. (a) Two parameter fit at NNLO and (b) NLO. The dots represent points of 
minimal for a large number of random choices for the scales and sets of ?i's. Experimen- 
tal errors are not displayed. (c) Single parameter fit at NNLO and (d) NLO taking as as an 
input. The dots represent points of minimal for ^ large number of random choices for the 
scales, sets of ti's and . Experimcntal errors at 95% CL are displayed as vertical bars. 



uncertainty. Hc obtained TO(m) 4.20±0.06 GeV for as{Mz) = 0.118±0.004. 
The uncertainty is larger than for the IS mass because of the perturbative un- 
certainty in the conversion (40 MeV) and the uncertainty in as{Mz) (40 MeV), 
which arises from a relatively large one-loop correction in the relation between 



IS and MS mass, see Eq. (100). In Ref. [[73[ Hoang repeated the analysis 
taking into account also the effects of the non-zero charm quark mass. For 
řňcijřic) = 1-4 ± 0.3 GeV he obtained 

Mis = 4-69 ± 0.03 GeV 

for the IS mass for the same fitting proceduře. The result was converted to the 
MS mass using the fuU three-loop relation leading to m{řří) = 4.17 ±0.05 GeV 
for as{Mz) = 0.118 ±0.003. 



In the analysis by Beneke and Signer (BS)lil3 the PS mass íTipg(2 GeV) was 



determined. The theoretical moments were computed in the PS mass scheme 
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using time-independent perturbation theory to solve the Schrodinger equation 
of Eq. ( |110| ). BS used dimensional regularization and identified the form of 
the potentials and coefhcients using the threshold expansion, i.e. they did not 
carry out an exphcit matching calculation. The moments of BS depend on /lí, 
the scale of as and the factorization scale, which was set to the bottom quark 



mass. The moment integration was solved in the form of Eq. (138), but all 
NLO and NNLO corrections were strictly treated perturbatively by expanding 
out the resonance delta functions in i?" in order to achieve an exphcit ana- 
lytical cancellation of the large infrared-sensitive corrections associated with 
the pole mass definition. BS also carried out a summation of some NLL order 
logarithms in ci. The prescription they used, however, is not consistcnt with 
the known NLL order anomalous dimensions in pNRQCD or vNRQCD. BS 
fitted the n = 10 moment taking as{Mz) — 0.118 ± 0.003 as an input and 
varying ^ between 2 GcV and 2mpg(2 GeV)/\/TO. At NNLO they obtained 

Mps(2 GeV) = 4.60 ±0.11 GeV , 

where the scale variation is the dominant source of uncertainty (100 MeV). The 
uncertainty from the strong coupling (30 MeV) and from the experimental 
data (20 MeV) are smaU. At NLO BS obtained mps(2 GeV) = 4.44 GeV 
as the centrál value with a smaller scale variation. No explicit uncertainties 
for the NLO result were quoted. A graphical representation of the results is 
shown in Fig. Similar to the analysis of Melnikov and Yelkhovsky, which 
was also based on fits of individual moments, the NLO and NNLO results 
are not compatible numerically. BS considered the NNLO result given above 
as the main result of the analysis. The discrepancy with the NLO number 
was not included in the error estimate. The PS mass at 2 GeV has a stronger 
correlation to the input value oí as{Mz) than for example the IS mass. Taking 
into account this correlation to as and fj, BS obtained m(m) = 4.26 ±0.09 GeV 
for the MS mass using two-loop conversion. The mean value also contains a 
— 10 MeV shift from a large- /3o estimate for the (at the time of the analysis 
unknown) three-loop correction. The uncertainty from the strong coupling 
(10 MeV) and from the conversion (20 MeV) were small. The mean value for 
the MS mass obtained by BS is highcr than the one by MY. This is because 
BS did not take into account their NLO results in the determination of the 
mean value. 

The analyses by MY, Hoang and BS are conceptually equivalent with 
respect to the computations of the moments in the framework of NRQCD. The 
differences in the results are a consequence of differences in the fitting methods 
and the treatment of perturbative uncertainties. Because the uncertainties are 
dominated by theory, the error estimate is, naturally, a quite delicate task. As 
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Ml 



PS 




Figuře 3i The value of Mps(2GeV) in units of GeV obtained by Beneke and Signer 

in Ref. 1 14C ] from the lOth moment as a function of the renormalization scale in NLO 
and NNLO and for Qs{Mz) = 0.118. The dark region specifies the variation due to the 
experimental error on the moment. The two outer lineš display the scale variation from 
which the perturbative error was computed. 



such none of the analyses can a priori be considered better or more reahstic 
than another one. Nevertheless, the quoted uncertainties have been subject to 
controversial discussions, which shaU be not repeated in this review. From a 
(at least for me) conservative point of view it is fair to say all three analyses 
contain elements that aUow to consider the quoted uncertainties as optimistic. 
The analysis by MY and BS showed that the resuhs for the threshold masses at 
NLO and NNLO obtained from fits of individual moments are not compatible. 
MY deah with this problém by using heuristic arguments on alternating series 
and BS by ignoring the NLO results in the finál result. The origin of the 
difference between NLO and NNLO results is that the R-ratio that is contained 
in the moments of Eq. (138) has quite large NNLO corrections. We face this 
problém also in fixed order computations of the tt production cross section 
close to threshold, which is discussed in Sec. 9.3, see Fig. In the fitting 
proceduře by Hoang the normalization problém of the R-ratio was treated by 
also using Information of the relative size of the moments for different n, similar 
to taking ratios of the moments. 

Clearly, the situation is not completely satisfactory and further improve- 
ments in the computation of theoretical moments will be necessary to settle the 
question of the uncertainties in the bottom quark determination from T sum 
rules. One promising path might be to use renormalization- group-improved 
perturbation theory for the computation of the moments. In the case of top 
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production this has indccd improved the stability of the R-ratio (Sec. p^ )- At 
the present time a renormalization-group-improved sum rule analysis has not 
yet been carried out. 



11 Heavy Quarkonium Spectrum 

Theoretical dcscriptions of the heavy quarkonium spectrum represent the clas- 
sic application of non-relativistic quarkonium physics. Early work in QCD- 
inspired potential models provided very precise determinations of charmonium 
and bottomonium Sftectra and physical insights into the long-distance aspects 
of interquark forcesB However, potential models cannot be derived quantita- 
tively from first principles in QCD, which makes them less suitable for deter- 
minations of QCD parameters, most notably the heavy quark masses. With 
the advent of effective theories for heavy quarkonium systems (Secs. ||, |[ |]) 
a first principle description of the spectrum has become possible conceptually, 
but precise computations are difRcult, because a reliable technology only ex- 
ists for systems where the hierarchy m 3> mv 3> mv'^ ^ Aqcd is valid. This 
is the case that is discussed generically in Sec. |ll.l[ . In this casc the predic- 
tions of the effective theories can be computed perturbatively supplemented by 
local condensate contributions, which incorporate averaged non-perturbative 
effects. As a matter of principle, at most the bottomonium ground states can 
be expected to be described with this technology. Determinations of the bot- 



tom quark mass from the mass of the T(1S) statě are described in Sec. 11.2 
Considerations of the toponium spectrum are irrelevant, due to the large top 
decay width Tt « 1.5 GeV > Aqcd (see Sec. |). 

For most quarkonium systems mv'^ is of order Aqcd or even smaller. In 
the franji£.work o f pN RQCD an approach for these systems has been roughly 
outlinedllj (Sec. [4.6|) , but no systematic prescription to carry out quantita- 
tive computations of the spectrum has yet been devised. In the framework of 
vNRQCD such systems can a priori not be treated because ultrasoft modes 
are by construction already present at the hard matching scale and because 
soft and ultrasoft scales are correlated at all times (Sec. ^). In the framework 
of NRQCD (Sec. on the other hand, a systematic treatment of the spec- 
trum of these systems is unknown with analytic methods, but is possible with 
lattice simulationsLj On the lattice, theoretical precision is limited by system- 
atic uncertainties, such as from unquenching and extrapolating and matching 
to the continuum limit. No further discussion on lattice methods for heavy 
quarkonium spectra is given in this review. In a recent analysis Brambilla 
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et al.li£3'LŽJ exaniincd the charmonium, bottomonium and Bc spcctruni using 
only the perturbative contributions for the case m 3> mv 3> mv'^ 3> Aqcd- 
Setting the renormahzation scale using a minimal sensitivity prescription 
they found remarkable agreement of the bottomonium sprectrum up to the 
n = 3 radiál excitation without accounting for any non-perturbative eíTects. 
Brambilla et al. concluded that perturbative methods devised for the case 
m 3> mv 3> mv"^ 3> Aqcd might have a wider range of apphcabihty than 
previously expected. For subsequent discussions foUowing the same interpre- 
tation see also Ref. |l33|| . The results of BrambiUa et al., however, appear to 
be in sharp contradiction to results on the size of non-perturbative corrections 
obtained for higher bottomonium excitations in Ref. |86|. 



11.1 Calculation of the QQ Spectrum 

For the case m 3> rav 3> mv^ 3> Aqcd the computation of the heavy quarko- 
nium spectrum foUows the presentation given in Sec. ^ on heavy quark pro- 
duction close to threshold. At NNLL order in the renormalization-group- 
improved approach of pNRQCD (Sec. |) or vNRQCD (Sec. |), or at NNLO 
in the fixed order NRQCD (Sec. ||) approach the perturbative contribution of 
QQ dynamics can be described by a common text-book two-body Schrodinger 
equation. In momentům space representation and in a threshold mass scheme 
the Schrodinger equation has the generic form 



*(P) 



^F(p,k)vl/(k)= 0. (142) 

The perturbative contribution of the masses of the QQ states are given by the 
energy eigenvalues of Eq. ( |142[ ) , 

iV/P;^'* = 2M + E. (143) 

The classification of the solutions according to quantum numbers is standard 
text-book knowledge and shall not be repeated here. The renormalization- 
group-improved computation for the energy eigenvalues at NNLL order takés 
the generic form 



E = Mal^{as\i^asý |l(LL);a„ (NLL); (NNLL)| 



(144) 
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In the fixed order approach logarithms of as are not summed and the energy 
eigenvalues at NNLO have the generic form 

E = Mal |l (LO); a^, (NLO); al (NNLO)| . (145) 



In Eq. (142) the difference between renormalization-group-improved and fixed 
order computation is just the form of the coefficients of the potential V (see 
Secs. 8^ a nd^.2| ). From the technical point of view the computations needed to 



solve Eq. (142) are equivalent in both approaches. In the fixed order approach 
the first compiitation of the NNL O Q Q spectrum was carried out by Pineda 
and YnduraintžJ (see also Refs. [|l55| , ^ ). Renormahzation-group-improved 



computations at NNLL order were carried out by Hoang et al. in Ref. |55| and 
by Pineda in Ref. ||3|] . Technically, the energy eigenvalues are determined us- 
ing Rayleigh-Schrodinger time-independent perturbation theory starting from 
the well known leading order non-relativistic Coulomb solution and includ- 
ing successively higher order terms from the potential and the kinetic energy 
term. To achieve an explicit cancellation of the unphysical large-order behav- 
ior associated with the pole mass it is mandatory to consistently include the 



contributions of 5M at each order of perturbation theory (Sec. 7.2). Since at 



NNLL order (or NNLO) the corresponding integrals do not contain any UV 



divergences, no regularization prescription for Eq. (142) necds to be specificd 
at this order. 

At N'^LL order (or N^^LO) the computations of the perturbative part of 
the QQ spectrum becomes more complicated. Apart from N'^LL order (or 
N'^LO) corrections to the potential V , which include for example two-loop 
corrections to the l/TO|k| potential or three-loop corrections to the Coulomb 
potential and the corresponding three- and four-loop anomalous dimensions, 
also non-logarithmic corrections from ultrasoft gluons need to be considered 
(see for example Figs. ^ and [l8|) . The latter corrections involve UV divergences 
and depend on the regularization scheme. The corresponding UV divergences 
are canceled by conterterms of the l/m?, l/TO|k| and potentials (Secs. 



4.3 



5.7 and |6|), which become renormalized. The sum of the ultrasoft cor- 
rections and the corrections from the potentials are scheme- independent. The 
ultrasoft corrections are the QCD analogue of the Bethe-log corrections known 
from QED and cannot be represented by a two-body Schrodinger equation, 
because they correspond to a higher Fock QQ-Gluon statě. Some corrections 
at N'^LO in the fixed order approach are already known. The non-logarithmic 



corrections from ultrasoft gluons were determined in Ref. |11£], and two-loop 
corrections to the l/TO|k| potential were determined in Ref. ||34|. The N^^LO 
corrections in the large-/?o approximation were determined in Refs. 
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From thc NNLL order rcnormalization-group-improvcd computation (see 
also Rcf. pl]) all N^LO corrections proportional to powers of Ina^ are also 



known. In Ref. |15É] Kiyo and Sumino claimed that additional NNLO cor- 
rections oc Ma^ and Ma^ In ag are found when Bethe-Salpeter rather than 
effective theory techniques are employed with dimensional regularization. The 
contributions found in Ref. [|l56| originate from an incomplete expansion of 
the small fco component of the one-loop corrections to the Coulomb potential 
in the potential region (Eq. (^l])) and are unphysical. 

In the case m ^ mv 3> mv'^ ^ Aqcd the non-perturbative contributions 
to the spectrum can be exnceased as an expansion in terms of local ghion 
and light quark condensatesQŮ At leading order in the multipole expansion 
the first non-perturbative correction arises from radiation and absorption of 
a gluon with cncrgy of order Aqcd from the QQ pair. Conceptually this 
contribution is dosely related to the ultrasoft correction mentioned before. 
For mu^ ^ Aqcd, both contributions can be separated by an expansion in 
Aqcd/™^^, see Refs. [pp^. Physically, this expansion means that the time 
span between emission and absorption of the gluon is much smaller than the 
non-perturbative correlation time. In configuration space representation the 
coupling of the gluon to the QQ pair is just the r.E interaction in pNRQCD. 
In momentům space representation the interaction arises from the diagrams in 



The first non-perturbative correction is in analogy to Eq. (132) and 



where ais — 1.652, a2s = 1.783, 035 = 1.869, a2P — 1.123, etc. The next 
order terms in the Aqcd jmv^ expansion correspond to dimension-six conden- 
sates and have been determined by Pineda in Ref. [ 157 1 . It should be noted that 
the coefficients of the operátor expansion in local condensatcs are at present 
only known at LO (or LL order) in thc non-relativistic expansion and still have 
a rather large scheme-dependence. 

The summation of all orders in K QQXi/m v^ leads formally to non-local con- 
densate contributions (see e.g. Refs. [EoLl5q|) describing the non-perturbative 
corrections for Aqcd ~ mv^. A reliable quantitative treatment of non- 
perturbative corrections for this case in quarkonium systems has not yet been 
developed. However, for Aqcd ^ also the perturbative contributions are 
affected because gluons with ultrasoft energies and energies of order Aqcd are 
identical and cannot be separated, and because as{mv^) has to be considered 
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of order one. This means that thc non-perturbativc/ultrasoft corrections have 
to be counted as NNLO contributions. 



11.2 Bottom Mass Determinations from the T{1S) Mass 



There are a number of analyses wherc thc bottom mass was determined from 
the T(1S) mass using fixed order calculations at NNLO as described in the 
previous section. From the technical point of view, in aU analyses it was 
assumed that the hierarchy m 3> mv 3> mv^ 3> Aqcd is vahd and that 
non-perturbative corrections can be expressed as an operátor expansion in 
terms of local condensates. The expression for the perturbative contribution 
of the T(1S) mass at NNLO in a threshold mass scheme (Sec. 7.2) with 5M — 
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where, assuming massless Hght quarks, 
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The one- and two-loop coefRcients of the beta-function, /3o and f3i, and the 
constants oi and 02 are from the corrections to the Coulomb potential in 
Eq. In the IS mass scheme <5A/lo,nlo = -Af A^^-^^^O and SM^^^^^ = 

—MA™^'-' — Af (A'^^)^, and, by construction, p ertur bative corrections are 
zero. The non-perturbative contribution from Eq. (146) reads 



pnll.p. _ 

^10 ~ 



0.465 TT 

Af 3 ai 



(149) 
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Typical values for thc giuon condensate used in the literatuře are in the range 
(OlasG^M^^lO) = 0.05 ± 0.03 GeV''. The thpa|iold mass is extracted from fit- 



ting AfP^fgj +-E"oP- to the experimental valuelll3 M^^^^g^ = 9460.30±0.26 MeV. 

In the foUowing I review the bottom mass determinations bascd on fixed 
order NNLO calculations in historical order. The results obtained in the anal- 
yses are cohected in Tab. ^j. To simpHfy the presentation on the detcrmination 
of uncertainties I discuss only the dominant source of uncertainties as obtained 
in the different analyses. 



Pineda and Yndurain (PY)ll£j determined the pole mass. For the strong 
coupling PY used a scheme where part of the corrections to the Coulomb po- 
tential are absorbed into a^. In this scheme the size of the non-perturbative 
contribution was at the level of 10 MeV. PY obtained m^^i^ = S.OOÍg GeV, 



where the uncertainty was dominated by the error in the strong coupling, 
as{Mz) = 0.114+o;°o4. PY also determined the MŠ mass using two-loop 
conversion, mim) = 4.44lg'Q3 GeV. Converting to the MS mass at three 
loops, which is necessary to achieve the numerical cancellation of the large 
unphysical corrections associated with the pole mass, the centrál value reads 
nT{m) = 4.32 GeV for as{Mz) = 0.114, using the formula for ra^^^^/rm in 
Eq. (^). 

HoangE^ determined the IS mass. In the IS scheme there is, by construc- 
tion, no correlation to a^. Perturbative uncertainties cannot be estimated in 
the usual way by changing the renormalization scale ^. A conservative estimate 
would be, for example, to account for the difference in the result when Annlq is 
replaced by its large-/3o approximation. Hoang assigned 100 MeV for the size of 
the non-perturbative contributions and treated them as an overall uncertainty 
without including them in the detcrmination of the centrál value. He obtained 
Mj^g = 4.73 ± 0.05 GeV. Perturbative uncertainties were not included. Using 
three-loop conversion to the MS mass he obtained m{m) = 4.21 ±0.07 GeV for 
as{Mz) — 0.118±0.004. The uncertainty is larger than the one of the IS mass 
because the MS -IS mass relation has a rather strong dependence on as com- 
ing from the one-loop correction in the MS -pole mass relation. Since the IS 
mass (as well as the PS and the RS mass mentioned below) are short-distance 
masses the conversion to the MS mass does not have to be carried out at a 
specific order. However, one should go for the highest order conversion to keep 
perturbative uncertainties small. A compact approximation of the three-loop 
conversion formula for the MS -IS mass relation, which also accounts for charm 
mass effects and the scale /i in the strong coupling, was given in Ref. [ |73[ |, 

m{m) = 4.169 GeV - 0.01 {m^imc) - 1.4 Gev) + 0.925 (m^ - 4.69 Gev) 
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-9.1 (a(5HA//2)-0.118) GeV + 0.0057 (/i- 4.69) GcV . (150) 

For nTc{nic) > 0.4 GeV and fj, > 2.5 GeV the difFerence between this approxi- 
mation formula and the exacL-cesult is less than 3 MeV. 

Beneke and Signer (BS)Ež3 determined the PS mass. They estimated 
the uncertainty in the PS mass coming from the non-perturbative contri- 
butions as ±70 MeV and obtained Mps(2 GeV) = 4.58 ± 0.08 GeV, where 
non-perturbative contributions were not included in the centrál value and 
as{Mz) — 0.118 ± 0.004 was taken as an input. For the MS mass BS used 
three-loop mnversion, and they obtained Wí(jřT) = 4.24 ± 0.09 GeV. 

PinedaE3 determined the RS and RS' masses assuming a 75 MeV uncer- 
tainty from the non-perturbative and N^^LO ultrasoft contributions. He ob- 
tained Mp,g(2 GeV) = 4.39 ± 0.11 GeV and Mp,g,(2 GeV) = 4.78 ± 0.08 GeV. 
The error in the RS mass is significantly larger than the error from the non- 
perturbative and N^^LO ultrasoft contributions due to a theoretical ambigu- 
ity in the subtraction term that defines the RS mass. For the MS mass 
he used three-loop conversion, and obtained a combined result of m(m) = 
4.21 ± 0.09 GeV. |_. 

Brambilla, Sumino and Vairo (BSV)ll^ determined the MS mass directly 
from the T(1S) mass without taking into account non-perturbative contribu- 



tions, which technically means that they used Eq. (150) inserting half the T(1S) 
mass as the IS mass value. They fixed the renormalization scale n using a min- 
imal sensitivity requirement, which corresponds to /i « 2.5 GeV in Eq. ( |150| ). 
BSV obtained m{m) = 4.19 ± 0.03 GeV taking as{Mz) = 0.118 ± 0.002 as an 
input. 



12 Conclusions 

In the last few years the conceptual understanding of the perturbative treat- 
ment of heavy quarkonium systems with m 3> mv 3> mv'^ 3> Aqcd, w and v 
being the heavy quark mass and velocity, respectively, has improved consider- 
ably. New effective theories have been developed that are more sophisticated 
than NRQCD and that allow systematic QCD calculations using fixed order or 
renormalization-group-improved perturbation theory. Starting from solutions 
of the non-relativistic Schrodinger equation, perturbative corrections can now 
be computed in dimensional regularization to any order using a fixed set of 
rules. At the same time, the required technology has been developed to carry 
out computations at next-to-next-to-leading and higher orders, and the role 



Heavy Quarkonium Dynamics 



113 



of the infrared sensitivity of quark mass definitions in the non-relativistic con- 

text has been undcrstood. In this rcvicw I havc tricd to give a comprehensive 
overview of the developments and the apphcations of the new technologies to 
bottom quark mass determinations and top quark threshold physics. What 
comes next? 

The main field of apphcation of the new technologies is top quark pair pro- 
duction close to threshold, which is a major part of the top physics program 
at the Linear Collider. This is because the top mass and the top quark widths 
are so much larger than the typical hardonization scale. Practically all aspects 
of the non-relativistic top-antitop dynamics can be described predominantly 
with perturbativc methods based on the new developments. Thus. top thresh- 
old measurements will very likely provide the insights into the interplay of 
perturbativc and non-perturbative aspects of QCD one hoped for already after 
the discovcry of J/ijj. Still, a lot of work rcmains to bc dono to cstablish the 
theoretical tools that will make this possible. At present, the new effective the- 
ory methods have exclusively been applied to computations of the total cross 
soction, whcrc non-pcrturbativc and clcctrowcak cfFccts are small. The most 
important outcome of these calculations is that top mass measurements from 
a threshold scan will have a precision at the level of 100 MeV, which will be 
unsurpassed by any othcr mcthod. But calculations should also be attempted 
with the same theoretical precision for differential observables. At certain kine- 
matic endpoints the impact of non-perturbative and also electroweak effects 
can bc enhanccd. At present, systematic studios of such effects using the new 
technologies have not yet been carried out. Measured by the progress that 
already happened in the last few years I am optimistic that developments in 
this direction will take pláce in short time. 
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